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SUMMARY 


This  report  represents  work  done  to  support  electro-optical  programs  at  the 
Dahlgren  Laboratory  of  the  Naval  Surface  Weapons  Center.  It  presents  analytical 
methods  which  were  developed  to  predict  ship  signatures  which  include  weather 
effects.  These  are  given  as  ship  contrast  temperatures  that  are  used  in  computer 
programs  to  predict  the  performance  of  present  optical  systems  used  for  infrared 
detection  and  tracking  (8-12  m  waveband).  This  report  also  presents  contrast 
temperature  predictions  for  weather  data  taken  by  ship  J  during  the  years  1964-71. 
Future  efforts  will  center  on  signature  predictions  for  other  locations  and  on  the 
expansion  of  the  analytical  models  to  include  other  ship  targets  and  background 
conditions.  This  work  was  sponsored  by  Code  N54  of  the  Naval  Surface  Weapons 
Center . 
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I.  INTRODUCTION 


The  Navy  has  long  been  Interested  In  the  possibility  of  using  infrared  sensors 
to  detect  and  track  ships  at  sea.  furthermore,  the  threat  of  anti-ship  missiles 
launched  from  ships  at  sea  has  heightened  the  Navy  interest  in  passive  surveillance 
methods.  The  continuing  development  of  detectors,  optical  materials,  and  scanning 
devices  has  resulted  in  electro-optical  systems  for  gunnery  improvement  and  passive 
detection  of  targets.  These  early  systems  established  the  feasibility  of  infrared 
detection  and  lead  to  Interest  in  a  future  program  to  develop  a  modular  Electro- 
optical  Fire  Control  Subsystem  for  present  and  future  naval  guns.  This  would  seek 
to  provide  passive  detection  and  tracking  along  with  laser  illumination  for  range¬ 
finding  for  targets  primarily  encountered  by  ships  at  sea.  The  sensor  suite  for 
this  system  includes  both  a  television  sensor  and  a  thermal  imaging  (FLIR)  sensor 
which  will  operate  in  the  8-12  micrometer  region. 


The  infrared  signatures  of  surface  ships  depend  upon  their  construction, 
operation,  and  existing  weather  conditions.  The  location  of  power  plants,  exhaust 
stacks  and  areas  of  different  construction  causes  spacial  variations  of  signature 
(emitted  radiation)  over  the  target.  However,  to  simplify  ongoing  performance 
studies,  it  was  decided  to  approximate  the  target  signature  by  an  average  signature 
which  is  constant  over  the  target  area.  This  average  signature  is  given  as  the 
ship  contrast  temperature  which  is  the  effective  temperature  difft  *ence  between  the 
ship  and  its  infrared  background.  The  ship  contrast  temperatures  must  include 
weather  effects  for  system  performance  studies.  In  these,  the  statistical  varia¬ 
tion  of  ship  contrast  temperature  with  weather  conditions  provided  by  weather  ship 
data  is  first,  found.  These  contrast  temperatures  are  then  used  in  the  system  per¬ 
formance  study  to  determine  the  percentage  of  time  that  the  ship  target  will  be 
detected. 


The  importance  of  infrared  ship  signatures  has  led  to  the  recent  development  of 
computer  codes-,'  for  their  calculation.  These  codes  divide  a  ship  into  a  large 
number  of  sections  and  numerically  compute  the  temperature  of  each  section.  The 
ship  temperatures  are  used  to  calculate  ship  radiances.  A  separate  calculation  re¬ 
quiring  atmospheric  temperatures  is  made  to  determine  the  corresponding  background 
radiance.  From  these,  contrast  radiance  or  signature  is  defined  as  the  sum  of  the 
radiance  contribution  from  each  ship  section  minus  the  background  radiance.  This 


’Batley,  P.  E.,  "Ship  Infrared  Signatures  (SIRS)  Computer  Model -Technical  Overview 
and  User's  Manual,'  NAVSHIPRANDCEN  Rept.  SME-78-37  (1978) 

;"Fleet  Signature  Computer  Model  Program  Manual,"  Westinghouse  Defense  and  Space 
Center,  Baltimore,  Maryland,  Rept.  3796A,  November  1958. 
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approach  results  in  large,  rather  slow  running  computer  programs  and  it  would  be 
costly  to  produce  signatures  for  more  than  a  few  days  of  at-sea  weather  conditions. 
The  objective  of  the  present  work  is  to  compute  ship  signatures  for  a  large  number 
of  weather  observations  so  that  the  statistical  variation  of  signature  Including 
weather  effects  can  be  given. 

A  fast-running  approximate  method  of  ship  signature  computation  based  on  the 
full-scale  codes  was  developed  to  utilize  the  large  weather  ship  data  base. 
Basically,  this  method  approximates  a  ship  by  replacing  it  with  a  single  element 
and  applies  3  correction  factors  to  average  the  spacial  variations.  The  correction 
factors  are  developed  for:  (1)  differences  in  internal  temperatures  found  through¬ 
out  the  ship;  (2)  differences  in  thermal  capacity  of  plating  on  different  sections 
of  the  ship;  and  (3)  the  presence  of  a  visible  stack  for  exhausting  products  of 
combustion.  The  technique  of  modeling  a  ship  by  a  single  element  and  3  correction 
factors  is  termed  the  Single  Element  Method.  It  is  designed  to  be  used  with 
weuther  ship  meteorological  data.  Weather  ship  data,  for  ships  operating  in  both 
the  Atlantic  and  Pacific  Oceans  were  obtained  from  the  Navy  Weather  Service.  These 
ships  provide  both  surface  and  upper  air  data  which  was  taken  either  hourly  or 
every  3  hours  in  well-defined  locations.  This  data  was  used  to  make  hourly  calcu¬ 
lations  of  ship  contrast  temperatures  for  use  in  the  statistical  compilations.  The 
Single  Element  Method  of  ship  signature  calculation  reduces  computer  running  time 
to  the  point  where  several  years  of  hourly  weather  data  cun  be  economically  run  and 
statistically  analyzed. 

The  purpose  of  this  report  is  to  describe  the  Single  Element  Method,  its 
associated  computer  program,  and  its  application  in  computing  ship  contrast  temper¬ 
atures.  Ship  signatures  were  calculated  for  a  model  of  the  U.S.  Navy  Patrol 
Frigate  located  at  one  North  Atlantic  weather  ship  location.  These  will  be  pre¬ 
sented  for  8  years  of  hourly  meteorological  observations.  Statistical  results  are 
given  for  ship,  sky,  background  and  ship  contrast  temperatures.  Ship  contrast 
temperatures  are  specified  for  day,  night  and  overall  operation  and  for  midday 
hours  where  they  are  the  largest.  Additional  computer  runs  were  made  to  determine 
the  effect  of  target  ship  heading  on  its  signature.  The  methods  and  results 
presented  herein  were  developed  to  aid  in  present  systems  performance  studies. 
However,  they  can  be  used  in  providing  estimates  of  ship  contrast  temperatures  for 
other  purposes. 
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BACKGROUND 


This  section  describes  the  methods  previously  used  to  computer  ship  signatures 
from  which  much  of  the  present  work  will  be  derived.  !t  will  cover  signature 
models  and  definitions  as  an  aid  to  understanding  the  approximate  method  described 
in  this  report.  However,  this  section  can  be  skipped  b  those  readers  mainly 
interested  in  the  statistical  results  or  tnose  already  familiar  with  the  previous 
computer  codes.  The  present  work  will  use  the  basic  models  of  the-e  codes  with 
the  important  difference  that  signatures  will  be  given  as  a  ship  contrast  tempera¬ 
ture  instead  of  as  the  radiative  power  induced  in  a  remote  sensor.  This  section 
will  briefly  describe  these  models  and  give  the  signature  related  definitions  while 
also  pointing  out  which  parts  of  the  model  will  be  of  future  use.  The  infrared 
Signature  of  a  ship  target  was  originally  defined3  as  the  signal  power  it  induces 
in  a  remote  sensor  as  the  sensor's  Instantaneous  fleld-of-view  scans  across  the 
ship  and  its  background.  This  definition  is  expressed  mathematically  by  the 
following  equation  for  radiative  signa1  power. 


P  *  ff  (Nt-  Nb) 


d\  •  da 


Here,  the  signature  is  really  the  change  i r.  irradiance  between  instantaneous  fields- 
of-view  containing  the  ship  superimposed  on  its  background,  Nt,  and  its  background, 
Nb,  respectively.  The  more  general  case  of  the  first  field-of-view  containing  the 
ship  plus  portions  of  the  background  is  not  covered  by  equation  (1).  However,  this 
case  will  not  be  considered  in  the  present  work  which  defines  the  signature  as  a 
contrast  temperature  occurring  at  the  ship  and  leaves  it  to  the  user  to  relate  this 
to  his  detector.  Furthmore,  the  detector  normalized  response  and  the  transmission 
of  the  air  path  between  target  and  detector  are  not  needed  in  the  present  analysis 
of  target  signature.  The  difference  between  the  target  surface  radiance,  Nt  and 
the  sea  surface  or  background  radiance,  Nb,  defines  the  spectral  radiance  contrast 
at  the  target.  This  quantity  will  be  later  related  to  the  contrast  temperature 
definition  of  ship  signature. 

The  target  surface  radiance  is  the  radiant  energy  that  the  target  emits  due  to 
its  surface  temperature  plus  what  it  reflects  from  its  background.  This  quantity 
is  given  by  equation  (2)  as  follows: 

Nt  =  R$t  $  ♦  R$  Af-  (1- cs)  +  Ng(l-rs)  (2) 

where  all  radiative  quantities  are  a  function  of  wavelength.  The  first  term  in 
equation  (2)  is  the  radiance  emitted  by  the  shiD's  surface  which  is  assumed  to  be  a 


3 Ibid ,  p.  2-2 
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graybody.  That  is,  this  radiance  is  the  product  of  blackbody  radiance  at  the 
ship's  surface  temperature  and  the  emissivity  of  Its  surface  coating  which  is 
assumed  to  be  constant  over  the  waveband  of  interest.  The  spectral  blackbody 
radiance  is  given  by  the  familiar  Planck's  law  as  follows: 

3.74  ^  10*  Watts 

R  - - -  (3) 

~\Mexp{l  .439  *  10V‘T-1)  nr -Sr  -  um 

The  middle  term  accounts  for  the  reflection  of  the  sun  off  of  the  target  and 
towards  the  detector.  This  quantity  is  governed  by  the  angle  factor,  Af,  which 
has  been  empirically  determined  to  be  about  10“*.  The  sun  reflection  was  found  to 
be  only  about  one  percent  of  the  remaining  terms  ir  the  8-1 2um  waveband  hence  It 
was  omitted  in  the  approximate  method.  In  the  final  term,  the  apparent  background 
radiance,  Nq,  is  the  radiance  the  target  receives  and  subsequently  reflects  (a  por¬ 
tion)  from  its  background.  The  background  seen  by  a  target  element  can  be  the  sky 
(upper  atmosphere) ,  the  air  (sea  level  atmosphere),  or  the  sea  depending  on  the 
orientation  of  the  element.  In  the  Single  Element  Method,  the  element  is  taken  to 
be  vertical  and  its  apparent  background  is  radiance  received  from  a  long  horizontal 
air  path.  This  background  is  represented  as  a  blackbody  at  sea  level  air  tempera¬ 
ture. 

Sensor  background  radiance  is  the  radiant  energy  felt  by  a  sensor  if  the 
target  were  suddenly  removed  from  its  field-of-view.  This  quantity  is  composed  of 
radiance  from  the  sea  surface  plus  radiance  reflected  from  the  sky  (upper  atmos¬ 
phere).  Radiance  is  received  from  the  sky  even  if  it  Is  not  in  the  field-of-view 
because  of  reflecting  facets  caused  by  sea  waves.  The  sea  surface  radiance  is 
given  by  equation  (4)  as  follows, 

Nb  -  R^t w+  Ms(l-tw)  (4) 

where  all  radiative  quantities  are  a  function  of  wavelength.  The  first  term  in 
equation  (4)  is  the  radiance  of  the  sea  surface.  The  graybody  assumption  is  made 
so  that  this  quantity  can  be  given  by  the  product  of  blackbody  radiation  at  the  sea 
surface  temperature  and  the  sea  surface  emissivity  which  is  assumed  constant  with 
wavelength.  However,  it  is  well-known  that  water  in  the  Infrared  (2-15um)  exhibits 
a  rapid  decrease  in  emissivity''  as  the  angle  of  incidence  is  increased  above  50c . 

A  curve  fit  of  measured  data  and  theoretical  calculations  of  unpolarized  reflec¬ 
tions  was  made  to  determine  the  change  in  effective  sea  surface  emissivity  with 
viewing  angle.  The  following  equation  results7: 

=  0.98(1  -  (1  -  cos  0)5)  (5) 


"Kreith,  F.,  "Principles  of  Heat  Transfer,"  International  Textbook  Company, 
Scranton,  PA,  Second  Edition,  p.  210,  April  1955 

5Hudson,  R.  0.,  Jr..,  "Infrared  System  Engineering,"  John  Wiley  and  Sons,  New  York 
page  35  (1965) 

6Wo1fe,  W.  1.,  "Handbook  of  Military  Infrared  Technology,"  Office  of  Naval 
Research,  Washington,  DC,  page  167  (1966) 

’"Fleet  Signature  Computer  Model  Program  Manual,”  Op.  Cit.,  p.  5-1 
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where  y  is  the  angle  between  the  1 ine-of-sigt t  from  the  target  to  the  detector  and 
the  normal  to  the  sea  surface  facets.  The  determination  of  sea  surface  emissivity 
depends  on  knowledge  of  sea  waves  to  predict  the  angle  of  the  sea  surface  facets. 
Surface  waves  may  reflect  light  from  the  sk>  over  a  wide  range  of  angles.  However, 
Hulburt6  concluded  from  his  data  that  the  T *j ght  reflected  by  a  "breezy  sea"  comes 
mainly  from  a  region  of  the  sky  that  is  25°  to  35°  above  the  horizon.  Hence,  it 
was  concluded  that  for  wind  velocities  between  5  and  25  knots,  the  sea  facets  at 
the  horizon  can  be  represented  by  smooth  water  tilted  15°  towards  the  observer. 

For  water  surfaces  before  the  horizon  the  average  tilt  angle  must  be  reduced  to 
account  for  negative  wave  slopes,  i.e.,  waves  tilted  away  from  the  detector  but  of 
sufficiently  small  slope  to  allow  light  to  be  reflected  towards  the  detector.  Ex¬ 
perimental  data  was  used  to  conclude  that  this  reduction  in  wave  slope  could  be 
given  by  half  the  target-detector  elevation  angle.  The  target  detector  elevation 
angle  is  given  as  90°  minus  the  target  zenith  angle,  Zp.  The  target  zenith  angle 
is  the  angle  perpendicular  to  the  sea  at  the  ship  and  the  1 ine-of-sight  to  the 
detector.  The  angles  which  determine  the  effective  sea  surface  emissivity  are 
shown  in  Figure  1. 


Figure  1.  Angles  Determining  Effective-  Sea  Surface  Emissivity 


I 

« 


The  second  term  in  equation  (4)  is  the  radiance  reflected  by  the  sea  surface 
facets  towards  the  detector.  The  sky  radiance.  Ns,  depends  on  the  amount  of  cloud 
cover.  Hence,  it  is  assumed  that  sky  radiance  is  composed  on  the  average  of  a 
portion  from  the  sky  which  is  dear  and  a  portion  from  sky  that  is  cloud  covered. 
This  is  expressed  in  equation  (6)  as  follows, 

Ns  »  {1-  f)  •  Nc+  f  •  Sa  (6) 


"Hulburt,  E.  0.,  "The  Polarization  of  Light  at  Sea,"  J.  Opt.  Soc.  Am.,  Vol .  24 
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where  all  radiative  quantities  are  a  function  of  wavelengths.  The  fractional  part 
of  clear  sky,  f,  is  relatable  to  the  cloud  cover  which  is  Included  in  meteorologi¬ 
cal  data.  The  cloudy  sky  radiance,  Nc ,  is  composed  of  two  parts.  The  first  part 
is  the  emission  of  the  cloud,  Rc,  taken  to  be  a  blackbody  at  the  air  temperature  at 
the  cloud  height.  This  radiance  is  attenuated  by  the  atmospheric  transmission,  Ua, 
of  the  path  between  the  clouds  and  the  target.  The  second  part  is  the  spectral 
radiance.  Sc,  of  the  atmosphere  contained  in  the  path  between  the  clouds  and  the 
target.  Thus,  the  cloudy  sky  radiance,  Nc,  is  the  sum  Rc  •  Ua  +Sc.  The  clear  sky 
radiative  emission,  Sa ,  is  the  spectral  radiance  of  the  atmosphere  in  a  clear  path 
from  the  target  to  the  top  of  the  atmosphere.  The  height  of  the  clear  sky  is 
arbitrarily  taken  to  be  30  km  in  calculating  clear-sky  radiance.  The  empirical 
result  that  most  of  the  sky  radiance  from  a  reflecting  sea  surface  comes  from  a 
region  about  30c  above  the  horizon  is  used  to  locate  sky  radiance.  It  is  assumed 
that  sky  radiance  is  reflected  from  wave  facets  at  the  sea  surface  which  are 
inclined  15:  from  the  horizontal.  This  wave  slope  was  modified  slightly  for  targets 
before  the  horizon  as  previously  explained.  A  graphical  description  of  the  com¬ 
ponents  of  sea  surface  radiance  is  given  in  Figure  2. 
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equation  (6)  generally  requires  a  radiance  calculation  along  two  paths  of  different 
lengths  but  similar  direction.  The  general  equation  for  the  apparent  spectral 
radiance  observed  from  one  end  of  a  column  of  atmospheric  gases  Is  given  by  the 
following  equation. 

t 

S  =  /  Ua  •  R  •  c  •  dt  (7) 

¥  da 

0 

This  integration  is  more  complicated  over  paths  of  changing  altitude  since  all  three 
terms  in  the  integrand  vary  with  altitude.  The  integral  was  simplified  by  a  trans¬ 
formation  of  variables  to  radiance  as  the  independent  variable10.  The  new  Integral 
was  solved  numerically  using  Guass's  quadrature  formula  in  the  fifth  degree.  This 
formulation  requires  that  air  radiance  (a  blackbody  at  air  temperature)  be  calcu¬ 
lated  at  five  altitudes.  In  calculating  air  radiance,  the  following  relationship 
between  air  temperature  and  altitude  was  assumed. 

H  ■.  1 1  km  Tfl  3  To  -  72 

H  '  11  km  T  3  To  -  72  •  H/ll 

3 

The  spectral  radiances  defined  by  equations  ( 2 ) - ( 7 )  were  calculated  at 
equally  spaced  wavelength  intervals  3  0.4i.m)  in  the  8-l£.m  region.  They  were 
then  converted  to  a  band  average  radiance  by  simply  multiplying  by  the  wavelength 
interval  and  summing  the  spectral  quantities. 

The  optical  transmissions  along  the  appropriate  air  paths  are  needed  in  the 
computation  of  both  clear  and  cloudy  sky  radiances.  In  the  8-12um  waveband,  thermal 
radiation  is  attenuated  by  discrete  absorption  lines  of  water  vapor,  ozone,  and  a 
combination  of  uniformly  mixed  gases  (CO,  CO?,  N;0,  etc.)  and  by  a  continuum 
absorption  due  to  water  vapor.  Additional  energy  is  absorbed  by  naturally  occurring 
aerosols  which  are  introduced  into  the  air  by  continental  and  surface  sources.  The 
present  work  updated  the  existing  ship  signature  codes  by  Improving  the  calculation 
of  atmospheric  transmission  In  the  8-12. m  region.  Transmission  assocaited  with 
molecular  absorption  due  to  water  vapor,  ozone  and  uniformly  mixed  gases  were  cal¬ 
culated  using  tables  of  absorption  coefficient  versus  wavenumber  taken  from  the 
LOWTRAN  III  Computer  Code’1'*.  Water  vapor  continuum  absorption  was  calculated  from 
an  empirical  law1'  derived  for  the  newer  Jll-b  version  of  LOWTRAN.  The  reduction 
in  transmission  due  to  the  presence  of  aerosols  was  recently  modeled13  by  combining 
predictions  of  marine  and  continental  aerosol  particle  distributions.  The  total 
transmission  for  a  given  wavelength  is  the  product  of  the  above  Individual  trans- 
mi  ssions. 


(5F) 

C’F) 


(8) 


? 


•Ibid,  p,  4-2 

'Selby,  0.  E.  and  HcClatchey, 


Computer  Code  LOWTRAN 
1JSelby,  J.  E. ,  et.  al . , 
LOWTRAN  38  (1976),"  Air 
1 ’Katz,  8.  S.  and  Hepfer, 
Environments,"  Proc.  of 
April  1978 


R.  A.,  "Atmospheric  Transmittance  from  0.25  to  28.5pm: 
Air  Force  Cambridge  Res.  lab.,  AFCRl-TR-75-0255,  May  1975 
'Atmospheric  Transmittance  from  0.25  to  28.5pm:  Supplement 
Force  Geophysics  Lab.,  AFGl-TR-76-0258,  November  1976 
K.,  "Electro-optics  Systems  Performance  In  Selected  Marine 
AGARO  Symp.  Model.  Aerospace  Prop.  Envir.,  Ottawa,  Canada, 
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111.  HEAT  TRANSFER  RATES  USED  IN  SINGLE  ELEMENT  METHOD 


Ship  signatures  are  defined  in  this  report  to  be  the  contrast  temperature  dif¬ 
ference  between  the  ship  and  its  infrared  background.  This  section  will  present 
methods  of  calculating  the  heat  transfer  rates  which  produce  the  temperatures  of 
the  Single  Element  Method.  The  Single  Element  Method  consists  of  calculating  the 
temperature  in  a  single  element  which  represents  a  ship  and  applying  three  cor¬ 
rection  factors  to  account  for  differences  between  the  element  and  the  ship.  The 
basic  thermal  model  of  this  method  for  predictions  of  this  report  is  a  single 
vertical  plate  of  unit  surface  area.  This  element  is  heated  by  solar  energy  and  by 
internal  heat  sources  such  as  the  power  plant  or  auxiliary  machinery.  The  element 
is  cooled  by  air  flow  due  to  ship  and  wind  motion  and  by  occasional  rainfall. 

Solar  heating  is  the  largest  heating  rate  and  is  Important  because  it  dominates 
ship  temperatures  during  the  day  and  may  influence  these  temperatures  for  several 
hours  after  sunset.  The  total  solar  heat  load  incident  on  a  given  surface  is  the 
sum  of  direct  and  diffuse  solar  energy  multiplied  by  the  solar  absorptivity  of  its 
surface  coating.  The  small  amount  of  reflected  solar  energy  due  to  the  albedo  of 
the  sea  will  be  neglected.  Direct  solar  radiance  is  simply  the  solar  constant  or 
solar  radiation  received  outside  the  atmosphere  (1350  watts/m7)  multiplied  by  the 
transmittance  of  solar  energy  in  traversing  the  atmosphere,  i.e., 

Qo  -  1350  Ut  (9) 

The  total  transmission,  Ut,  consists  of  4  components  which  account  for,  (1)  absorp¬ 
tion  of  water  vapor,  (2)  molecular  scattering,  (3)  aerosol  scattering,  and  (4)  the 
masking  of  solar  rays  by  clouds.  Equations  for  these  components  are  given  in  the 
references 1 11 .  The  diffuse  or  indirect  solar  radiance  is  due  to  the  aerosol  and 
molecular  scattering  of  direct  solar  radiation  out  of  the  1 ine-of-sight  between  the 
sun  and  a  surface  on  the  earth.  It  is  assumed  that  half  of  what  is  scattered 
eventually  reaches  this  surface  on  the  average.  From  this  assumption,  the  following 
equation  for  diffuse  solar  radiance  is  derived. 


4 


Qf  *  Qo  (1  -  Urns  +  Uas)/2)  (10) 

The  amount  of  direct  or  diffuse  solar  energy  that  impinges  on  a  surface  depends 
on  the  angular  relationship  between  the  1 ine-of-sight  from  surface  to  sun.  This 
relationship  is  different  for  direct  and  diffuse  radiance  but  both  depend  on  the 
position  of  the  sun  relative  to  the  earth.  The  sun’s  location  relative  to  the  earth 


1 '“’Fleet  Signature  Computer  Model  Program  Manual,"  Op.  Cit.,  p.  6-2,  6-3 
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•is  described  by  the  sun's  zenith  and  azimuth  angles.  The  sun's  zenith  is  defined  I 
relative  to  a  perpendicular  to  the  earth's  surface  and  its  azimuth  is  defined  1 
relative  to  North.  The  following  equations  define  the  solar  zenith  and  azimuth  i 
angles1 s. 


cos  z  =  (sin  c)(sin  d)+  (cos  c)(cos  d)(cos  b)  (11) 

cos  Az  =  (cos  c)(sin  d)-  (cos  d)(sin  c)(cos  b)  (12) 

UTrTzT^ 

In  these  equations,  the  solar  declension,  d,  is  a  function  of  time  of  the  year  and 
the  hour  angle,  b,  is  15°  for  every  hour  of  time  away  from  either  side  of  solar 
noon.  Once  the  sun  has  been  located,  it  is  only  an  exercise  in  geometry  to  find 
the  direct  solar  radiance  which  is  given  by  the  component  of  sunshine  normal  to 
the  surface.  Figure  3  shows  the  geometrical  considerations  in  determing  the  normal 
component  of  sunshine  for  a  vertical  element. 


I 


f 

k 


'Threlkeld,  J.  L.,  "Thermal  Environmental  Engineering," 
Englewood  Cliffs,  New  Jersey,  p.  319  (1962) 
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The  following  equation  for  direct  solar  radiance  on  a  vertical  element  repre¬ 
senting  a  ship  of  heading,  M  is  derived  following  Figure  3. 

Qd+  Qo  ■  sin  z  •  cos(M  *  90°  -  Az)  (13) 

In  equation  (13),  +90°  is  used  for  the  portside  and  -90°  for  the  starboard  side  of 
the  ship.  The  heading  M  is  measured  clockwise  (eastward)  from  north  (heading  is 
shown  counterclockwise  in  Figure  3  to  avoid  cluttering). 

The  diffuse  component  of  solar  radiance  is  Incident  on  a  surface  from  all 
directions  and  is  found  by  integrating  the  dot  product  of  the  surface  normal  and  a 
distribution  function  for  the  intenstiy  of  diffuse  energy.  The  following  result 
was  dervied* 

Qi  ,  Qf  (1+  2-  sin  z  •  cos  t  90  -  Az!)  (14) 

Here,  +0''  is  for  a  portside  element  and  -90  for  a  starboard  element  in  agreement 
with  elation  (13).  In  summary,  the  total  solar  energy  on  a  vertical  surface  is 
given  by  the  following  equation 

Qs  «  a(Qd  *  Qi )  (15) 

Heating  due  to  internal  sources  is  described  by  providing  the  internal  tempera¬ 
ture  and  a  convective  heat  transfer  coefficient  for  each  compartment  of  the  ship. 
These  temperatures  and  coefficients  are  found  in  ship  design  specifications  and 
typical  values  are  chosen  for  use  in  the  single  element  method. 

Convective  cooling  iates  due  to  air  flow  over  the  ship  are  generally  lower 
than  those  due  to  rainfall  but  are  more  important  because  they  always  contribute  to 
the  ship  temperature  calculated  using  weather  data.  Cooling  due  to  wind  is  de¬ 
scribed  by  a  convective  coefficient  which  was  derived  from  measurements  made  on 
ships  at  sea.  The  following  equation  results  for  vertical  surfaces** 

ha  =  3.56  •  V0-6  (16) 

Here,  V  is  the  relative  wind  speed  which  is  the  vector  sum  of  the  wind  velocity 
vector  (found  in  weather  data)  and  the  ships  speed  and  heading  vector. 

Cooling  rates  for  rainfall  are  also  described  by  a  convective  coefficient 
which  was  derived  for  this  work  using  boundary  layer  analysis.  Before  proceeding 
with  this  analysis,  it  is  necessary  to  classify  rain  conditions  based  on  the  rate 
of  rainfall.  Light  rain  was  defined  to  be  0.1  in/hr  or  less  and  heavy  rain  0.3 
in/hr  or  greater  in  reference16.  The  boundary  layer  analysis  was  done  for  each 
of  these  rain  rates  in  the  following  manner,  first,  it  was  assumed  that  these  rain 


*Th!s  equation  was  derived  by  Hr.  0.  Friedman  of  the  Naval  Research  Laboratory  who 
provided  much  of  the  analysis  going  into  the  Ship  Signature  Computer  Codes  pre¬ 
viously  mentioned. 

**TMs  equation  also  provided  by  Mr.  D.  Friedman 

16"Visible  and  Infrared  Transmission  Through  Clouds,  Fog,  and  Rain,"  Raytheon  Missile 
Systems  Div.,  3R-4029,  December  30,  1966 
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rates  fall  onto  and  flow  along  a  3-foot-wide  deck  and  collects  on  Its  edge.  This 
forms  the  Initial  conditions  to  a  laminar  flow  boundary  layer  computer  program17 
which  calculates  the  velocity  profile  for  water  flowing  down  a  vertical  plate.  An 
average  velocity  is  found  from  the  fully-developed  velocity  profile  and  this  is 
used  to  compute  an  average  convective  coefficient18.  In  these  calculations  the 
ship  side  was  taken  to  be  10  feet  in  height,  Coefficients  were  found  for  the 
maximum  light  rain  and  minimum  heavy  rain  conditions  assumed  above.  The  minimum 
heavy  rain  coefficient  was  Increased  10%  and  the  maximum  light  rain  coefficient  was 
decreased  10%  to  realize  coefficients  for  representative  light  and  heavy  rain.  An 
intermediate  value  was  chosen  as  a  medium  rain  condition.  Many  weather  conditions 
occur  as  drizzle  which  have  less  rainfall  than  light  rain.  Convective  coefficients 
for  drizzles  are  provided  by  assuming  these  are  15%  of  the  corresponding  rain 
values.  This  arbitrary  assumption  Is  based  on  the  fac*  that  free  convective 
coefficients  are  generally  10-20%  of  forced  convective  .oefficlents,  The  following 
table  presents  the  rain  coefficients  derived  for  this  work. 


Table  1.  Convective  Cooling  Coefficients  for  Rain  Falling  on  a  Vertical  Element 


Intensity 

Drizzle 

Rain 

Light 

12  watts/m" -rC 

80  watts/m:'-°C 

Med i urn 

20 

135 

Heavy 

30 

o 

o 

c\j 

’Wilson,  0.  M,  and  Katz,  8.  $.,  "The  Use  of  Water  Cooling  for  Protection  Against 
Thermal  Radiation  from  a  Nuclear  Weapon  Detonation N0LTR  74-59,  23  Apr  1974 
’Kreith,  K.,  Op.  Cit.,  p.  296 


NSWC/WOl  TR  78-187 


IV.  SINGLE  ELEMENT  METHOD  FOR  CALCULATION  OF  SHIP  TEMPERATURE 

Equations  (9)-(16)  for  heat  transfer  rates  will  be  used  in  the  Single  Element 
Method  to  calculate  the  average  temperature  of  a  target  ship.  In  this  section,  the 
basic  temperature  equation  and  the  equations  for  the  3  correction  factors  which 
define  this  method  will  be  developed.  The  Single  Element  Method  will  be  applied  to 
a  specific  ship,  the  U.S.  Navy  Patrol  Frigate,  and  used  with  weather  ship  data  to 
calculate  contrast  temperature  differences. 

The  Single  Element  Method  breaks  a  ship  into  sections  having  common  internal 
temperatures  and  subsequently  into  sections  having  comnon  thermal  capacity  of  outer 
walls,  A  basic  thermal  element  is  chosen  to  be  a  vertical  element  of  unit  surface 
area  and  the  lowest  internal  temperature,  Tnm^n,  and  lowest  thermal  capacity,  Wmjn, 
of  all  ship  sections.  The  basic  element  temperature  Is  calculated  and  correction 
factors  are  derived  from  the  differences  of  the  remaining  sections.  Correction 
factors  account  for  the  variations  in  internal  temperature  and  thermal  capacity  and 
for  the  presence  of  a  hot  stack.  An  equation  describing  the  temperature  of  the 
basic  element  is  found  from  a  heat  balance  on  this  element.  This  balance,  which 
equates  the  heat  stored  in  the  element  to  the  net  heat  transferred  to  it,  is  given 
by  equation  (17). 

Wmingy  5  Os~ha(T-Ta)-hr(T-Ta)-hi(T-Tnmin)  (17) 

Here,  Wmin  is  the  thermal  capacity  of  the  basic  element  which  is  the  product  of  its 
density,  specific  heat,  and  thickness.  Heat  conduction  rates  are  not  present  in 
equation  (17)  because  the  basic  element  is  thermally- th in ,  i.e,,  heat  is  assumed  to 
be  uniformly  distributed  from  its  front  to  its  rear  surface.  Also,  it  is  assumed 
that  each  section  is  thermally  isolated  from  its  neighbors,  i.e.,  there  is  no  thermal 
conduction  between  adjacent  elements. 

Equation  (17)  is  solved  for  the  basic  element  temperature  by  first  replacing 
the  derivative,  dT/d:  by  its  finite  difference  equivalent,  AT/A-.  Here,  AT  is  the 
change  in  element  temperature  from  an  initial  value,  To,  to  a  final  value,  Tb,  which 
results  after  a  time  interval,  \t.  The  following  non-dimensional  parameters  are 
defined  for  convenience. 


HA  » 

(18) 

HR  -  hru“' 

(19) 

HI  =  hi  -  At 

(20) 

W 
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A  final  grouping  of  variables  has  the  dimensions  of  temperature. 


(21) 


The  basic  element,  temperature  is  found  by  solving  equation  (17)  using  the  definitions 
expressed  in  equations  (18)-(21).  The  result  is 


Tb 


To  +  QS  +  (HA  +  HR)Ta  +HI  •  Tn 
1  +HA  +  HR  +  HI 


(22) 


Equation  (22)  is  an  equation  for  the  stepwise  numerical  computation  of  the  basic 
element  temperature.  That  is,  this  temperature  is  found  from  a  starting  temperature, 
Tc,  a  time  interval ,  and  weather  data  from  which  the  heat  transfer  rates  appearing 
in  equations  (i8)-(21)  can  be  calculated.  The  computed  temperature,  Tb,  then 
becomes  the  starting  temperature  for  a  subsequent  calculation  based  on  equation  (22) 
and  weather  data  for  the  succeeding  time  interval.  Equation  (22)  provides  a  stable 
calculation  of  temperature  regardless  of  choice  of  time  step.  This  result  was  not 
achieved  in  the  previous  ship  signature  computer  programs  where  a  time  step  of  1/16 
hour  was  typically  chosen  to  avoid  stability  problems. 

The  temperature  correction  factor  to  account  for  the  different  internal  tem¬ 
peratures  found  throughout  the  ship  is  developed  by  subdividing  a  side  of  the  ship 
into  sections  having  common  internal  temperatures.  The  side  of  the  ship  chosen 
corresponds  to  the  signature  while  viewing  that  side.  Next,  the  average  or  overall 
ship  temperature  is  the  average  temperature  of  all  of  these  sections  which  have 
different  internal  temperatures.  This  fact  is  expressed  by  the  following  equatior. 

T  =  A;  T  s  +  A,  T.  + _ —  AmTm  (23) 

^  +TTT+  — ' 


Here,  A,  is  the  area  fraction  or  fraction  of  the  total  visible  ship  surface  area 
which  has  a  temperature,  T  due  to  the  first  internal  temperature,  Tn).  A?  is  the 
fraction  of  visible  surface  having  a  temperature  T;  due  to  the  next  internal 
temperature,  Tn..,  and  so  forth  for  a  total  of  m  such  sections.  The  sum  of  the  area 
fractions.  A)  +  A.  *■  — A^  is  1.0  hence  the  denominator  of  equation  (22)  was 
included  for  clarity  only.  These  sections  are  identical  except  for  their  area  and 
internal  temperatures  hence  their  individual  temperatures  can  be  calculated  by 
equation  (22).  Equation  (22)  for  calculation  of  overall  section  temperature  is 
rewritten  below  in  a  form  which  isolates  the  internal  temperature  effect  which  is  to 
be  studied,  i  .e. ,  let, 


where 


Tj  =  D  +  E  •  Tn j 

(24) 

To  +  QS  +  (HA  *  HR)Ta 

(25) 

1  +HA  +  HR  +H1 

E  =  HI 

(26) 

1  +  HA  ♦  HR  +  HI 

and  j  is  an  index  running  from  1  tc  m  (m  is  the  total  number  of  sections).  Now, 
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when  equation  (24)  is  substituted  in  equation  (23)  for  the  section  temperatures, 
these  are  eliminated  and  the  average  ship  temperature  is  derived  In  terms  of  inter¬ 
nal  temperatures. 


T  *  D  +  E  In,  +  A  E(Tn,  -  Tn,)+  —  AnE(Tnm-  Tn, )  (27) 

basic  temperature  correction  factor 

element  for  internal  temperature  variation 

temperature 

Equation  (27)  was  designed  to  divide  into  a  basic  element  temperature  plus  some 
terms  depending  on  area  fractions  and  internal  temperature  differences  of  the 
various  ship  sections.  The  first  two  terms  define  the  average  temperature  of  the 
entire  ship  (area  fraction  =  1.0)  if  it  had  internal  temperature,  Tn!t  which  is 
the  lowest  internal  temperature.  The  remaining  terms  represent  the  temperature 
correction  factor  for  internal  temperature  variations.  An  illustration  of  the  cal¬ 
culation  of  this  correction  factor  will  be  given  'ater  for  the  Patrol  Frigate. 

The  temperature  correction  factor  to  account  for  the  different  thermal  capa¬ 
cities,  W,  of  outside  ship  walls  are  developed  in  a  similar  manner  but  Independently 
of  the  internal  temperature  correction  method.  First,  the  side  of  the  ship  in  view 
is  subdivided  into  f  number  of  sections  each  having  a  different  thermal  capacity. 
Each  of  these  sections  will  have  a  common  internal  temperature  which  is  the  value 
of  the  section  which  has  the  lowest  internal  temperature.  The  average  or  overall 
ship  temperature  for  this  situation  is  expressed  by  the  following  equation. 

T  =  8. T-  +  B  T  +  —  BT  (28) 

•  •  •  •  p  p 

Here,  5.  is  the  area  fraction  having  thermal  capacity  W.  which  results  in  the  over¬ 
all  section  temperature  T..  8  is  the  fraction  of  visible  area  having  temperature 
T,  due  to  thermal  capacity,  W  ,  and  so  forth  for  a  total  of  p  such  sections.  The 
individual  temperature  for  each  of  these  sections  is  also  given  by  equation  (22). 
Equation  (22)  is  rewritten  below  in  a  form  which  isolates  the  thermal  capacity 
effect  to  be  studied,  i.e.,  let 


Tj  =■ 

J  t*  g  rvj 

(29) 

where 

Vj  =  Wj 

(30) 

F  =  Qs  ■*  (ha  +  hr)Ta  +  hi  Tn  mjn 

(31) 

G  =  ha  +  hr  +  hi 

(32) 

and  j  is  an  index  running  from  1  to  p  (p  is  the  total  number  of 
equation  (29)  is  substituted  into  equation  (28)  for  the  section 
are  eliminated  and  tne  following  equation  in  terms  of  variable 
resul ts . 

sections).  Now,  if 
temperatures,  these 
thermal  capacity 
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T  .  lojLU  + 

1  1  +  V  G 


basic  element 
temperature 


o  ,  iV;-V.)(F-  G-  To), 

b?-iV5(v.  +T.T+  g*v,v: 


p  .  (Vp~  V) H  F-  G  •  To) ,  temperature  cor-  (33) 
pir+GtVi  +  Vp)  +  Y^pG2 !  rection  factor  for 


thermal  capacity 
variation 


Equation  (33)  was  also  designed  to  divide  into  a  basic  element  temperature  plus 
terms  representing  the  temperature  correction  factor  for  thermal  capacity  variations 
of  the  remaining  sections.  The  basic  element  temperature  is  the  average  temperature 
of  the  entire  ship  if  it  had  the  thermal  capacity  function  Vj.  V)  represents  the 
lowest  thermal  capacity  found  on  ship  sections.  The  temperature  correction  factors 
for  internal  temperature  variations  [ equation  (22)!  and  for  thermal  capacity 
variations  (equation  (33)1  were  developed  independently  of  each  other  but  both  will 
be  added  to  the  basic  element  temperature  to  find  the  overall  ship  temperature. 

The  justification  for  doing  this  is  that  each  correction  is  small  compared  to  the 
basic  element  temperature  hence  no  appreciable  error  should  result  when  internal 
temperature  and  thermal  capacities  different  from  basic  element  values  exist 
simultaneously  on  some  sections.  The  size  of  these  correction  factors  will  be 
estimated  later  for  the  Patrol  Frigate. 


The  temperature  correction  factor  to  account  for  the  presence  of  one  or  more 
hot  stacks  is  developed  using  the  same  procedure  as  in  the  preceding  correction 
factors.  First,  the  overall  ship  temperature  is  taken  to  be  the  area  fraction 
multiplied  by  the  overall  temperature  of  the  ship  without  stacks  plus  a  contribution 
consisting  of  the  area  fraction  multiplied  by  the  overall  temperature  of  each  stack. 
This  first  temperature  is  taken  to  be  the  basic  element  temperature  (see  equation 
(22);  .  The  overall  ship  temperature  is  given  in  the  following  convenient  form  by 
eliminating  the  area  fraction  of  the  ship  without  stacks  from  a  preceding  equation 
s  imilar  tc  equation  (28) 


T  3  Tb  +  C,  ( T.  -  Tb)  +  C;(T;  -  Tb)  +  - Cq(Tq-  Tb) 


Here,  Tb  is  the  basic  element  temperature  for  an  area  fraction  of  1.0  and  Tlf  T ; , 
et'-.,  are  average  stack  temperatures  at  stack  area  fractions  C),  C; ,  etc.,  for  a 
total  of  q  stacks.  The  average  temperature  of  stacks  is  given  by  equation  (22) 
with  the  important  difference  that  these  elements  have  a  high  Internal  heat  trans¬ 
fer  rate,  hk,  and  a  high  interna!  temperature,  Tk,  replacing  hi  and  Tnm)n  of  the 
basic  element.  This  requires  that  the  internal  heating  parameter,  HI  (see  equation 
(20) !  be  replaced  by  the  stack  heating  parameter,  HK,  which  is  defined  as  follows: 


HK  *  -~ 


The  internal  temperature,  Tn,  in  equation  (22)  is  replaced  by  the  temperature,  TK, 
which  is  the  average  temperature  of  flue  gases  leaving  the  stack.  A  temperature 
correction  factor  for  stack  presence  can  be  derived  by  substituting  equation  (22) 
with  stack  modifications  into  equation  (34)  for  each  stack  present.  This  results 
in  a  complex  cumbersome  equation  for  stack  correction  which  can  be  simplified  by 
making  the  following  assumptions.  These  are:  (i)  the  thermal  capacity  difference 
between  the  stack  and  basic  element  is  neglected;  (2)  terms  involving  the  internal 
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heat  transfer  coefficient,  hi,  will  be  neglected  in  expressions  involving  terms 
containing  the  stack  heating  coefficient,  hK;  and  (3)  product  terms  containing  hi 
will  be  neglected  in  expressions  vhich  also  have  product  terms  containing  hK.  These 
assumptions  were  justified  by  calculations  based  on  the  Patrol  Frigate.  Using  these 
assumptions,  the  following  equation  for  the  average  ship  temperature  including  the 
presence  of  stacks  is  derived. 


T 


Tb  +  C1 


•  HK 


,  (HC  +  HR  + 1 ) ( T K i  -  Ta )  -  Qs  , 

1  '  (li  HR  ~+Ty  (H  C+HR~r  HK  yT  1 ' 


uv  ,(HC  +  HR  +  1  )  (  TK„  -  Ta)  -( 
ThTThTT  +  1)  (  H C  +  HR  +  H K, 


A  final  simplification  can  be  made  which  will  allow  the  stack  effect  to  be  rapidly 
estimated.  This  is  that  the  contribution  to  stack  temperature  by  solar  heating  is 
likely  to  be  small  in  comparison  to  the  heat  transferred  by  stack  gases.  In  this 
case,  the  solar  heating  parameter,  QS,  can  be  neglected  in  equation  (36)  which 
reduces  to  the  following 


.  r  .  HK(Tk-Ta)  , 
C:  (HC  ♦  HR  ♦  HKT 


Equation  (37)  is  the  estimated  average  temperature  of  a  ship  having  one  stack  or 
several  identical  stacks  and  it  can  be  used  to  initially  judge  what  effect  the 
stacks  will  have  on  the  total  ship  signature.  The  stack  correction  factor  considers 
the  effect  of  stack  gas  heating  but  does  not  consider  the  signature  contribution  of 
the  hot  gases  (plume)  leaving  the  stack.  This  is  justified  in  the  present  work 
for  SEAFIRE  systems  which  operate  in  the  8-12,.m  waveband  where  radiance  from  the 
predominately  carbon  dioxide  gas  in  the  exhaust  plume  is  rapidly  absorbed  by  the 
atmosphere.  A  system  operating  in  the  3-5.  m  waveband  would  receive  radiance  from 
the  plume  and  this  would  have  to  be  included  in  the  stack  correction  factor  for 
comp! eteness . 


The  correction  factors  for  internal  temperature  differences,  thermal  capacity 
differences  and  for  the  presence  of  stacks  are  given  in  equations  (27),  (33).  and 
(36),  respectively,  along  with  the  basic  element  temperature,  Tb.  The  single 
element  method  for  calculating  ship  temperature  consists  of  taking  the  average  ship 
temperature  to  be  the  basic  element  temperature,  Tb,  plus  the  sum  of  the  .3 
correction  factors.  The  basic  element  temperature  is  the  average  temperature  of  a 
given  plane  view  of  a  ship  if  its  entire  surface  area  had  the  lowest  internal 
temperature  and  lowest  thermal  capacity  of  all  major  sections  of  the  ship.  It 
should  be  remembered  that  an  average  ship  temperature  is  meaningful  in  ship  sig¬ 
nature  work  only  if  no  hot  spots  dominate  the  signature.  It  will  be  shown  later 
using  the  Patrol  Frigate  Model  that  internal  temperature  and  thermal  capacity 
differences  do  not  dominate  the  average  temperature  calculation.  However,  large 
stack  internal  heat  transfer  coefficients  and  high  stack  gas  temperatures  do  have 
the  potential  of  contributing  as  much  or  more  to  the  signature  as  the  remaining 
ship  for  stacks  of  even  relatively  small  surface  area.  If  this  is  the  case,  it 
makes  sense  to  apply  the  single  element  method  to  thn  stack  alone  and  neglect  the 
remaining  ship  surface.  It  also  can  be  applied  to  the  stacks  alone  for  the  case 
where  the  main  body  of  the  ship  is  obscured  by  the  horizon.  The  Single  Element 
Method  is  applied  by  using  the  stack  variables  in  equation  (22)  while  remembering 
that  the  resulting  temperature  applies  to  the  small  stack  surface  area  only.  The 
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correction  factors  do  not  apply  in  this  case.  The  stack  correction  is  usually 
small  when  signature  suppression  techniques  have  been  used  to  cool  the  exhaust 
gases  or  surrounding  stack  casing.  Equation  (37)  can  be  used  to  gauge  the  stack 
contribution  by  using  an  average  value  of  5.0  for  the  sum,  "MC  +  HR,"  with  a  time 
step,  of  one  hour  used  in  the  calculation  of  HK  i  see  equation  (35)  J. 

The  single  element  method  requires  a  ship  model  for  the  calculation  of  the 
basic  element  temperature  and  the  3  correction  factors  which  define  the  average 
ship  temperature.  A  model  of  the  U.S.  Navy  Patrol  Frigate  was  obtained  along  with 
the  SIRS  computer  code  ir.  which  this  model  was  used*.  This  model  breaks  the  ship 
into  92  sections  and  provides  the  internal  temperature  and  heat  transfer  coefficient, 
the  thermal  capacity,  emissivity,  absorptivity,  surface  area,  and  the  orientation 
of  each  section.  The  application  of  the  Patrol  Frigate  Model  in  the  Single  Element 
Method  is  given  in  Appendix  A.  This  appendix  shows  how  the  ship  is  broken  into 
sections  for  the  calculation  of  internal  temperature,  and  thermal  capacity  cor¬ 
rection  factors.  Also,  the  stack  is  analyzed  and  the  size  of  its  correction 
factor  estimated.  The  work  of  Appendix  A  will  later  be  used  with  weather  ship 
data  in  a  computer  code  to  calculate  the  ship  contrast  temperatures  and  their 
statistics. 

The  equations  developed  in  this  section  allow  the  actual  ship  temperature  to 
be  computed.  However,  for  ship  signature  work,  an  effective  temperature  corres¬ 
ponding  to  the  ship  radi'  ,ce  is  desired.  This  temperature  is  defined  as  the 
temperature  that  yields  the  same  radiance  in  a  graybody  calculation  as  the  target 
surface  radiance  calculated  by  equation  (2).  Since  sun  reflection  is  neglected 
and  the  apparent  background  is  the  air,  the  effective  ship  temperature  is  the 
fol lowing 


IS  =  c  T  +  ( 1  -  c, )Ta  (38) 

This  definition  is  based  on  the  expectation  that  the  Single  Element  Method 
Temperature,  T,  and  the  air  temperature,  Ta,  are  somewhat  the  same  so  that  the  sum 
of  radiances  in  Equation  (2)  is  approximately  reproduced  by  the  effective  ship 
temperature.  However,  the  emissivity  in  the  8-i2..m  waveband  of  Navy  Gray  paint 
used  with  the  Patrol  Frigate  Model  is  0.94, hence  the  effect  of  the  air  temperature 
ter^  is  small,  thereby  lessening  the  requirement  that  it  be  near  the  ship  temperature. 
The  problem  of  representing  sum  of  radiances  with  sum  of  temperatures  is  discussed 
in  the  next  section. 


*This  model,  along  with  the  SIRS  Computer  Program  and  instructions  on  its  use, 
was  obtained  from  Mr.  R.  Burns  of  the  Naval  Ship  Research  and  Development  Center 
(Annapolis). 
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V .  CONTRAST  TEMPERATURE  CALCULATION 

The  effective  ship  temperature  calculated  using  the  Single  Element  Method  Is 
combined  with  a  background  temperature  to  define  the  target  contrast  temperature 
that  Is  required  In  SEAF1RE  systems  analysis.  Ship  contrast  temperature  Is  based 
on  the  radiometric  definition  of  ship  signature  given  In  Equation  (2).  That  Is, 
the  contrast  temperature  at  the  target  ship  Is  equivalent  to  the  radiant  contrast, 
Nt-Nb,  which  appears  In  Equation  (1).  The  ship  contrast  temperature  Is  defined  to 
be  the  effective  ship  temperature  minus  the  background  temperature  as  given  by  the 
following  equation. 


TC  *  :S  -  TG  (39) 

The  ship  temperature,  TS,  Is  the  effective  temperature  which  approximately  produces 
the  target  Surface  radiance  when  used  In  a  blackbody  calculation  for  the  8-lZpm 
waveband.  Background  temperature  will  be  similarly  defined  such  that  It  will 
approximately  produce  the  background  radiance.  The  background  temperature  that  Is 
analogous  to  the  background  radiance  defined  In  Equation  (4)  Is  given  as  follows 

TG  -  c*  •  Tw*  (1  -i  J  -T y  (40) 

Mere,  Tw  and  Ty  are  the  sea  surface  temperature  and  the  effective  sky  temperature 
at  the  appropriate  reflection  angle.  the  sea  surface  emlsslvlty  defined  In 

Equation  (5).  ew  Is  about  0.24  for  viewing  one  nhlp  from  another.  Hence,  the 
background  raolance  and  therefore  Its  effective  temperature  Is  composed  of  about 
one  part  water  temperature  and  3  parts  sky  temperature.  Equation  (4)  applies  to 
the  case  where  the  target  ship  Is  slightly  before  the  horizon  so  that  sky  radiation 
will  be  reflected  off  of  the  wave  slopes.  If  the  ship  were  on  the  horizon,  the 
background  would  be  atmosphere  at  Its  sea  level  temperature.  Calculations  based  on 
Weather  Ship  J  data  showed  that  sky  temperatures  range  from  about  -40#C  to  +10°C 
while  sea  temperatures  are  close  to  +10bC.  It  will  be  shown  that  the  background 
radiance  computed  from  the  background  temperature  of  Equation  (40)  Is  not 
significantly  different  from  the  sea  surface  radiance  computed  from  Equation  (4), 
This  was  show:,  by  computing  radiances  by  each  method  and  comparing  them  for  sky 
tempera r.v.  •*.  r.?tween  -40°C  and  10°C  and  a  water  temperature  of  10®C.  The  result 
of  these  (vulations  which  were  made  for  the  8-12um  waveband  are  given  In  Figure 
4. 

Fi.yve  4  shows  ,i  maximum  difference  of  5  percent  between  radiances  at  -40°C 
and  this  : • ce  contiguously  decreases  as  sky  temperature  Increases.  It  will 
b r  shown  iaf  Hut  sky  temperatures  less  than  -40#C  art  rarely  encountered  end 
that  median  v  temperatures  are  In  the  0-5°C  region  for  the  weather  conditions  of 
this  report  .or.cu,  very  little  error  will  result  from  using  equet1on(40)  to 
define  the  background  temperature.  The  above  justification  of  background  Umpera- 
ture  depends  on  defining  a  sky  tempera ture .  Ty,  which  will  reproduce  the  sty  radiance 
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Ns,  which  was  defined  by  Equation  (6).  This  is  done  in  a  straight  forward  manner 
bel ow. 

Sky  temperature  is  determined  directly  from  a  sky  radiance  calculation.  The 
sky  radiance  is  composed  of  clear  sky  and  cloudy  sky  components  whose  calculation 
was  previously  outlined.  Sky  radiance  for  the  8-12i.m  region  was  found  by  summing 
the  product  of  spectral  sky  radiance  given  by  Equation  (6)  and  a  wavelength  in¬ 
crement  of  (M;,m  for  wavelengths  separated  by  0.4;;?n  in  the  8-12'im  waveband.  Sky 
temperature  is  defined  to  be  the  temperature  which  yields  this  radiance  in  a  black- 
body  calculation  made  in  a  like  manner.  The  relationship  between  temperature  and 
radiance  for  the  blackbody  in  the  8-12;:m  waveband  is  shown  in  Figure  5.  Sky  tempera¬ 
ture  could  be  found  by  computing  the  sky  radiance  and  reading  the  corresponding 
temperature  from  Figure  5  whenever  a  ship  signature  is  required.  The  previous 
ship  signature  codes  compute  a  signature  after  making  a  sky  radiance  calolation. 
However  this  procedure  will  not  be  used  here  because  sky  radiance  calculation  is  a 
lengthy  numerical  procedure.  This  procedure  as  previously  outlined  would  be  pro¬ 
hibitively  expensive  in  computer  time  for  the  thousands  of  contrast  temperatures 
calculated  herein  using  weather  data.  Instead,  sky  radiance  will  be  computed 
beforehand  and  provided  in  tabular  form  to  a  computer  code  which  calculates  contrast 
temperatures.  The  development  of  these  tables  is  described  next. 

Sky  radiance  is  composed  of  a  clear  sky  component  and  a  cloudy  sky  component 
which  depend  on  the  fraction  of  clear  sky,  f,  as  indicated  in  Equation  (6).  The 
calculation  of  both  of  these  components  require  that  the  atmospheric  emission  over 
the  paths  shown  in  Figure  2  be  determined.  This  requires  the  calculation  atmo¬ 
spheric  transmission  by  the  methods  given  previously.  The  above  review  is  given  to 
point  out  that  the  following  inputs  are  needed  in  making  up  the  sky  radiance  tables. 
These  are:  (1)  A  model  atmosphere  to  specify  conditions  of  the  atmosphere  above 
sea  level;  (2)  cloud  heights  to  determine  the  length  of  cloud  radiance  paths,  and 
(3)  wind  speeds  which  are  required  in  aerosol  transmission  calculations.  The  re¬ 
lative  humidity  is  important  to  transmission  due  to  water  vapor  and  water  vapor 
continuum  hence  this  variable  should  influence  sky  radiance  calucation.  However, 
it  was  found  that  specifying  a  sea  level  relative  humidity  has  only  a  very  small 
effect  on  sky  radiances  which  are  governed  by  water  vapor  amounts  specified  by  the 
model  atmospheres.  Hence,  relative  humidity  cannot  be  included  in  the  sky  radiance 
tables  until  a  method  of  relating  sea  level  and  higher  altitude  relative  humidities 
is  developed.  Two  model  atmospheres  were  chosen  in  anticipation  of  the  weather 
data  which  will  be  taken  from  a  representative  weather  ship.  This  ship,  Weather 
Ship  J,  operated  off  the  coast  of  Ireland  (latitude  52CN,  Longitude  20^)  hence 
model  atmospheres  for  a  midlatitude  summer  and  a  midlatitude  winter  were  chosen 
from  the  LOWTRAN  transmission  work’*.  Five  cloud  heights  and  six  wind  speeds  were 
chosen  to  cover  the  range  of  expected  weather condi tions .  The  clear  and  cloudy 
sky  radiances  which  make  up  the  sky  radiance  tables  are  given  in  Table  2.  The 
radiances  in  Table  2  are  calculated  using  a  sea  level  air  temperature  of  1 0° C . 

This  value  approximates  the  year-round  air  temperatures  observed  by  Weather  Ship  J. 


A  computer  program  to  compute  ship  contrast  temperature  was  written  using  the 
methods  of  this  report.  This  program  starts  with  weather  ship  data  and  computes 
ship  temperature,  sky  temperature,  background  temperature  and  contrast  temperature 
for  many  data  points.  It  then  computes  the  probability  of  occurrence  of  given 
values  of  temperature  for  all  of  these  ship  signature  variables.  This  program  was 

’’Selby,  R.  E.  and  McClatchey,  R.  A.,  Op.  Cit. 
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designed  to  analyze  large  amounts  of  data  economically  with  respect  to  computer 
running  time.  For  example,  six  months  of  hourly  weather  data  is  analyzed  in  180 
seconds  of  COC  6*500  computer  time.  This  includes  calculating  and  printing  data 
hourly  and  calculating,  printing  and  plotting  the  statistical  data  on  a  monthly 
basis. 


Table  2 .  Sky  Radiance  Tables  for  Surface  Temperature  of  10°C. 


Wind 

Clear; 

Cloudy  Sky  Radiance* 

Speed 

=  1750m  I' 

=  1250m 

H  *  800m  H 

*  400m 

H  «  150m 

r 

0. 

1 

16.25 i 

j 

26.75 

27.73 

28.86 

29.93 

30.78 

,  5.0  m/sec 

16.46 

26.77 

27.80 

28.87 

29.98 

30,78 

10.0  m/sec 

17.22 ' 

26.84 

27.84 

28.89 

29.99 

30.78 

1 S. 0  m/sec 
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VI.  RESULTS 
DaJLd  Jenna  t 

The  AS1RCT  Computer  Program,  described  in  Appendix  8,  was  run  with  weather  data 
from  Weather  Ship  J  to  compute  contrast  temperatures  using  the  method  of  this  report. 
These  temperatures  were  analyzed  statistically  to  find  cumulative  probabilities  of 
occurrence  of  temperature  which  was  the  form  required  for  SEAFIRE  systems  analysis. 
The  contrast  temperatures  reported  are  for  the  8-l?,.m  region  of  the  infrared  since 
the  SEAFIRE  systems  as  well  as  most  present  modular  FLIR  systems  operate  in  this 
waveband.  Ship,  sky  and  background  temperature  are  the  component  temperatures 
which  define  the  contrast  temperatures.  These  are  computed  to  aid  understanding 
of  contrast  temperature  results  and  to  tie  in  the  contrast  temperature  and  contrast 
radiance  definitions  of  ship  signature  All  calculations  were  made  using  the  model 
of  the  U.S.  Navy  Patrol  Frigate  which  is  given  in  Appendix  A.  It  was  also  shown  in 
Appendix  A  that  results  computed  using  this  model  will  provide  a  good  approximation 
for  all  similarly  constructed  ships  having  cooled  stacks.  The  results  of  contrast 
and  component  temperatures  will  be  presented  for  8  years  of  Weather  Ship  J  operation. 
These  temperatures  are  given  primarily  for  a  target  ship  heading  which  results  in  a 
maximum  contrast  (due  cast  or  west).  However,  contrast  temperatures  for  other 
headings  will  be  given  and  the  effect  of  heading  on  contrast  will  be  discussed. 

This  section  will  also  provide  the  rationale  for  defining  probabilities  and  results 
will  be  given  in  terms  of  these  probabilities.  A  method  of  relating  the  ship 
contrast  temperature  and  ship  radiance  signature  will  also  be  given.  Finally,  ship 
contrast  temperatures  will  be  given  for  a  background  consisting  of  air  at  sea  level 
temperature.  This  case  night  be  expected  for  a  target  ship  slightly  above  the 
horizon  where  there  is  no  reflection  from  the  sea  in  the  background. 

Statistical  analysis  of  temperatures  calculated  in  this  report  consists  of 
finding  the  cumulative  probability  of  occurrence  of  all  temperatures  that  are  equal 
to  or  less  than  a  given  value.  This  definition  is  useful  in  system  studies  where 
a  given  system  is  assumed  to  be  operational  at  a  given  contrast  temperature.  Then, 
the  probability  of  system  success  is  equal  to  the  probability  that  the  contrast 
temperature  is  above  the  critical  value.  The  results  are  produced  by  compiling  the 
individually  computed  ship  contrast  temperatures  into  a  number  of  standard  tempera¬ 
ture  values.  The  standard  values  are  provided  in  a  list  of  temperatures  separated 
by  a  0.2:C  interval  between  ‘3.0:C  and  slightly  larger  intervals  outside  this  range. 
This  list  can  be  found  in  Figure  8-2  which  is  a  listing  of  the  ASIRCT  computer 
program  (see  TliST).  A  computed  temperature  is  compared  sequentially  to  this  list 
until  the  first  smaller  temperature  is  encountered.  The  computed  temperature  is 
stored  at  this  list  temperature.  This  process  is  repeated  for  many  calculated 
temperatures  until  a  large  number  of  points  have  been  stored  in  various  positions 
on  the  standard  list.  A  probability  density  function  for  ship  temperature  is 
produced  by  dividing  the  number  in  each  temperature  of  the  list  by  the  total 
number  of  temperatures  stored.  The  probability  that  temperature  is  equal  to  or 
lower  than  the  listed  temperature  is  computed  sequentially  by  use  of  the  probability 
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density  function  starting  with  the  highest  temperature.  Calculation  is  continued 
until  probabil ities  for  all  temperatures  in  the  standard  list  have  been  computed. 
These  cumulative  probabilities  are  printed  as  tables  in  the  computer  output  for 
each  of  the  contrast  and  component  temperatures  mentioned  in  this  report.  However, 
to  simplify  the  presentation  of  data, temperatures  corresponding  to  standard 
probabilities  from  0.10  to  0.95  were  computed  by  interpolating  the  computer  printed 
tables.  The  standard  probabilities  and  their  corresponding  temperatures  were  also 
printed  in  *.he  computer  output  and  are  reproduced  in  this  report. 

Ship  Heading  and  Solar  Insolation 

The  ship's  heading  influences  the  amount  of  solar  energy  absorbed  on  a  vertical 
element  see  Figure  3  and  equations  (13)-(15)j.  It  has  a  significant  effect  on 
daytime  signatures  since  solar  heating  is  a  major  factor  in  rasising  the  bulk  ship 
temperature  above  its  ambient  air  temperature  value.  Daytime  contrast  temperatures 
ar*  very  important  in  SEAFIRE  systems  studies  where  FL I R  systems  are  compared  to 
television  systems  for  passive  daytime  operation.  Solar  Insolation  is  the  rate  of 
heating  received  at  the  earth  from  the  sun  hence  an  element  aligned  perpendicularly 
to  the  1 ine-of-sight  to  the  sun  receives  the  maximum  possible  heating.  The  line- 
of-sight  from  earth  to  sun  also  depends  on  the  time  of  year  due  to  the  variation  in 
the  sun's  declination.  Hence,  the  heating  of  a  ship  which  is  represented  by  a 
vertical  element  depends  on  the  orientation  of  the  element  with  respect  to  the  1 ine- 
of-sight  to  the  sun.  The  maximum  solar  insolation  occurs  on  June  21  when  the  sun's 
declination  is  a  maximum  and  minimum  insolation  occurs  on  December  21.  The  sun 
rises  in  the  east  and  traverses  the  southern  sky  in  a  westward  direction.  Hence,  a 
ship  heading  due  east  will  receive  maximum  or  all  day  direct  heating  on  its  star¬ 
board  side  and  a  ship  heading  due  west  will  receive  maximum  heating  on  its  port- 
side.  Ships  heading  due  south  (starboard)  or  due  north  will  receive  direct  solar 
heating  in  the  morning  and  afternoon,  respectively,  and  indirect  or  diffuse  solar 
heating  the  remaining  time  when  its  sides  are  shaded.  Ships  heading  due  west 
(starboard)  or  due  east  (portside)  are  always  shaded  hence  receive  diffuse  solar 
heating  only.  These  3  headings  for  an  element  on  the  starboard  side,  due  east, 
due  south,  and  due  west,  which  have  full,  half,  and  no  direct  solar  heating  will 
provide  the  reference  headings  for  this  report.  The  daily  incident  solar  energy  on 
a  starboard  vertical  element  for  each  of  these  headings  was  computed  by  summing 
hourly  the  direct  and  indirect  heating  rates  given  by  equations  (13)  and  (14). 

These  are  shown  in  Figure  6  for  the  52  N  longitude  of  Weather  Ship  J.  Figure  6 
shows  that  a  factor  of  3  or  more  difference  in  daily  incident  solar  energy  exists 
between  the  minimum  (270  heading)  and  maximum  (90'  heading)  values.  The  dip  in 
solar  heating  received  by  the  vertical  element  when  solar  insolation  is  a  maximum 
(June  21)  is  due  to  the  relatively  large  angle  between  the  perpendicular  of  the 
element  and  the  1 ine-of-sight  to  the  sun.  Contrast  temperature  calculations  given 
later  will  show  that  solar  heating  differences  are  important  to  daytime  contrast 
temperatures.  However,  the  bulk  of  the  contrast  temperature  calculations  of  this 
report  are  for  the  heading  producing  maximum  solar  heating  (90°). 

Temporal  Division  of  Data 

The  variation  of  incident  solar  heating  with  time  of  the  year  suggests 
breaking  the  year  into  seasons.  Three  seasons  were  defined.  These  are:  (1)  a 
winter  season,  January,  February,  November  and  December  having  minimum  solar 
heating;  (2)  a  spring-fall  season,  March,  April,  September  and  October  having 
intermediate  solar  heating;  and  (3)  a  summer  season.  May,  June,  July  and  August 
having  maximum  solar  heating.  Solar  heating  is  seen  to  decrease  in  Figure  6  over 
the  summer  month  for  the  due  east  heading  due  to  a  loss  in  angular  efficiency  caused 
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by  high  solar  altitude  angles.  But,  even  here  there  is  a  characteristic  dif¬ 
ference  between  the  summer  and  spring-fall  seasons  even  though  their  total  solar 
heating  is  about  the  same.  That  is,  the  spring-fall  season  has  higher  mid-day 
heating  due  to  the  better  angular  alignment  between  the  sun  and  the  vertical 
element.  It  will  be  shown  that  this  will  result  in  significantly  higher  daytime 
contrast  temperatures. 

Ship  contrast  and  its  component  temperatures  were  computed  using  data  gathered 
by  Weather  Ship  J  during  the  years  1964-71.  These  will  initially  be  presented  and 
discussed  for  the  due  east  ship  heading  and  for  the  3  seasons  defined  above. 
Subsequential Iy ,  representative  results  will  be  given  for  the  due  south  and  due 
west  headings  shown  in  Figure  6.  For  the  due  east  heading,  tables  of  the  standard 
probabilities  will  be  given  for  each  temperature,  each  season,  and  each  of  the  8 
years.  A  representative  average  value  of  the  8  years  will  be  plotted  for  each 
season  and  for  each  temperature. 

Effective  Ship  Temperature  Definition 

The  effective  ship  temperature  is  defined  to  be  a  thermodynamic  measure  of  the 
ship  target  surface  radiance.  It  is  the  temperature  of  a  blackbody  [equation  (3)1 
that  radiates  the  same  amount  of  energy  in  the  8-12. .tn  region  as  that  calculated  by 
summing  the  radiation  emitted  and  reflected  from  components  of  the  ship's  structure. 
The  reflected  component  is  from  surrounding  air  at  its  sea  level  temperature.  This 
effective  ship  temperature  was  given  by  equation  (38)  in  terms  of  the  actual  ship 
temperature,  T,  and  the  apparent  background  air  temperature,  Ta.  The  actual  ship 
temperature  is  calculated  by  the  Single  Element  Method  developed  in  Section  III. 
Calculated  effected  ship  temperatures  are  given  in  Table  3  for  the  years  1964-71 
and  the  3  seasons.  Table  3  lists  the  standard  probabilities.  For  example,  a 
temperature  of  25' C  at  a  probability  of  0.90  means  that  90  percent  of  all  calculated 
temperatures  are  25  C  or  lower.  Representative  effective  ship  temperatures  are 
given  in  Table  3  and  are  shown  in  Figure  7. 

In  this  and  in  future  graphs  and  tables,  representative  values  are  defined  to 
be  the  average  of  the  highest  and  the  lowest  values  of  8  yearly  temperatures  given 
in  the  tables.  Figure  7  indicates  that  the  effective  ship  temperatures  are 
generally  highest  in  the  summer  and  lowest  in  the  winter  as  might  be  expected. 
However  a  closer  examination  of  Table  3  shows  that  some  spring- fall  temperatures 
are  higher  than  their  corresponding  surrener  values.  This  occurs  at  higher  ship 
temperatures  which  correspond  to  mid-day  hours  where  spring-fall  solar  heating 
rates  are  the  larger.  Also,  the  fact  that  more  spring-fall  temperatures  are  not 
noticeably  higher  is  explained  by  an  analysis  of  the  weather  data  which  shows  that 
wind  speeds  are  about  15  percent  lower  in  the  summer  than  in  the  spring- fall.  High 
wind  speeds  are  especially  effective  in  reducing  the  high  ship  temperatures. 

&'*ckSLroijnc!  Temperatures 

The  background  temperature  was  defined  to  be  a  thermodynamic  measure  of  the 
sea  surface  radiance.  The  use  of  the  sea  surface  radiarce  as  the  background  is  a 
consequence  of  the  receiver-target  geometry  chosen  for  this  study,  i.e.,  a  ship  at 
or  near  the  hori2on.  This  background  was  defined  in  equation  (40)  as  int'insically 
derivable  from  a  mixture  of  sky  and  sea  radiances.  The  weather  data  show  that  sea 
temperature  is  always  within  a  few  degrees  centigrade  of  1 0" C .  But  sky  temperature 
which  depends  on  cloud  conditions  in  addition  to  air  temperature  varies  widely 
daily  and  will  also  depend  on  season  of  the  year.  Sky  temperature  calculation  is 
handicapped  by  the  availability  of  only  two  model  atmospheres,  one  for  summer  and 
one  for  winter.  This  causes  a  discontinuity  in  the  clear  sky  component  when  going 
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from  one  model  to  the  other.  That  is,  the  Midlatitude  Winter  Model  predicts  about 
-40  C3  and  the  Midlatitude  Summer  about  -20°C  for  the  clear  sky  component  of  tempera¬ 
ture.  Fortunately,  it  was  found  that  sky  temperatures  in  this  report  are  largely 
governed  by  its  cloudy  sky  component  because  of  the  perponderance  of  cloudy  skies 
encountered.  That  is,  an  analysis  of  the  weather  data  showed  that  median  cloud 
covers  are  90  percent  or  more  for  each  season.  Cloudy  sky  radiances  and  hence 
temperatures  are  only  slightly  influenced  by  choice  of  model  atmosphere  because 
sea  level  temperatures  are  taken  from  the  weather  data  and  air  temperatures  supplied 
by  the  model  are  used  for  low  altitudes  only  (median  cloud  heights  less  than  1  km). 
These  low  altitude  air  temperatures  do  not  vary  much  with  season.  Probabilities 
of  sky  temperatures  calculated  using  the  Weather  Ship  J  data  are  given  in  Table  4 
for  the  years  1964-71  and  the  3  seasons.  Representative  values  of  sky  temperature 
are  aiso  shown  in  Figure  8  for  each  season. 


Figure  8.  Probability  of  Observing  Reflected  Sky  Temperature 
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Table  4  and  Figure  8  show  that  summer  sky  temperatures  are  the  warmest  and 
winter  the  coldest  as  was  expected.  Most  of  the  variation  of  sky  temperature  with 
season  is  due  tc  differences  in  air  temperature  near  sea  level  because  of  the  low 
clouds.  This  difference  is  only  about  $°C  between  winter  and  summer  seasons.  How¬ 
ever,  the  lower  sky  temperatures  show  the  influence  of  some  relatively  clear  skies 
and  here  there  is  a  much  larger  difference  between  winter  and  summer  sky  temperature. 
The  spring-fall  season  contains  two  months  in  each  model  atmosphere  hence  it  has 
intermediate  sky  temperatures  in  the  low  temperature  region  which  is  influenced  by 
the  choice  of  model  atmosphere. 

Background  temperatures  are  computed  by  adding  the  sea  temperature  component 
to  the  sky  temperature.  These  were  calculated  for  a  representative  year  (1966)  and 
the  cumulative  statistics  are  shown  in  Figure  9  for  the  3  seasons.  Background 
temperatures  are  not  given  in  the  tables  because  their  calculation  was  not  printed 
in  the  ASIRCT  computer  program  output  when  the  8  years  of  weather  data  were  run. 

The  background  temperature  consists  of  approximately  1  part  sky  temperature  and  3 
parts  sea  temperature  for  viewing  one  ship  from  another.  Hence,  the  background 
temperatures  have  less  variability  than  the  sky  temperature  because  sea  temperatures 
vary  little  with  season.  The  fact  that  sea  temperatures  are  warmer  in  summer  con¬ 
tributes  to  the  larger  summer  background  temperatures. 
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Ship  Signatures  -  Contrast  Temperatures 

The  ship  contrast  temperature  was  defined  to  be  simply  the  effective  ship 
temperature  minus  the  background  temperature  Tsee  equation  (39)1.  A  zero  contrast 
temperature  is  associated  with  zero  radiance  difference  between  the  ship  structure 
and  the  combined  air-sea  background.  Contrast  temperatures  were  computed  by  the 
ASIRT  computer  program  for  daytime,  nighttime  and  total  time  of  operation. 

Daytime  Contrast  Temperatures 

Daytime  contrast  temperatures  were  computed  for  all  daylight  hours  and  for  the 
mid-day  hours  between  10  a.m.  and  3  p.m.  Cumulative  probabilities  of  daytime  con¬ 
trast  temperatures  are  given  in  Table  5  for  the  years  1964-71  and  for  the  three 
seasons.  Representative  values  of  daytime  contrast  temperature  which  appear  in 
Table  5  are  also  plotted  in  Figure  10  for  each  season. 
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Figure  10.  Probability  of  Observing  Daytime  Contrast  Temperatures 


Figure  10  and  Table  5  show  that  statistically  the  differences  between  winter 
and  sunmer  daytime  contrast  temperatures  are  not  large  if  the  contrast  tempera¬ 
tures  are  less  than  about  10°C.  This  closeness  can  be  explained  by  comparing  the 
ship  temperatures  plotted  in  Figure  9.  There,  it  is  seen  that  the  difference 
between  summer  and  winter  temperatures  is  about  the  same  (within  0.5cC)  for  80%  of 
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the  time.  Thus,  the  near  agreement  of  suirner  and  winter  contrast  temperatures  is 
largely  a  fortuitous  combination  of  ^'-.p  and  background  temperatures.  It  will  be 
shown  that  the  agreement  in  contrast  temperature  is  caused  by  similar  early  morning 
and  later  afternoon  values  where  solar  heating  effects  are  minimized  due  to  the  un¬ 
favorable  angular  position  of  the  sun.  The  spring-fall  daytime  contrast  tempera¬ 
tures  are  significantly  larger  than  the  other  seasons  due  to  the  relatively  larger 
ship  temperatures  and  relatively  low  background  temperatures.  That  is,  the  solar 
heating  is  maximized  as  a  consequence  of  favorable  solar  altitude  angles  while  the 
background  radiance  and  therefore  the  effective  background  temperature  is  reduced 
by  the  smaller  probability  of  cloud  cover.  High  spring-fall  and  summer  contrast 
temperatures  are  evident  in  the  mid-day  calculations.  These  are  given  in  Tables 
6-11  for  the  hours  10  a.m.  to  3  p.m.  ,  the  years  1964-71  and  the  3  seasons.  Repre¬ 
sentative  average  values  are  given  in  each  table  and  these  are  plotted  in  Figure  11 
for  the  maximum  daytime  contrast  temperatures  and  the  3  seasons.  The  maximum 
temperatures  are  found  to  occur  at  1  p.m. 


Figure  11.  Probability  of  Observing  Maximum  Contrast  Temperatures 


Figure  11  and  Tables  6-11  show  that  the  mid-day  contrast  temperatures  are  the 
highest  in  the  spring- fall  season  due  primarily  to  the  more  favorable  solar  alti¬ 
tude  angles.  However,  this  data  also  shows  th,.  the  mid-day  summer  contrast 
temperatures  are  significantly  higher  than  the  corresponding  winter  values.  These 
higher  mid-day  contrasts  of  summer  are  not  obvious  in  Figure  10  which  showed  similar 
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winter  and  summer  overall  daytime  contrasts.  Hence,  it  can  be  concluded  that  it  is 
the  early  morning  and  later  afternoon  contrast  temperatures  which  are  similar  in 
winter  and  summer  The  mid-day  contrast  temperatures  are  the  highest,  temperatures 
occuring  in  Figure  10  and  these  are  seen  to  be  higher  in  summer.  Mid-day  contrasts 
are  significantly  larger  in  sutnner  than  winter  due  to  the  relatively  larger  ship 
temperatures.  These  can  be  seen  in  Figure  7  which  shows  that  the  difference 
between  summer  and  winter  ship  temperatures  increases  dramatically  for  probabilities 
greater  than  0.8. 


Nighttime  Cont rast  Temperatures 

Cumulative  probability  of  contrast  temperature  for  nighttime  operation  are 
given  in  Table  12  for  the  years  1964-71  and  the  3  seasons.  Representative  values 
of  nighttime  contrast  temperatures  appear  in  Table  12  and  are  also  plotted  in  Figure 
12  for  the  3  seasons. 


Figure  12.  Probability  of  Observing  Nightt'.me  Contrast  Temperatures 

Figure  12  shows  that  nighttime  contrast  temperatures  are  about  a  factor  of  2 
lower  than  daytime  values  and  have  less  variability  between  seasons.  These  facts 
show  the  importance  of  solar  heating  in  enhancing  daytime  FLIR  performance.  It  may 
be  surprising  that  summer  nighttime  contrasts  are  generally  the  lowest.  This  be¬ 
havior  is  due  to  the  relatively  high  summer  background  temperatures  without  off¬ 
setting  high  ship  temperatures. 
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TABLE  3  SHIP  TEMPERATURE  FOR  SHIP  J(1964-71) 


PROS 

1964 

1965 

1966 

1967 

1968 

1969 

1970 

1971 

REP 

.96 

17.28 

15.99 

19.49 

17.91 

18.98 

17.27 

19.38 

18.98 

17.74 

.90 

14.97 

14.30 

16.23 

15.13 

15.54 

14.54 

15.71 

15.48 

16.27 

.80 

13.74 

12.90 

13.61 

13.42 

13.36 

12.82 

13.c6 

13.92 

13.37 

.70 

12.98 

12.32 

12.71 

12.72 

12.61 

11.92 

12.33 

13.01 

12.47 

.60 

12.44 

11.74 

12.02 

12  24 

12.00 

11.19 

11.58 

12.41 

11.82 

.50 

11.91 

11.18 

11.51 

11  73 

11.48 

10.54 

10.93 

11.82 

11.23 

.40 

11.43 

10.66 

11.01 

11.19 

10.96 

9,90 

10.36 

11.27 

10.52 

.30 

10.92 

10.17 

10.45 

10.62 

10.46 

9.17 

9.76 

10.64 

10.05 

.20 

10.14 

9.66 

9.86 

10.04 

9.93 

8.36 

9,15 

9.83 

9.25 

.10 

9,50 

8.97 

9  04 

9.22 

8.87 

7.31 

8,30 

8.84 

8.41 

WINTER  (JAN..  FEB.,  NOV..  DEC.) 


PROB 

1964 

1965 

1966 

1967 

1968 

1969 

1970 

1971 

REP 

.95 

28.15 

30.54 

28.35 

25.87 

28.08 

26.92 

24.86 

28.44 

27.70 

.90 

22.81 

24.82 

22.68 

21.21 

22.24 

22.45 

20.52 

23.06 

22.67 

.80 

17.88 

18.26 

17.64 

16  45 

17  02 

17.61 

16.93 

17.84 

17.36 

.70 

16.18 

15.70 

15.66 

14  92 

15.05 

15.52 

15.40 

16.70 

15.81 

.60 

14.85 

14.72 

14.62 

13.73 

14.42 

14.48 

14.48 

15.62 

14.68 

50 

13.82 

13.95 

13.60 

12.89 

13.72 

13.69 

13.58 

14.61 

13.75 

.40 

12.74 

13.10 

12.68 

12  28 

12  78 

12.96 

12.60 

13.45 

12.87 

.30 

11.76 

12.15 

11.87 

11.65 

11.76 

12.02 

11.68 

12.25 

12.01 

.20 

10.95 

11.34 

11.15 

11  02 

10.66 

11.34 

10.79 

11.36 

11.01 

.10 

9.95 

10.25 

10.32 

9.69 

9.32 

10.47 

9.70 

10.25 

9.90 

SPRING-FALL  (MAR..  APR.,  SEPT.,  OCT) 


PROB 

1964 

1965 

1366 

1967 

1968 

1969 

1970 

1971 

REP 

.95 

26.63 

28.38 

28.82 

27.09 

28.41 

26.09 

25.70 

26.68 

27.26 

.90 

22.56 

23.53 

23.33 

22  66 

23  47 

22.14 

21.77 

22.67 

22.65 

.80 

18.12 

19.49 

18.89 

18.41 

19.28 

17.85 

17.83 

18.75 

18.66 

.70 

16.91 

17.38 

17.07 

16.55 

17.46 

15.86 

16.51 

17.11 

16.62 

.60 

15.88 

16.13 

15.89 

15.53 

18.49 

15.13 

15.50 

15.96 

15.81 

.50 

15.30 

15.55 

15.14 

14.87 

15  51 

14.55 

14.72 

15.37 

15.06 

.40 

14.67 

15.06 

14.41 

14.34 

14.44 

14.01 

14.02 

14.66 

14.54 

.30 

13.95 

14.34 

13.61 

13.74 

13.46 

13.42 

13.34 

13.77 

13.6% 

.20 

12.87 

13.53 

12.70 

13.02 

1?  65 

12.74 

12.58 

12.78 

13.06 

.10 

11.58 

12.58 

11.62 

11.31 

11.77 

11.71 

11.54 

11.77 

11.96 

SUMMER  (MAY,  JUNE.  JULY,  AUG.) 
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J  u&ta  exhibit  large  amounts  of  cloud  cf  er  which  effects  both  the  daytime  ship 
temperatures  and  all  background  temper,  jres.  Hence,  these  contrast  temperatures 
may  show  mere  variation  with  season  in  a  more  cloud-free  location. 

It  was  previously  stated  that  target  ship  heading  controls  the  contrast 
temperature  of  a  ship  target  because  its  solar  heating  depends  on  its  heading. 
Further,  figure  6  defined  3  basic  headings  depending  on  whether  the  vertical  eluent 
receives  direct  solar  heat  all,  one-half  or  none  of  the  time.  The  ASIRCT  computer 
program  was  run  for  the  3  basic  headings  and  a  representative  year  (1966),  It  was 
found  that  the  daytime  contrast  temperatures  are  very  sensitive  to  heading  due  to 
the  great  differences  in  incident  solar  heating  as  shown  in  Figure  6.  These  are 
shown  in  Figure  14  for  the  3  seasons.  Figure  14  shows  that  there  is  about  a  factor 
of  2  difference  in  daytime  contrast  temperature  between  a  due  east  (90°)  and  a  due 
west  (270  )  heading  for  the  starboard  vertical  element.  Also,  the  coincidence  of 
some  winter  and  summer  daytime  contrast  temperature  observed  in  Figure  10  for  the 
due  east  heading  is  not  present  in  the  remaining  headings.  Daytime  contrast  tem¬ 
peratures  for  the  due  south  and  due  west  headings  show  spring- fall  values  are  the 
largest  in  agreement  with  the  results  found  for  the  due  east  heading.  This  result 
remains  true  due  to  the  relatively  lower  background  temperatures  calculated  in  the 
spring- fall  season  (see  Figure  9).  The  maximum  (1  p.m.)  contrast  temperatures  show 
an  even  greater  dependence  on  heading  than  the  overall  daytime  contrasts.  Here, 
tne  due  south  heading  has  maximum  contrast  temperatures  which  are  about  40%  of  the 
due  east  heading  when  probabilities  of  observing  contrast  temperatures  or  less  are 
compared.  This  drops  to  about  30%  for  comparing  due  west  to  the  due  east  heading. 
The  night  contrast  temperatures  are  nearly  the  same  for  the  various  headings 
although  the  due  east  contrasts  are  slightly  higher  (e.g.,  0.25-0.50  C)  due  to  the 
lasting  effect  of  solar  heating.  An  interesting  observation  is  that  the  daytime 
contrast  temperatu-es  for  the  due  west  heading  are  nearly  equal  to  the  nighttime 
contrast  temperatures  seen  in  Figure  12  and  Table  12.  Thus,  the  relatively  low, 
totally  diffuse  so’ar  energy  received  in  this  heading  and  at  t.iis  location  is  almost 
totally  ineffective  in  building  daytime  contrast  temperatures. 

Si gnatures  with  Air  Backgrounds 

The  previously  given  contrast  temperatures  calculations  assume  that  the  back¬ 
ground  of  the  target  ship  consist  of  radiant  emission  from  the  sea  and  reflected 
radiance  from  the  sky  and  off  of  the  sea  surface  wave  slopes.  In  this  background 
r-del ,  the  target  must  be  s'ightlv  before  the  horizon  so  that  no  radiance  from  an 
air  path  along  the  1 ine-of-sight  to  the  target  but  behind  it  appears  in  the  back¬ 
ground.  If  part  of  the  ship  is  above  the  horizon,  radiance  from  this  air  path  will 
contribute  to  tne  background.  Contrast  temperatures  were  calculated  for  backgrounds 
consisting  of  sea  level  air  temperature  only.  This  is  the  limiting  case  where  the 
target  ship  is  at  or  slightly  over  horizon  and  no  sea  emission  or  sky  ••eflection 
reaches  the  infrared  detector.  In  the  general  case  of  a  ship  near  the  horizon,  the 
background  radiance  will  be  generally  less  than  the  all  air  emission  but  greater 
than  the  all  sea  emission  plus  sky  reflection  case.  This  is  due  to  the  fact  that 
the  apparent  temperature  of  the  sky  is  sometimes  quite  cold  (see  Figure  3)  and  is 
always  less  than  the  sea  level  air  temperature. 

Contrast  temperatures  for  backgrounds  consisting  of  sea  level  air  are  given  in 
Figure  15.  Figure  15  presents  the  night,  day  and  total  contrast  temperature  for 
each  of  the  three  seasons  and  weather  data  taken  by  Ship  J  in  1966.  This  figur° 
shows  much  smaller  ship  contrast  temperatures  in  comparison  with  those  having  a  sea 
surface  background  (see  Figures  10-13).  flight  contrast  temperatures  are  especially 
reduced  and  do  not  vary  significantly  with  season.  Nighttime  contrasts  are  low 


Figure  and  Table  13  show  that  the  total  contrast  temperatures  do  not  vary 
significantly  with  season.  This  is  understandabl e  from  the  previous  results  where 
it  was  shown  that  only  the  mid-day  contrast  temperatures  show  large  seaonsal  dif¬ 


ferences  due  primarily  to  differences 
day  hours  occur  only  257-  of  each  day, 
distributions  of  Figure  13.  However, 


in  received  solar  heating.  Since  the  mid- 
their  effect  is  largely  lost  in  the  cumulative 
it  should  be  remembered  that  the  Weather  Ship 
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SUMMER  DAY  CONTRAST  TEMPERATURES.  C 


SPRING-FALL  DAY  CONTRAST  TEMPERATURES.  "C 


Figure  14.  Representative  Probahil  it  ies  of  Observing  Daytime  Contrast 
Temperatures  for  Various  Ship  Headings 
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because  they  are  governed  by  the  difference  between  ambient  air  and  background 
temperatures.  That  is,  the  ship  temperatures  increase  due  to  internal  sources 
generally  contribute  only  slightly  to  night  contrast  temperatures  (see  Appendix  A). 
However,  in  the  present  case  of  air  backgrounds,  internal  sources  are  a  major 
contributor  to  night  contrast  temperature.  Day  contrasts  are  somewhat  reduced  for 
the  air  background  but  spring-fall  values  are  the  highest  as  they  were  for  a  sea 
surface  background  (see  Figure  10).  This  also  occurs  because  spring-fall  ship 
temperature  are  about  the  same  as  those  of  summer  (see  Table  3)  but  background 
temperatures  are  lower  due  to  the  colder  season.  The  total  contrast  temperatures 
of  Figure  15  show  lower  winter  values  due  to  the  much  lower  winter  ship  temperatures. 
This  behavior  did  not  occur  with  sea  surface  backgrounds  because  these  have 
relatively  lower  background  temperatures  than  the  present  case  (see  Figure  9). 

Single  Element  Method 

The  results  of  contrast  temperature  calculations  presented  in  Figures  7-13  and 
Tables  3-13  were  based  on  the  Patrol  Frigate  Model  discussed  in  Appendix  A.  This 
model  treats  the  ship  as  a  single  vertical  element,  and  applies  three  correction 
factors.  These  corrections  account  for:  (1)  variations  in  ship  internal  tempera¬ 
tures;  (?)  variations  in  thermal  capacities  of  ship  construction;  and  (3)  the 
presence  of  one  or  more  exhaust  stacks.  It  was  shown  in  Appendix  A  that  the 
individual  correction  factors  are  likely  to  be  small  for  a  ship  having  a  cooled 
stack.  The  implication  of  the  correction  factors  being  small  is  that  the  features 
calculated  using  the  Patrol  Frigate  Model  would  approximately  apply  to  any  similarly 
constructed  ship.  To  judge  the  influence  of  ship  model  on  contrast  temperature, 
the  3  correction  factors  used  in  the  Single  Element  Method  were  sunned  and 
statistically  analyzed.  That  is,  the  probability  of  occurrence  of  the  sum  was 
computed  for  each  of  the  seasons  in  the  1964-71  time  period.  These  were  found  to 
be  remarkably  similar  for  each  season  and  for  each  year.  The  1966  spring-fall 
season  was  -*hosen  as  representative  and  the  sum  of  the  correction  factors  is 
plotted  in  Figure  16  for  this  season. 

Figure  16  shows  that  the  total  temperature  correction  factor  is  rarely  above 
1 . 0 : C  and  has  a  median  value  of  almost  0.5  C.  An  examination  of  the  ship  tempera¬ 
tures  given  in  Tabic  3  for  8  years  of  weather  data  shows  that  the  variation  In  ship 
temperature  between  years  is  comonly  greater  than  1.0WC.  Furthermore,  the  tem¬ 
perature  additions  shown  in  Figure  15  are  small  compared  to  the  total  temperature 
rises  above  ambient  (8-10cC)  shown  in  Table  3.  Hence,  it  is  concluded  that  the 
choice  of  ship  model  is  not  critical  and  that  the  results  predicted  herein  apply 
to  all  similarly  constructed  ships.  However,  note  that  these  results  use  a  cooled 
stack  and  this  conclusion  would  not  hold  for  a  model  containing  a  hot  stack. 

The  work  cf  this  report  is  directed  towards  finding  the  target  contrast  sig¬ 
nature  that  exists  at  the  ship  target.  This  ship  target  signature  is  useful  in 
system  studies  which  define  system  performance  from  a  target  contrast  temperature 
such  as  those  performed  for  the  StAFIRE  Program.  However,  for  many  other  purposes 
the  signal  power  definition  of  ship  signature  is  needed.  The  signal  power  is 
defined  in  equation  (1)  and  this  yields  the  radiant  power  received  which  is  useful 
in  providing  a  signal  to  an  infrared  system.  A  consideration  of  the  apparent 
signal  requires  that  the  attentuation  of  radiance  and  the  emission  over  air  path 
between  the  detector  and  the  target  be  known.  Attenuation  due  to  the  response  of 
the  detector  must  also  be  known.  Both  of  these  are  outside  the  scope  of  the  present 
work.  However,  it  is  possible  to  relate  the  contrast  temperature  of  this  work  to 
the  radiant  contrast  which  exists  at  the  target.  Target  radiant  contrast  is  simply 
the  difference  between  target  surface  radiance  and  background  radiance  integrated 
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TABLE  4  SKY  TEMPERATURE  FOR  SHIP  JI1964  71) 


PROB 

1964 

1965 

1966 

1967 

1968 

1969 

1970 

1971 

REP 

.96 

10.76 

10.06 

10.24 

10.55 

10.90 

9.48 

”la56~ 

10.81 

10.15 

.90 

9.88 

8.83 

9  13 

9.82 

9.98 

8.18 

9.43 

9.78 

9.03 

.80 

8,39 

7.22 

7.62 

7.99 

7.61 

5.97 

6.89 

8.33 

7.81 

.70 

7.02 

6  78 

6.06 

6.38 

5.35 

4.03 

4.65 

6.64 

5.53 

.60 

5.52 

3.80 

4  58 

4  60 

3.35 

1.90 

2.76 

5.00 

3.71 

.50 

3.63 

2  11 

2  80 

2.37 

1.39 

-  .07 

.81 

2.90 

1.78 

.40 

1  11 

-  93 

70 

.07 

-  1.54 

-  3.11 

-  2.30 

.49 

-  1.00 

.30 

-  4.26 

-  6  94 

-  3  19 

-  4.68 

-  7.15 

-  7  73 

-  8.67 

-  3.09 

-  5.88 

.20 

-12  71 

-13  29 

-  9  91 

-11.5? 

-14.61 

-14.39 

-15.52 

-  8.92 

-12.22 

10 

-22.83 

-22  96 

-18,26 

-20.72 

-23.66 

-22.71 

-28  40 

-18.34 

-23.33 

WINTER  (JAN..  FEB  .  NOV..  DEC.) 


PROB 

1964 

1965 

1966 

1967 

1968 

1969 

1970 

1971 

REP 

.95 

13.44 

12  14 

11  75 

11.80 

11  70 

12.64 

12.89 

14.44 

13.10 

.90 

11,88 

10  83 

10.38 

10.42 

10.49 

11.73 

11.55 

13.52 

11.95 

.80 

9.46 

9.55 

8.87 

8.67 

8.62 

9.64 

9.68 

10.48 

9.55 

70 

7  69 

7  95 

7.50 

6.97 

6.78 

8.25 

7.62 

8.37 

7.58 

60 

5  92 

6.04 

5  84 

5  44 

4.93 

6.41 

5  77 

6.17 

6.67 

.50 

3.37 

3  84 

3  92 

3  76 

3.06 

4  12 

3.26 

3.62 

3.59 

40 

47 

1.53 

1  01 

1  71 

79 

1  09 

43 

1.04 

1.08 

30 

»  3.55 

-  2  75 

-  2  85 

-  1  56 

-  2  85 

-  3  53 

-  3.49 

-  2.52 

-  2.56 

.20 

-  9.87 

-  7  85 

8  48 

-  6  57 

6.82 

-  9  46 

-  995 

-  7.69 

-  8.26 

.10 

-17  85 

-14  41 

-15  50 

-15  22 

-15.03 

-18.78 

-17.44 

-14.58 

-16.60 

SPRING  FAIL  (MAR  .  APR  .  SEPT  .  OCT  I 


PROB 

1964 

1965 

1966 

1967 

1968 

1969 

1970 

1971 

REP 

.95 

13.64 

13  15 

13.40 

13  31 

13  20 

13  00 

14.21 

14.16 

13.61 

90 

12.85 

12.38 

1266 

12  47 

11  66 

12.47 

12.86 

13.13 

12.40 

80 

11.57 

11  08 

11.51 

10.98 

10.07 

11.26 

11  56 

11.39 

10.82 

.70 

10  18 

958 

10  17 

954 

8.65 

9  61 

10.33 

1005 

9.49 

.60 

9  05 

8  01 

8  89 

7  96 

7.33 

7  88 

9.31 

8.95 

8.32 

50 

7  79 

6  52 

7  b8 

6  76 

5.96 

6.46 

0.14 

7.37 

7.05 

40 

6  33 

4.66 

6  15 

5  25 

4.39 

4.64 

6  87 

5.52 

5.63 

.30 

4  18 

1.66 

3  87 

3.44 

1.23 

2.07 

5.27 

2.88 

3.25 

20 

-  1.56 

-  2  40 

-  1  29 

-  76 

-  3.02 

-  3.38 

2.43 

-  .74 

-  .48 

10 

-  7.84 

-  7  72 

-  7  86 

-  6.71 

-  8.07 

-  8.78 

-  4.70 

-  6.09 

-  6.74 

SUMMER  (MAY.  JUNE,  JULY.  AUG.) 
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TABLE  5  DAY  CONTRAST  TEMPERATURE  FOR  SHIP  Jtl964-71) 


PROB 

1964 

1965_ 

JJ966 

1967 

1968__ 

1969  „ 

1970 

1971 

jEP 

.95 

16.87 

13.52 

18.78 

17.56 

19.23 

18.48 

18.71 

18.51 

16.16 

.90 

12.45 

10.41 

14.24 

13.66 

15.31 

13.15 

14.46 

13.01 

12.86 

.80 

8.34 

7.69 

10.04 

9.18 

10.76 

8.54 

10.55 

9.09 

9.12 

.70 

6.29 

5.89 

7.58 

7.27 

7.84 

6.23 

8.01 

6.86 

6.95 

.60 

4.97 

4.92 

5.86 

5.87 

6.10 

4.94 

5.92 

5.38 

6.61 

.50 

3.81 

3.69 

4.30 

4.44 

4.75 

3.81 

4.54 

3.94 

4.22 

.40 

2.86 

2.89 

3.03 

3.44 

3.47 

2.93 

3.36 

2.90 

3.18 

.30 

2  07 

2.10 

2.17 

2.67 

2.26 

1.98 

2.41 

2.19 

2.33 

.20 

1.36 

1.33 

1.50 

1.96 

1.41 

1.06 

1.45 

1.32 

1.61 

10 

.57 

.39 

.39 

1.01 

.33 

.01 

.28 

.31 

.52 

WINTER  UAN.,  FEB.,  NOV.,  DEC.) 


PROB 

1964 

1965 

1966 

1967 

1968 

1969 

1970 

1971 

REP 

95 

26.85 

27.95 

25.82 

24  41 

24.83 

24  94 

22.62 

26.79 

25.29 

.90 

20.60 

22.32 

20.57 

17.94 

20.03 

19.78 

17.24 

20.44 

19.78 

.80 

14.04 

15.91 

14.00 

12.46 

13.43 

14.04 

11.77 

13.97 

14.19 

.70 

10.44 

11.53 

10.19 

9.53 

9.51 

10.51 

9.18 

10.30 

10.36 

.60 

7  73 

8.57 

7.91 

7.11 

6.87 

8.17 

7.30 

7.56 

7.72 

50 

5,74 

6.21 

5.83 

5  16 

5,12 

6.35 

5.34 

5.60 

6.74 

.40 

3.88 

4  42 

429 

3.64 

3.80 

4.69 

3.96 

4.03 

4.17 

.30 

2.71 

3.07 

3.04 

260 

2  72 

3.39 

2.89 

3.00 

3.00 

.20 

1.81 

1.99 

1  85 

1.90 

1.81 

2.15 

1.78 

2.12 

1.98 

.10 

.95 

99 

73 

.91 

.85 

1.20 

.56 

1.13 

.88 

SPRING  FALL  (MAR..  APR..  SEPT.,  OCT.) 


PROB 

1964 

1965 

1966 

1967 

1968 

1969 

1970 

1971 

REP 

.95 

19.67 

21.90 

22.78 

20.75 

20.86 

19  86 

17.41 

18.91 

20.10 

.90 

14.61 

16.84 

16.51 

15.38 

15.86 

1541 

13.15 

14.62 

15.00 

.80 

9.86 

11.04 

10.51 

10.34 

10.72 

iy.t>2 

8  71 

9.52 

9.72 

.70 

7.29 

8.18 

741 

7.67 

8.16 

7.88 

6.06 

7.14 

7.12 

.60 

5.44 

6.30 

5.39 

5.74 

6.35 

5.89 

4.49 

6.47 

6.87 

.50 

4.10 

4.82 

4.04 

4.20 

4.71 

4.33 

3.37 

4.03 

4.10 

.40 

3.22 

3.66 

3.14 

3.12 

3.58 

3.26 

2.70 

3.16 

3.18 

.30 

2.53 

2.74 

2.41 

2.38 

2.67 

2.46 

2.11 

2.42 

2.43 

.20 

1.94 

2.05 

1.79 

1.81 

1.96 

1.93 

1.54 

1.72 

1.89 

.10 

1.32 

1.33 

1.12 

1.17 

1.2C 

1.19 

.73 

.86 

1.03 
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TABLE  6  10  A.M.  CONTRAST  TEMPERATURE  FOR  SHIP  JU964-71) 


PROS 

19G4 

1965 

1966 

1967 

1968 

1969 

1970 

1971 

REP 

.95 

11.48 

9.66 

11.25 

13.75 

14.83 

12.20 

14.10 

13.30 

12.25 

.90 

9.39 

8.33 

9.19 

11.00 

12.43 

7.60 

13.12 

10.10 

10.73 

.80 

7.40 

5.77 

6.90 

7.50 

8.30 

5.85 

8.60 

6.45 

7.19 

.70 

5.09 

4.59 

5.81 

5.50 

5  49 

5.01 

6.13 

5.66 

5.36 

.60 

4.06 

3.70 

4.60 

4.60 

4.68 

3.92 

4.31 

4.71 

4.15 

.50 

3.35 

3.17 

3  70 

3.76 

4.16 

3.25 

3.46 

3.64 

3.71 

.40 

3.04 

2.25 

2.48 

3.33 

3  38 

2.74 

2.73 

2.72 

2  82 

.30 

2.53 

1.27 

2.01 

2.60 

2.36 

2.13 

2.16 

2.16 

1.94 

.20 

1.81 

.54 

1.07 

2.00 

1.49 

1.22 

1.54 

1.32 

1.27 

.10 

.81 

-  06 

10 

60 

87 

.36 

.49 

.45 

.36 

WINTER  ( 

JAN..  FE8..  NOV 

.  DEC.) 

PROS 

1964 

1965 

1966 

1967 

1968 

1969 

1970 

1971 

REP 

.95 

20.85 

25.35 

20.55 

21.81 

19.35 

21.41 

14.90 

18.75 

20.13 

.90 

16.70 

17.90 

18.10 

14.90 

14.80 

18.08 

13.20 

16.26 

15.65 

.80 

11.60 

13.40 

12.80 

11  13 

11.90 

14.55 

10.95 

13.50 

12.75 

70 

9.03 

11  35 

10  40 

9.52 

8  70 

10.33 

9.22 

9.68 

10.03 

.60 

6.98 

8.6f 

8.98 

8.20 

6.87 

8.30 

6.43 

8.30 

7.71 

.50 

6.05 

6.96 

7.80 

6.38 

5.88 

6.95 

5.29 

6.45 

6.56 

.40 

4.24 

5.57 

6.39 

5  28 

4  76 

5.43 

4.34 

5.38 

6.31 

.30 

3.39 

3.52 

5  40 

3.38 

3.66 

4.54 

3  37 

3.90 

4.39 

.20 

1.64 

1.82 

3  32 

1.92 

2.44 

2.28 

1.92 

2.15 

2.48 

.10 

1.03 

1.02 

1.68 

73 

1  26 

1.19 

.96 

1.30 

1.21 

PROB 

1964 

SPRING  FALL  (MAR., 

1965  1966  1967 

.  APR  .  SEPT  .  OCT.) 

1968  1969  1970 

1971 

REP 

.95 

18.80 

18.45 

21  60 

16.40 

19.76 

17.7/ 

14.45 

16.50 

18.03 

.90 

14.60 

16.53 

16.80 

14.40 

18  09 

16.64 

10.98 

13.00 

14.54 

.80 

11.64 

12.20 

12.20 

10.92 

12.80 

12.95 

8.13 

10.00 

10.17 

.70 

9.20 

9.97 

7.80 

8.68 

9.77 

10.18 

6.62 

8.50 

8.90 

.60 

6.85 

7.93 

6.55 

6  87 

8.63 

8.10 

5.33 

6.83 

6.72 

.50 

5.40 

6.67 

4.80 

5  75 

6.00 

6.54 

4.40 

5.90 

6.54 

.40 

4.46 

5.45 

3.66 

4.17 

4.36 

5.54 

3.68 

5.08 

4.60 

.30 

3.42 

4.22 

2.81 

2.92 

3.34 

3,84 

2.82 

4.07 

3.52 

.20 

2  16 

2.74 

1.96 

2.10 

2.44 

242 

1.96 

2.35 

2.35 

.10 

1.49 

1.37 

1.23 

1.41 

1.36 

1.56 

1.21 

1.20 

1.38 

S 

t 

I 
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TABLE  7  11  A.M.  CONTRAST  TEMPERATURE  FOR  SHIP  J(1964  71) 


PROB 

1964 

1965 

1968 

1967 

1968 

1969 

1970 

1971 

REP 

.95 

17.40 

11.95 

15.40 

15.65 

17.63 

13  94 

19.95 

17.75 

15.95 

.90 

13.90 

8.45 

12.96 

12.30 

13.97 

10.90 

16.15 

12.50 

12.30 

.80 

9.40 

6.63 

11.12 

8.87 

10.4/ 

8.75 

11.80 

9.13 

9.22 

.70 

6.85 

5.48 

8.93 

7.47 

8.67 

6.44 

8.15 

8.25 

7.21 

.60 

5.28 

4.96 

7  62 

6.07 

7.10 

5.58 

6.53 

6.63 

6.29 

50 

4.20 

4.10 

6.58 

5.13 

5.92 

4.55 

5.19 

4.67 

5.34 

40 

3.18 

3.24 

4.52 

3.92 

4.63 

3.67 

3.98 

3.13 

3.83 

.30 

2.06 

2.29 

2.46 

2.68 

3.06 

2.45 

2.82 

2.25 

2.56 

.20 

1.28 

1.62 

1.87 

2.04 

1.71 

1.27 

1.34 

1.16 

1.60 

.10 

.52 

.22 

48 

1.24 

.42 

.30 

-  .35 

.30 

.45 

PROB 

1964 

WINTER  (JAN.,  FEB.,  NOV 

1965  1966  1967  1968 

,  DEC.) 

1969 

1970 

1971 

REP 

.95 

33.75 

27.60 

26.30 

26.25 

28.65 

26.65 

26.10 

25.98 

29.87 

.90 

25.50 

24.00 

22.20 

19.50 

22.55 

22.73 

22.08 

23.40 

22.50 

.80 

18.75 

21.00 

16.08 

14.33 

13.10 

16.80 

17.03 

16.20 

17.67 

70 

14.00 

17.20 

12,80 

12.00 

15.83 

12.30 

13.60 

12.93 

14.60 

.60 

12  40 

14.50 

10  07 

10.25 

11.90 

10.24 

9.90 

10.93 

12.20 

.50 

10.20 

11.25 

9.25 

8.33 

9.92 

8.65 

8.70 

S.75 

9.79 

4C 

8.50 

8.88 

690 

5.17 

7.65 

6.13 

6.97 

8.15 

7.03 

.30 

5  75 

5.50 

5  10 

2  75 

5.95 

4.48 

4.81 

5.30 

4.35 

.20 

2  60 

2.20 

2.28 

1.80 

3.32 

3.08 

3.48 

3.91 

2.86 

10 

1.30 

1  33 

1  42 

.80 

1.53 

1.47 

1.24 

1.64 

1.22 

PR08 

1964 

SPRING  FALL  (MAR.. 

1965  1966  1967 

APR  .  SEPT..  OCT. 

1968  1969 

) 

1970 

1571 

REP 

.95 

23.33 

25.95 

26.85 

20.48 

22.65 

22.50 

17.13 

20.38 

21.99 

90 

1/.55 

22.92 

22.28 

17.07 

19.70 

19.50 

15.50 

18.47 

19.21 

80 

12.55 

16.27 

17  70 

13.55 

14.63 

16.57 

11.75 

14.73 

14.73 

.70 

10  33 

13  93 

13  03 

11.35 

12.04 

13.00 

9.75 

11.06 

11.84 

.60 

8.07 

11,95 

10.45 

9.64 

8.70 

950 

7.42 

8.40 

9.69 

50 

658 

9  17 

8.58 

7.80 

7.10 

7.50 

6.25 

7.44 

7.71 

.40 

5  36 

7.40 

6.65 

6.60 

6.28 

6.50 

5.00 

5.83 

6.20 

.30 

4.16 

5.60 

4.98 

4.68 

4.96 

3.30 

3.33 

3.57 

4,45 

.20 

2.79 

3.46 

3.24 

2.36 

2.64 

2.20 

2.04 

2.58 

2.75 

.10 

1.79 

1.85 

1.49 

1.88 

1.66 

1.50 

.75 

1.09 

1.32 
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TA8LE8  NOON  CONTRAST  TEMPERATURE  FOR  SHIP  J(1964-71) 


PROB 

1964 

1965 

1966 

1967 

1968 

1969 

1970 

1971 

REP 

.96 

18.00 

15.62 

19.01 

18.90 

23.33 

16.70 

23.03 

26.55 

21.09 

.90 

13.60 

12.85 

16.90 

14  80 

18  10 

13  80 

1935 

19.40 

16.13 

.80 

11.00 

9  07 

13.68 

11.40 

13.44 

10.93 

14.20 

12.95 

11.64 

.70 

8.20 

7.89 

11.03 

9.13 

12.04 

8.60 

11.62 

10.73 

9.97 

.60 

6.38 

6  48 

8.73 

7.43 

10.20 

7.13 

8.63 

7.20 

8.29 

.50 

5.26 

5  39 

7.19 

6.30 

7  44 

5.67 

6.75 

5.75 

6.35 

.40 

4  60 

*5  52 

5.90 

5.22 

5  83 

4.68 

4.86 

4.12 

5.01 

.30 

3.20 

3  25 

3.98 

3.32 

4.54 

3.85 

3.95 

3.24 

3.87 

.20 

1  90 

2,31 

2.12 

2  48 

2  16 

2.26 

1.84 

2.15 

2  16 

.10 

90 

82 

,86 

1  62 

79 

.36 

33 

.57 

.98 

WINTER  (JAN  .  FEB..  NOV..  DEC.i 


PROB 

1964 

1965 

1966 

1967 

1968 

1969 

1970 

1971 

REP 

.95 

39.75 

34.46 

33  25 

28.95 

34.67 

33.75 

32.17 

29.45 

34.36 

90 

32  83 

30  58 

27  96 

23  85 

28  60 

25.50 

23.76 

25.43 

28.20 

.80 

19.70 

22  65 

22  10 

17.52 

21  72 

21.00 

15.20 

20.80 

18.93 

.70 

16.78 

19.16 

17.60 

14.57 

1760 

17.17 

13.20 

14.93 

16.18 

60 

12.70 

15.70 

12.47 

11  32 

12.70 

13  50 

10.68 

12.90 

13.19 

.50 

9  13 

1208 

10.50 

9,25 

9  33 

11  13 

9.17 

10.50 

10.63 

.40 

7.22 

7  94 

8.90 

6  70 

7  60 

9.25 

8.05 

8.57 

7.80 

.30 

4  56 

5  43 

5  80 

4  06 

4  63 

7  69 

6  05 

6.06 

5.88 

20 

2  38 

3  52 

3.91 

2.42 

2  76 

4  60 

2.56 

3.52 

3.58 

.10 

1.46 

1  23 

1  63 

91 

94 

3  20 

.94 

1.91 

2.06 

SPRING  FALL  (MAR  .  APR  ,  SEPT  .  OCT  I 


PROB 

1964 

1965 

1966 

1967 

1968 

1969 

1970 

1971 

REP 

95 

27  53 

~29A3~ 

32  13 

28  50 

29  15 

26.20 

20.36 

25.45 

26.25 

.90 

20.61 

27  35 

28  05 

23.75 

26.40 

23  03 

18.60 

21.68 

23.33 

80 

17.40 

20  31 

18.60 

18  38 

18.94 

17  94 

15.24 

18.76 

17.78 

70 

12  85 

16  07 

13.60 

13  67 

14.47 

14.60 

11.76 

14.90 

13.92 

.60 

10.70 

12  54 

996 

11  00 

11  48 

11.80 

8.85 

11.43 

10.70 

.50 

8.63 

9.92 

8.60 

9  25 

8.80 

9.00 

6.50 

8.88 

8.21 

.40 

6.65 

8  53 

5.70 

7.25 

7  43 

7.73 

4.84 

7.23 

6.69 

.30 

3.78 

6  36 

4.44 

5  50 

5  80 

6.45 

3.04 

3.92 

4.70 

,20 

2.68 

4.01 

2.64 

3.00 

3  32 

7.72 

2.12 

2  88 

3.07 

.10 

1.83 

1  75 

1.66 

1.65 

1.84 

1.37 

54 

1.51 

1.26 
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TABLE  9  1  P.M.  CGNTRAST  TEMPERATURE  FOR  SHIP  J(1964-71) 


PROS 

1964 

1965 

1966 

1967 

1968 

1969 

1970 

1971 

REP 

.95 

22.29 

18.40 

24.80 

21.68 

24.38 

24.20 

20.48 

25.60 

22.00 

.90 

18.30 

14.27 

19.54 

17.27 

22.00 

19.05 

17.69 

21.20 

18.14 

.80 

10.80 

9  80 

14.60 

11.70 

16.67 

13.45 

14.20 

12.10 

13.24 

.70 

8.29 

8  64 

11.63 

9.43 

12.00 

10.70 

10.60 

10.40 

10.15 

.60 

6.80 

7.53 

9.55 

7.83 

9.25 

8.15 

8.37 

7.48 

8.18 

.50 

5.66 

5.60 

7  63 

681 

6.92 

6.00 

6.33 

5.88 

6.62 

.40 

4.16 

4.36 

5.77 

5  73 

6.10 

4.54 

5.22 

4.28 

4.97 

.30 

2.82 

3.21 

3.44 

3.34 

4.73 

3.22 

4.22 

2.B8 

3.78 

.20 

1.84 

2.22 

2.38 

2  68 

2  00 

2  13 

2.49 

2.16 

2.26 

10 

.57 

.84 

1  34 

1  62 

90 

.72 

.88 

.38 

1.00 

WINTER  (JAN  ,  FEB..  NOV  .  DEC.) 


PROB 

1964 

1965 

1966 

1967 

1968 

1969 

1970 

1971 

REP 

95 

38.81 

36.00 

35.25 

31  94 

33  03 

36  05 

34.00 

36.06 

35.38 

.90 

32.63 

31  17 

29  35 

25.26 

30  50 

30.13 

22.95 

29.78 

27.79 

.80 

23.66 

25.68 

23.40 

19.37 

22.36 

22.20 

17.20 

24.17 

21.44 

70 

18.78 

20.66 

16.80 

15.10 

16.85 

17  33 

12.56 

19.10 

16.61 

.60 

15.80 

17  32 

12.93 

12  71 

11.87 

14.91 

10.68 

14.80 

14.00 

.50 

10.88 

13  80 

11.33 

11.17 

8  92 

13.13 

9.25 

12.25 

11.53 

.40 

8.43 

10  85 

8  80 

7  77 

6.20 

10.15 

8.64 

10.20 

8.53 

.30 

5  93 

7  13 

6  80 

3.68 

3  53 

8.43 

6.70 

7.48 

5.98 

.20 

2  58 

4  64 

4.28 

2.28 

2.12 

4.57 

4.56 

4.12 

3.38 

10 

1  33 

2  46 

1  12 

1  06 

1  27 

2  78 

1.48 

2.22 

1.92 

SPRING  FALL  (MAR.,  APR  .  SEPT..  OCT.) 


§ 

a 


PROB 

1964 

1365 

1966 

1967 

1 96S _ 

1969 

1970 

1971 

REP 

.95 

26.10 

32  00 

32  75 

29.31 

28.05 

28.88 

23.66 

26.60 

28.21 

90 

20.34 

24.45 

28  70 

26  03 

25.16 

25  80 

20  61 

20.76 

24.52 

.80 

17  07 

19.95 

20.10 

21,72 

19.20 

19.50 

16.13 

18.12 

18.04 

70 

13.13 

16.70 

14.64 

17.55 

14  70 

14.50 

12.58 

14.20 

15.07 

.60 

11.31 

13.47 

11.20 

12.27 

11.95 

12.22 

10.20 

11.93 

11.84 

.50 

900 

9.70 

8  90 

10.30 

9.75 

9.80 

7.58 

9.25 

8.94 

.40 

7.45 

7.65 

7.23 

8.62 

7.90 

7.25 

5,58 

7.87 

7.10 

.30 

4.84 

5.20 

5.47 

5.11 

6.24 

3  80 

4.34 

5.35 

5.02 

.20 

2.74 

3.38 

4.02 

3.86 

3.32 

2.72 

3.08 

3.35 

3.38 

,10 

1.97 

1.83 

2.38 

1.86 

2.06 

2.07 

1.38 

1.38 

1.88 
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TABLE  1C  2  P.M.  CONTRAST  TEMPERATURE  FOR  SHIP  J(1964-71> 


PROB 

1964 

1965 

1966 

1967 

1968 

1969 

1970 

1971 

REP 

.95 

21.80 

18.50 

22.73 

22  05 

23.00 

24.30 

19.35 

24.87 

21.69 

.90 

16.60 

13  75 

17.15 

18.66 

18.75 

19.05 

16.35 

10.36 

16.56 

.80 

11.10 

10.50 

12.52 

13  40 

14.50 

13  80 

10.93 

11.27 

12.50 

.70 

7.98 

7.75 

10  32 

9  58 

12.00 

9,55 

9.68 

8.93 

9.88 

.60 

6.44 

6.30 

8  53 

7.81 

8  50 

6.77 

7.70 

7.56 

7.42 

.50 

5.00 

5.31 

6.08 

6.94 

6.75 

5.50 

6.21 

5.75 

5.97 

.40 

3.75 

4.00 

4.04 

4  07 

6.00 

4.39 

5.03 

3.96 

4.88 

.30 

2.36 

3.10 

2.47 

3  97 

4  07 

3.18 

3.65 

2.73 

3.27 

.20 

1.52 

2.30 

2.01 

2.56 

1.80 

1.44 

2.32 

1.95 

2.00 

.10 

.44 

1  60 

1.27 

1  58 

.45 

44 

.66 

.47 

1.02 

PROB 

1964 

WINTER  (JAN.,  FEB 

1965  1966  1967 

.  NOV., 

1968 

DEC  ) 

1969 

1970 

1971 

REP 

.95 

30.38 

36  06 

31.83 

39.00 

28.30 

36  00 

31.63 

34.50 

32.15 

.90 

27.18 

31.42 

28  37 

28.65 

25  73 

27.40 

24.45 

30.00 

27.94 

.80 

19  98 

25  44 

23  13 

21.20 

20  70 

20.90 

18.30 

22.50 

21.87 

.70 

16.45 

21.37 

19  87 

14  87 

16  53 

15.90 

13.98 

18.43 

17.68 

60 

13  55 

15  80 

16  60 

13  31 

12.80 

14.40 

10.80 

15.40 

13.70 

.50 

11  08 

12  75 

13.00 

10  80 

7  83 

12.00 

9.54 

11.67 

10.42 

.40 

7.73 

9  70 

10  05 

8.20 

6  10 

9  73 

7.85 

7.26 

8.98 

30 

4.02 

7  48 

6.35 

5  23 

3  48 

7.35 

6.76 

5.88 

6  48 

.20 

254 

4  44 

3  32 

264 

1  % 

4  68 

4.12 

3.40 

3.32 

10 

1  19 

1  87 

1  38 

1  86 

63 

1  96 

2.27 

1.80 

1.46 

PROB 

1964 

SPRING  FALL  (MAR.,  APR  .  SEPT  .  OCT.) 

1965  1966  1967  1968  1969 

1970 

1971 

REP 

95 

24  45 

29  00 

31.69 

27  30 

28.61 

27  25 

27.20 

25.73 

28.07 

.90 

21  80 

24  75 

26  10 

23  85 

25.38 

22  50 

23.10 

21.60 

23.85 

.80 

17  47 

19  50 

19.80 

19  92 

19  54 

16.80 

17.45 

15.64 

17  78 

.70 

13.70 

16.00 

14.53 

15  80 

16.26 

14.13 

13.60 

13.47 

14.74 

.60 

10.47 

1343 

10.96 

11  34 

12.30 

11  50 

10.60 

11.10 

11.95 

.50 

7  17 

9.25 

8.38 

9.88 

10.36 

9.35 

9.13 

9.50 

8.77 

.40 

4.81 

8.00 

5.10 

7.30 

8.47 

7.63 

6.60 

7.93 

6.64 

30 

3.47 

5  92 

4  35 

4  44 

7  15 

4.53 

4.38 

6.10 

5.31 

.20 

2.31 

4.73 

3.08 

3  16 

5.57 

2.50 

3.04 

3.54 

3.94 

10 

1.82 

1.80 

2.33 

1.87 

2.54 

1  45 

2.12 

1.63 

2.00 
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TABLE  11  3  P.M,  CONTRAST  TEMPERATURE  FOR  SHIP  JU964-71) 


PROB  1964 


.96 

17.40 

.90 

13.40 

.80 

9.20 

.70 

6.48 

.60 

5.63 

.50 

5.00 

.40 

3.83 

.30 

1.97 

.20 

1.33 

.10 

.36 

PROB 

1964 

.96 

29.40 

.90 

23.10 

.80 

17.68 

.70 

14.20 

.60 

12.28 

.50 

9.50 

.40 

6.70 

.30 

4.68 

.20 

2.56 

.10 

.96 

1965  1966 


13.66 

25.80 

12.15 

19.20 

8.45 

11.90 

7.15 

9.03 

5.20 

5.93 

4.15 

5.00 

3.05 

3.59 

2.17 

2.76 

1.62 

2.02 

.85 

1.16 

1965 

1966 

29.25 

27.79 

24.75 

25.05 

21.30 

19.92 

17.63 

15.43 

14.25 

10.85 

11.50 

8.25 

9.00 

6.93 

6.18 

4.83 

3.00 

2.16 

1.48 

.89 

1967 

1968 

17.50 

18.85 

14.75 

16.28 

10.00 

10.73 

7.50 

8.05 

6.00 

6.23 

4.87 

5.04 

3.60 

3.96 

2.90 

2.38 

2.17 

1.46 

1.47 

.32 

1967 

1968 

32.75 

27.93 

27.30 

24.10 

17.30 

20.07 

13.73 

16.20 

11.20 

13.70 

9.58 

11.17 

7.73 

6.15 

5.68 

4.14 

2.68 

2.08 

1.74 

1.59 

1969  1970 


18.75 

21.23 

13.75 

14.65 

10.00 

11.15 

7.95 

9.19 

6.50 

7.05 

4.40 

4.80 

3.00 

4.24 

2.03 

3.01 

.87 

1.54 

.10 

.27 

1969 

1970 

29.70 

31.67 

23.10 

22.71 

17.20 

18.43 

13.80 

13.71 

11.60 

12.00 

9.00 

9.00 

7.45 

7.33 

5.65 

5.83 

3.28 

3.20 

1.66 

1.20 

1971 

REP 

17.26 

19.73 

14.00 

15.68 

9.42 

10.18 

8.60 

7.84 

6.33 

6.13 

5.64 

4.96 

4.20 

3.62 

2.70 

2.49 

1.70 

1.52 

.50 

.79 

1971 

REP 

29.00  30.27 

26.60  26.20 

20.60  19.26 

18.00  16.67 

12.26  12.66 

10.17  9.88 

6.88  7.68 

5.00  5.16 

3.60  2.79 

1.60  122 


WINTER  (JAN..  FEB..  NOV..  DEC.) 


SPRING  FALL  'MAR..  APR..  SEPT..  OCT.) 


PROB 

1964 

.95 

28.95 

.90 

23.16 

.80 

18.10 

.70 

12.95 

.60 

9.84 

.50 

7.50 

.40 

6.88 

.30 

5.55 

.20 

3.06 

10 

1.89 

1965 

1966 

26.65 

28.65 

24.30 

24.15 

18.12 

18.20 

14.90 

13,30 

11.08 

10.69 

8.67 

8.83 

6.93 

6.60 

5.37 

4.88 

2.98 

3.08 

1.67 

1.19 

1967  1968 


27.25 

26.43 

21.50 

23.33 

17.00 

18  09 

13.08 

13.58 

11.33 

11.20 

9.15 

9.25 

7.50 

7.15 

5.25 

6.31 

2.80 

4.56 

1.57 

2.58 

1969  1970 


?7.08 

23.70 

21.83 

18.90 

16.05 

14.43 

12.72 

12.49 

11.07 

8.U. 

8.79 

7 .7.0 

6.26 

4.9* 

3.97 

3.56 

zm 

1.5- 

1.50 

1971 

REP 

25.50 

26.33 

21.00 

21.60 

13.4E 

16.32 

t-33 

13.7C 

11,00 

d.94 

r.  "a 

a.23 

7.25 

'v.tJ 

6 

4.91 

3.0  r 

3.11 

1.90 

SUMMER  (MAY.  JUN?:.  JULY.  AUG.f 
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TABLE  12  NIGHT  CONTRAST  TEMPERATURE  FOR  SHIP  J(1964-71) 


PROB 

1964 

1965 

1966 

1967 

1968 

1969 

1970 

1971 

REP 

95 

“086” 

10.47 

8.90 

1C  15 

10,18 

10.13 

10.39 

10.06 

9.88 

.90 

8.83 

8.61 

7.31 

8.12 

8.70 

7.39 

9.18 

7.25 

8.22 

.80 

6,35 

6  19 

5.31 

5.60 

6.27 

5.47 

6.28 

4.96 

5.6S 

70 

4.59 

4.39 

3.93 

3.98 

4.19 

3.90 

4.33 

3.70 

4.15 

.60 

,.32 

2.97 

3.07 

3  14 

3.01 

2.90 

3.00 

2.92 

3.11 

.50 

2.51 

2  25 

2.35 

260 

2.37 

2.21 

2.21 

2.37 

2.41 

40 

206 

1.72 

1  80 

2.14 

1.84 

1.62 

1.60 

1.86 

1.87 

.30 

:  58 

1.29 

1.23 

1.71 

1.34 

1.06 

1.02 

1.41 

1.30 

.20 

VOS 

75 

.54 

1.21 

79 

.45 

.44 

.93 

.83 

.10 

.24 

-  .15 

-  28 

.47 

-  .30 

-  .62 

-  .40 

.02 

.19 

PROB 

1964 

WINTER  (JAN.,  FEB.,  NOV.,  DEC.) 

1965  1966  1967  1968  1969 

1970 

1971 

REP 

95 

10.04 

8.03 

9.88 

10.06 

7.97 

10.64 

9.80 

7.96 

9.30 

90 

8.09 

7.09 

7  46 

7.89 

6.88 

8.46 

7.75 

6.80 

7.63 

80 

6  40 

5.54 

5  86 

5  49 

5.30 

6.08 

6.07 

5.19 

5.86 

70 

A  0? 

A  A  ^ 

A  on 

4  06 

4.08 

4.43 

4.51 

3.90 

4.39 

.60 

3.62 

3.48 

3.75 

3.20 

3.16 

3.26 

3.43 

3.19 

3.38 

50 

2  77 

2  79 

3.02 

2.58 

2.40 

2.55 

2.74 

2.61 

2.71 

40 

2  29 

2.27 

2  47 

2  13 

1.95 

2.05 

2.16 

2.19 

2.21 

.30 

1  84 

1.77 

1.90 

1.73 

1.49 

1.53 

1.59 

1.68 

1.67 

.20 

1  33 

1  20 

1  29 

1  23 

88 

1.02 

.98 

1.26 

1.13 

10 

66 

39 

60 

.18 

.02 

36 

.17 

.55 

.34 

PROB 

1964 

SPRING  FALL  (MAR  . 

1965  1966  1967 

APR., SEPT.,  OCT.) 

1968  1969  1970 

1971 

REP_ 

95 

7.82 

7  31 

6  80 

6  39 

7.32 

7.19 

6.84 

6.07 

6.96 

90 

6  51 

6  21 

5  97 

5.66 

6.41 

5.85 

5  23 

4.99 

6.75 

.80 

4  71 

4.87 

4.55 

4.33 

5.10 

4.38 

3.63 

3.89 

4.37 

70 

3.55 

3  73 

3.48 

3.24 

4  06 

3.43 

2.86 

3.06 

3.46 

.60 

2.83 

2.95 

2  77 

2.64 

3.11 

2.78 

2.24 

2.46 

2.68 

.50 

•>  25 

2  51 

2  33 

2.21 

2.47 

2.31 

1.89 

2.12 

2.20 

40 

1  87 

2  12 

1  95 

1  87 

2.04 

2.00 

1.64 

1.77 

1.88 

.30 

1  54 

1.76 

1.60 

1.56 

1  76 

1.63 

1.31 

1.40 

1.64 

.20 

1  19 

1.31 

1.25 

1.20 

1  37 

1.31 

.91 

1.06 

1.14 

10 

.70 

.78 

62 

63 

.65 

.74 

.31 

.60 

.65 
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TABLE  13  TOTAL  CONTRAST  TEMPERATURE  FOR  SHIP  JU964-71) 


PROB 

1964 

1965 

1966 

1967 

1968 

1969 

1970 

1971 

REP 

.95 

12.09 

11.30 

12.41 

12.35 

13.11 

12.09 

12.87 

12.33 

12.21 

.90 

9.76 

9.28 

9.43 

9.67 

10.46 

9.65 

10.49 

9.41 

9.89 

.80 

7.15 

6.65 

6.76 

6.95 

7.49 

6.40 

7.88 

6.38 

6.94 

.70 

5.29 

5  14 

5.01 

5.14 

5.66 

4.75 

5.53 

4.61 

5.14 

.60 

3.88 

3.64 

3.69 

3.78 

3.93 

3.49 

4.03 

3.42 

3.76 

.50 

2.83 

2.67 

2.77 

3.00 

2.3'. 

2.68 

2.84 

2.72 

2.84 

.40 

2.26 

2.04 

2.12 

2.44 

2.16 

1.92 

1.99 

2.14 

2.18 

.30 

1.71 

1.48 

1.52 

1.95 

1.55 

1.27 

1.38 

1.60 

1.61 

.20 

1.15 

.91 

.82 

1.35 

.92 

.60 

.63 

1.02 

.98 

.10 

.38 

02 

-.16 

.62 

-  .00 

-  .39 

-  .24 

.10 

.00 

WINTER  MAN..  FEB..  NOV.,  DEC.I 


PROB 

1964 

1965 

1966 

1967 

1968 

1969 

1970 

1971 

REP 

.95 

20.38 

22.22 

20.42 

17.73 

19.84 

19.63 

16.88 

20.33 

19.56 

.90 

13.84 

15.72 

13.81 

13.00 

13.26 

13.94 

11.93 

13.91 

13.83 

.80 

9.22 

9.2S 

9.02 

8.61 

7.94 

9.71 

8.47 

8.48 

8.61 

.70 

6.82 

6.67 

6.60 

6.16 

5.90 

7.20 

6.45 

6.24 

6.56 

60 

5.18 

5.08 

5.17 

4.45 

4.49 

5.34 

4.86 

4.51 

4.90 

.50 

3.68 

3.82 

3.99 

3.32 

3.38 

3.86 

3.62 

3.51 

3.66 

.40 

2.72 

2.88 

3.02 

2.5? 

2.52 

2.79 

2.76 

2.75 

2.77 

.30 

2.12 

2  15 

2.26 

2.01 

1.90 

2.07 

1.99 

2.18 

2.06 

.20 

1.53 

1.50 

1.52 

1.51 

1.26 

1.42 

1.24 

1.61 

1.39 

.10 

.79 

.74 

.63 

.48 

.33 

.70 

.35 

.78 

.56 

SPRING-FALL  (MAR.,  APR..  SEPT  .  OCT.) 


PROB 

1964 

_1965 

1966 

1967 

1968 

1969 

1970 

1971 

REP 

.95 

16  81 

18.55 

18.89 

17.53 

18.19 

16.85 

1496 

15.83 

16.93 

.90 

IV-*- 

13.23 

12.47 

11.94 

12.47 

12.30 

10.56 

11.23 

11.90 

.80 

7  ij 

8.10 

7.37 

7.4F 

8.06 

7.84 

6.56 

7.15 

7.33 

.70 

5.67 

6.07 

f.43 

5.5* 

6  18 

5.7.: 

4.58 

5.18 

5.38 

.60 

4.18 

4.72 

1.10 

4.12 

4.75 

4.23 

3.40 

3.88 

4.08 

.50 

3.27 

3.59 

3.19 

3.08 

3.63 

3.2b 

2.72 

3.06 

3.18 

.40 

7.58 

2.78 

2.50 

2.45 

2.75 

2.53 

2.12 

2.39 

2.45 

.30 

2.03 

2.22 

1.96 

1.95 

2.12 

2.07 

1.71 

1.89 

1.97 

.20 

1.56 

1.73 

1.51 

1.51 

1.69 

1.56 

1.24 

1.37 

1.49 

.10 

1.00 

1.06 

.95 

.92 

.98 

1.02 

.53 

.72 

.80 

SUMMER  (MAY,  JUNE.  JULY,  AUG.) 


TEMPERATURE  ADDITION  OR  LESS 


NSWC/WOL  TR  78-187 


over  the  wavelength  region  of  interest  (8-12,-m).  The  ship  and  background  tempera¬ 
tures  were  designed  to  approximately  reproduce  the  target  surface  radiance  and  the 
background  radiance  respectively  when  used  in  the  blackbody  spectral  radiance 
(Planck's)  equation  Equation  (3)1.  Tiius,  these  temperatures  are  used  to  compute 
spectral  radiances  and  these  are  integrated  numerically  over  the  8-12pm  region  to 
find  the  ship  contrast  radiance.  Contrast  temperatures  are  composed  of  Independent 
ship  and  background  temperature  components.  While  the  radiance  from  each  of  these 
two  sources  is  uniquely  related  to  the  real  or  effective  temperature,  there  is  not 
a  unique  radiance  associated  with  their  temperature  difference.  This  is  illustra¬ 
ted  in  Figure  17  which  shows  a  range  of  radiant  contrasts  corresponding  to  back¬ 
ground  temperatures  of  5C  and  1 5 ' C . 


Figure  17.  Relationship  between  contrast  temperature  and 
radiant  contrast  ship  signatures 


These  curves  show  that  the  range  of  radiance  contrasts  corresponding  to  a  con¬ 
trast  temperature  increases  with  contrast  temperature.  However,  this  range  is  not 
large  for  in  the  background  and  contrast  temperatures  encountered  1r,  the  present 
work,  i.e.,  median  contrast  temperatures  of  5’C  or  less.  In  fact,  an  average  value 
chosen  from  Figure  17  will  adequately  relate  any  of  the  previously  given  contrast 
temperatures  to  the  target  radiant  contrast.  Recommended  average  values  will  be 
given  in  the  next  section  of  this  report. 
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CONCLUSION 


This  report  has  presented  the  results  of  a  study  whose  goal  was  to  provide 
ship  signatures  in  a  form  useful  to  SCAFIRf  system  perfonnance  studies.  Since  the 
SEAFIRE  thermal  imaging  sensor  operates  in  the  8-12, m  region,  all  work  of  the 
present  study  was  directed  towards  that  waveband.  The  results  presented  and  the 
conclusions  drawn  will  not  hold  for  the  3-5. ,m  or  any  other  waveband.  However,  the 
methods  developed  could  be  extended  to  the  3-5.  m  waveband  by  considering  the  effect 
of  the  ship  exhaust  plume  and  redefining  the  background  temperatures.  SEAFIRE 
system  studies  required  a  statistical  distribution  of  target/background  temperatures 
which  include  weather  effects  typically  encountered  at  sea.  It  was  found  that 
existing  methods  of  determining  ship  signatures  are  long  and  involved  computation¬ 
ally.  Hence,  a  relatively  large  amount  of  expensive  computer  time  would  be  required 
to  present  an  analysis  of  the  effect  of  weather  on  ship  signatures.  Therefore,  a 
simpler  method  of  computing  contrast  temperatures  was  devised  so  that  a  large 
number  of  signatures  using  weather  data  could  be  generated.  It  was  shown  that  this 
method,  named  the  Single  Element  Method,  provides  a  good  approximation  to  the 
average  ship  temperature.  This  method  will  accurately  include  the  effect  of  a  hot 
stack  on  average  ship  temperature.  However,  the  average  ship  temperature  is  not 
meaningful  in  cases  where  the  small  stack  surface  area  contributes  half  or  more  of 
the  total  temperature.  In  these  cases,  the  stack  will  probably  be  sensed  at  a  much 
greater  distance  than  the  remaining  ship.  Furthermore,  the  most  obvious  counter¬ 
measure  for  a  ship  signature  is  to  cool  the  stack  hence  it  nay  be  assumed  that 
future  warships  will  have  some  form  of  cooled  stack.  The  average  ship  temperatures 
and  resulting  contrast  temperatures  calculated  were  for  a  ship  model  having  a 
cooled  stack.  For  this  case,  it  was  shown  that  differences  in  ship  construction  or 
internal  temperature  distribution  do  not  signi f icantly  effect  its  average  tempera¬ 
ture  hence  the  results  given  will  have  general  appl ication.  Improvements  in  the 
Single  Element  Method  could  still  be  made.  The  most  noteworthy  improvement  being 
the  obtaining  of  a  model  atmosphere  that  varies  continually  with  day  of  the  year. 
Also,  any  advance  in  methods  of  computing  ship  signatures  could  be  incorporated  into 
this  method  as  they  become  available. 

It  is  difficult  to  summarise  the  ship  contrast  temperatures  given  in  this 
report  because  they  are  greatly  influenced  by  several  variables.  First,  solar 
heating  influences  ship  temperatures  hence  the  target  ships  heading  and  its 
location  on  the  earth  are  important.  Secondly,  weather  conditions  (i.e.,  rain, 
wind,  air  temperature  and  cloud  cover)  largely  govern  daily  signatures  and  these 
can  also  depend  on  location  or  season  of  the  year.  Finally,  the  relative 
positioning  between  the  target  and  the  sensor  could  be  important,  especially  in 
determining  backgrounds.  A  summary  of  day,  night  and  total  ship  contrast  tempera¬ 
tures  provided  in  the  previous  section  is  given  in  Table  M  for  a  sea  surface 
background.  This  summary  consists  of  the  probability  of  observing  contrast 
temperature  ranges  or  lower  values  for  year-round  conditions  when  the  effect  of 
season  and  ship's  heading  are  considered.  Here,  the  heading  of  27u  (due  west)  on 
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a  starboard  element  is  excluded  from  daytime  contrast  temperature  because  it  leads 
to  contrast  temperatures  very  similar  to  nighttime  values. 


Table  14. 
Probabll  ity 
0.90 


Range  of  Computed  Ship  Contrast  Temperatures 
Paytime  Range  Nighttime  Range  Total  Time  Range 
8-20'C  6-8°C  8-1 3°C 

3- 5 0  C  2-2. 5°C  2.5-3.5°C 

0.5-rc  0-0. 5!C  0-1  °c 


Table  14  can  be  used  to  provide  quic 
peratures  which  can  also  be  converted  to 
However,  the  limitations  imposed  by  using 
remembered.  First,  Ship  J  was  at  a  midla 
tures  are  expected  to  vary  with  latitude, 
signatures  were  computed  that  the  locatio 
amount  of  cloudy  weather.  Cloud  conditio 
increase  background  temperatures  thereby 
tures.  Hence  the  contrast  temperatures  1 
mates  for  the  midlatitude  location.  Fina 
ship  target  f’-om  another  ship  and  do  not 
or  elevated  ground  station. 


k  estimates  of  ship  target  contrast  tem- 
contrast  radiance  by  use  of  Figure  16. 
the  Ship  J  weather  data  should  be 
titude  location  ( 5 1 0 N )  and  these  signa- 
Second,  it  was  noticed  after  these 
n  of  Ship  J  contains  an  extraordinary 
ns  act  both  to  lower  ship  temperatures  and 
leading  to  low  values  of  contrast  tempera- 
i s tod  in  Table  14  are  conservative  esti- 
11 y ,  these  estimates  are  for  viewing  a 
hold  for  viewing  the  ship  from  an  aircraft 
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SYMBOLS 

a  solar  absorptance 

A  area  fraction  of  ship  surface  are  with  respect  to  internal  temperature 

AF  angle  factor  for  solar  reflection  off  of  a  ship 

AS  ship  surface  area,  m 

A z  solar  azimuth  angle,  degrees  from  north 

b  solar  hour  angle,  degrees 

B  area  fraction  of  ship  surface  area  with  respect  to  heat  capacity 

c  latitude,  degrees 

C  area  fraction  of  ship  surface  are  with  respect  to  ship  stack 

d  solar  declension,  degrees 

D  grouping  of  variables  for  internal  temperature  effects  !  see  equation  (25)) 

L  second  grouping  of  variables  for  internal  temperature  effects 

I  see  equation  (26) 

f  fraction  of  dear  sky 

F  grouping  of  variables  for  thermal  capacity  effects  [see  equation  (21)] 

G  second  grouping  of  variables  for  thermal  capacity  effects  [see  equation 

(22) 

ha  convective  coefficient  for  air  flowing  over  a  vertical  surface, 
watts/m'  -  C 

hi  internal  heat  transfer  coefficient,  watts/m;  -  9C 

hk  convective  coefficient  for  stack  gases,  watts/m;  -  °C 

hr  convective  coefficient  for  rain  on  a  vertical  surface,  watts/m5  -5C 

H  altitude  above  sea  level 

HA  dimensionless  parameter  for  air  cooling 

HI  dimensionless  parameter  for  interna)  heating 

HK  dimensionless  parameter  for  stack  gas  heating 

HR  dimensionless  parameter  for  rain  cooling 

•  length  of  air  path  from  target  to  sky  background 

M  target  ship  heading,  degrees  clockwise  from  north 


57 


NSWC/WOl  TP  73-187 

sea  surface  or  background  radiance,  watts/nr  -  sr  -  ;<m 

cloud  radiance,  watts/nr’  -  sr  -  ..m 

apparent  background  radiance,  watts/nr  -  sr  -  ;im 

sky  radiance,  watts/nr  -  sr  -  ..m 

target  surface  radiance,  watts/nr  -  sr  -  :,m 

radiative  signal  power,  watts/m* 

direct  solar  energy  incident  on  a  vertical  element,  watts/mJ 
diffuse  solar  energy,  watts/rr 

indirect  or  diffuse  solar  energy  on  a  vertical  element,  watts/nr 
direct  solar  energy,  watts/nr 

total  solar  energy  absorbed  on  a  vertical  element,  watts/nr 

spectral  hlackbody  ••adiance,  watts/nr  -  sr  -  ,.m 

spectral  radiance  of  sun,  watts/nr  -  sr  -  ,.ir. 

atmospheric  radiance  emission,  watts/m"  -  sr  -  ,.m 

atmospheric  emission  of  a  clear  sky.  watts/nr  -  sr  -  ..m 

atmospheric  emission  of  path  from  target  to  clouds,  watts/m?  -  sr  -  pm 

tempera  turc 

basic  element  or  average  ship  temperature,  cC 

ship  target  contrast  temperature  ,  sec-  equation  (39)  i 

internal  temperature  correction  factor,  :C 

thermal  capacity  temperature  correction  factor,  :C 

stack  presence  temperature  correction  factor,  C 

ship  target  contrast  temperature  ;  see  equation  (39)',  5C 

temperature  of  e<haust  gases  inside  ship  stack,  'C 

internal  temperature  of  ship  compartments,  C 

sea  level  temperature  or  initial  ship  temperature,  'C 

effective  average  ship  temperature  'see  equation  (30)! 

sea  surface  water  temperature.  C 

effective  sky  temperature,  C 

atmospheric  transmission 

loss  in  direct  solar  energy  due  to  aerosol  scattering 
loss  in  direct  solar  energy  due  to  molecular  scattering 
solar  energy  transmission 

relative  wind  speed  between  air  and  ship  velocities,  m/sec 
ratio  of  time  step  to  thermal  capacity  ! see  equation  (30)] 
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thermal  capacity  of  ship  surface  sections,  watt-hrs/m:'  -  ”C 
solar  zenith  angle 

target-detector  zenith  angle  ! see  Figure  (1)1 
solid  angle  of  target  area  in  detector  f iel d-of-view,  sr 
normalized  detector  response 
emissivity  in  8-12,.m  waveband 

time  step  in  numerically  computing  ship  temperatures,  hours 
wavelength,  ..m 

elevation  angle  of  sea  surface  facets,  degrees 
reflection  angle  of  sea  surface  facets,  degrees 


SUBSCRIPTS 


Clouds 

internal  temperature  divisions  of  ship 


thermal  capacity  divisions  of  ship 
stacks  of  ship 

ship 


sea  suridce 
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APPENDIX  A 


APPLICATION  OF  SINGLE  ELEMENT  METHOD 


This  appendix  will  describe  how  the  Single  Element  Method  is  applied  to  a 
given  ship,  the  Patrol  Frigate,  to  evaluate  the  correction  factors  referred  to  in 
equations  (27),  (33)  and  (36).  The  basic  element  of  this  method  is  simply  a 
vertical  rectangle  of  unit  area  which  has  the  lowest  internal  temperature  and 
lowest  thermal  capacity  of  all  ship  sections  considered.  This  element  is  basically 
the  same  for  all  similarly-constructed  ships.  The  magnitude  of  the  correction 
factors  vary  with  differences  in  ship  construction  and  operation  and  cause  dif¬ 
ferences  in  average  temperatures  between  ships.  However,  it  will  be  shown  that 
these  correction  factors  are  small  in  comparison  with  the  basic  element  temperatures 
as  computed  by  the  Single  Element  Method  for  the  Patrol  Frigate  are  a  good  approxi¬ 
mation  for  all  similar  ships.  Here,  similar  ships  are  those  having  aluminum  super¬ 
structures  and  steel  hulls  of  equivalent  thickness  and  relatively  small  cooled 
stacks  as  those  found  on  the  Patrol  Frigate  Model  . 

A  portion  of  the  Patrol  Frigate  Model  indicating  the  Internal  temperatures, 
neat  transfer  coefficients  and  thermal  capacities  is  given  in  Table  1.  Here,  a 
broadside  view  of  the  starboard  side  of  the  snip  is  divided  into  33  sections  and 
these  variables  are  specified  along  with  the  surface  area  for  each  section.  Table 
A-l  shows  that  the  internal  heat  transfer  coefficient  is  0.625  watts/m;'  -  °C  over 
most  of  the  ship  hence  this  value  is  used  in  the  single  element  method.  The 
starboard  view  is  arbitrarily  chosen  as  any  view  can  be  obtained  from  the  3- 
dimensional  Patrol  Frigate  Model.  In  specifying  the  starboard  view  in  Table  1, 
the  3-dimens ional ity  of  the  ship  is  neglected.  That  is.  all  of  the  sections  listed 
are  assumed  to  be  vertical  and  the  slight  tilt  of  some  sections  is  not  considered. 
Correction  factors  are  derived  from  Table  1  by  considering  the  ship  first  in  terms 
of  its  internal  temperatures  and  subsequently  in  terms  of  its  thermal  capacity. 

The  internal  temperatures  are  modeled  by  first  grouping  the  sections  listed 
into  3  temperature  regions.  These  are:  (1  )  80-95r,F;  (2)  1 00- 1 1 0° F ;  and  (3)  120- 
1 40° F .  The  area  fractions  for  each  of  these  regions  are  found  by  summing  their 
surface  areas  and  dividing  by  the  total  surface  area  of  the  starboard  view  (1130  m  ). 
The  following  area  fractions  result  and  are  applied  at  the  listed  intermediate  tem¬ 
peratures  . 


Temperature 
Area  Fraction 


0.441 


1051’  F 
0.453 


130°  F 
0.016 


The  temperature  correction  factor  for  internal  temperature  differences  can  be 
derived  from  the  above  table  and  equation  (27).  Let,  A:  =  .453,  A;  =  .106,  Tn? 
Tn.  =  20'F  (U.rc)  and  Tn,  -  Tn;  =  45 1 F  ( 2  5 0  C ) .  Equation  (27)  reduces  to 
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SECTION  IDENTIFICATION 

AS 

Ti 

hi 

W 

Peak  Tank 

35.1 

85 

2.660 

10.4 

Boson  Storeroom 

46.4 

105 

0.625 

10.4 

Hull  Storeroom 

93.5 

105 

0.625 

10.4 

Passage 

21.0 

105 

0.625 

8.65 

Passage 

32.4 

105 

0.625 

10.4 

Crew  Living  Complex 

68.3 

80 

0.625 

10.4 

Passage 

33.4 

95 

0.625 

13.9 

Laundry  Room 

31.4 

80 

0.625 

8.65 

Passage 

26.8 

95 

0.625 

17.3 

Dry  Provisions  Storeroom 

15.6 

85 

0.625 

8.65 

Service  of 

4.8 

105 

0.625 

8.65 

Passage 

31.6 

no 

0.625 

17.3 

Auxiliary  Machinery  Room  2 

21.1 

120 

0.625 

8.65 

Passage 

35.5 

130 

0.625 

13.9 

Engine  Room 

28.2 

140 

0.625 

10.4 

EOS  Room 

30.2 

80 

0.625 

10.4 

Auxiliary  Hachinery  Room  3 

21.6 

120 

0.625 

8.65 

Central  Office  Complex 

33.6 

80 

0.625 

8.65 

Supply  Storeroom 

28.2 

105 

0.625 

8.65 

Gas  Cylinders  ETC.  Room 

35.0 

105 

0.625 

8.65 

Steering  Gear  Room 

30.8 

85 

0.625 

8.65 

Auxiliary  Machinery  Room  1 

13.6 

120 

0.625 

8.65 

Helocopter  Hanger 

176.0 

105 

0.625 

4.09 

Ammunition  Magazine 

139.2 

80 

0.625 

4.09 

WR  Mess  Room 

14.1 

80 

0.625 

4.09 

WR  Pantry 

14.4 

90 

0.625 

4.09 

PI  enum 

23.8 

105 

0.625 

4.09 

Passage 

5.0 

90 

0.625 

4.09 

Helocopter  Shop 

8.4 

100 

0.625 

4.09 

Ammunition  Magazine 

14.2 

85 

0.625 

4.09 

Fan  Room 

10.8 

105 

0.625 

4.09 

OTOM 

5.8 

95 

0.625 

4.09 

Main  Stock 

12.6 

825 

26.0 

11.0 

AS  -  Surface  Area  in  M'  hi  =  Internal  Heat  Transfer  Coefficent 

in  watts/M?-°C 

Ti  =  Internal  Temperature  in  F  W  =  Thermal  Capacity  in  Watt-Hrs/M?-°C 


TABLE  A-l  Internal  Temperature  and  Thermal  Capacity  Model  for 
Starboard  Side  of  Patrol  Frigate 


4  A- 2 
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(A- 1 ) 


NSWC/HOL  TR  78-187 
T  -  Tb  +  6.16  •  E 

Equation  ( A- 1 )  shows  that  the  internal  temperature  corrections  depend  on  E  which  is 
composed  cf  the  thermal  capacity  of  the  basic  elements  W,  the  weather  coefficients, 
ha,  hr,  and  the  time  step,  A*.  i  see  equations  (18)-{20)1.  A  maximum  value  of  this 
correction  can  be  calculated  for  the  Patrol  Frigate  Model.  The  maximum  E  results 
for  minimum  convective  cool  i ng  coefficients .  Here,  hr  =  0  (no  rain)  and  ha  =  6.0 
watts/m2  -  v‘C  (free  convection,  no  relative  wind).  The  thermal  capacity  of  the 
basic  element  is  4.09  watt- hr/nr  -  °C  and  the  time  step  is  1.0  hour  (time  step  has 
only  a  small  influence  on  this  calculation).  With  these  values,  the  following 
value  for  maximum  temperature  correction  factor  results. 

TCI  =  0. 36°C  (A-2) 

It  was  shown  in  Table  3  that  ship  temperatures  up  to  20°C  above  ambient  were 
calculated  from  the  weather  data.  Since  the  higher  ship  temperatures  correspond 
to  the  minimum  cooling  rates  which  produced  the  maximum  value  of  temperature  cor¬ 
rection,  this  factor  will  be  small  compared  to  the  difference  between  ship  and 
ambient  temperature  during  daylight  hours.  At  night,  the  internal  temperatures  ma.y 
contribute  significantly  to  ship  temperature  but  will  not  dominate  the  resulting 
contrast  temperatures.  For  example,  median  nighttime  contrast  temperatures  for 
the  Weather  Ship  J  data  were  found  to  be  about  2.5°C  (see  Table  14).  Thus,  the 
contribution  of  internal  temperature  differences  to  contrast  is  likely  to  be  small 
compared  to  the  solar  heating  effect  in  daytime  or  the  ship  to  background  tempera¬ 
ture  differences  at  night.  Also,  the  rather  large  internal  temperatures  given  for 
some  sections  of  the  F-atrol  Frigate  Model  will  not  have  a  significant  influence  on 
the  results. 

The  Thermal  Capacity  Model  was  derived  by  first  grouping  the  ship  sections 
listed  in  Table  1  into  common  thermal  capacities.  It  was  decided  that  only  the 
three  lowest  thermal  capacities  would  be  considered  and  that  those  thermal  capa¬ 
cities  above  10.4  watt-hr/m'  -  °C  would  be  neglected  because  their  associated  sur¬ 
face  area  is  small.  The  following  area  fractions  were  found  for  the  three  lower 
capaci ties . 


Thermal  Capacity  (watt-hr/m  -  C)  4.09  8.65  10.4 

Area  Fraction  0.408  0.257  0.335 

Here  the  basic  element  thermal  capacity,  4.09  watt-hr/m-’  -  °C,  corresponds  to  the 
vinch  thick  aluminum  plating  found  on  the  superstructure  while  the  higher  values 
correspond  to  the  steel  hulls.  The  temperature  correction  factor  for  thermal 
capacity  differences  can  be  derived  from  the  above  table  and  equation  (33).  Let, 

B.  =  .256,  B,  =  .335  and  V;  =  .244  ,  V.-  =  .116  and  V->  =  .096  for  a  time  step.  Ax,  of 
one  hour.  Equation  (33)  can  be  written  as 

t  .  ib-  -08z5  ^Ve-rrsk5  (A-3) 

where  the  slight  difference  in  denominator  between  the  two  terms  of  the  temperature 
correction  factor  has  been  neglected.  The  temperature  correction  factor  in 
Equation  (A-3)  is  a  complex  combination  of  the  solar  heating,  internal  heating  and 
convective  cooling  rates  which  define  the  factors  F  and  G  !  see  equations  (31)  and 
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(32)  i.  An  estimate  of  this  correction  is  made  by  assuming  that  the  initial  ship 
temperature  is  20°C  for  the  minimum  cooling  conditions  previously  defined.  The 
following  thermal  capacity  correction  is  found  for  an  average  absorbed  solar  heating 
rate  of  75  watts/m:  and  an  average  air  temperature  of  10°C, 

TC2  =  -0. 42°C  (A-4) 

Thus,  the  contribution  of  thermal  capacity  differences  to  the  average  ship  tempera¬ 
ture  calculated  by  the  Single  Element  Method  is  also  likely  to  be  small  compared 
to  the  other  contributions. 

The  stack  correction  factor  for  the  single  stack  found  on  the  Patrol  Frigate 
is  estimated  using  Equation  (37)  and  noting  that  the  stack  area  fraction  is  .0106 
from  the  data  in  Table  1.  The  following  estimate  of  maximum  stack  correction 
results  by  using  the  minimum  cooling  rates. 

TC3  =  2 . 38°C  (A-5) 

If  the  previously  reccnmended  average  cooling  condition,  HC + HR  =  5.0,  is  used,  the 
stack  correction  factor  drops  to  1.50  C.  Table  14  indicates  that  median  contrast 
temperatures  for  the  Weather  Ship  J  data  are  about  2.5°C  at  night  and  5°C  during 
the  day.  Thus,  the  stack  correction  factor  is  likely  to  be  a  significant  portion 
of  the  overall  ship  temperatures  for  hot  stacks.  However,  the  stack  correction 
factor  is  much  reduced  for  ships  having  cooled  stack  or  stack  gases.  For  example, 
if  it  Is  assumed  that  the  stack  gases  are  cooled  to  1 50n F  (66°C),  the  corresponding 
maximum  stack  correction  is  ,363'C  and  the  average  correction  is  only  ,189°C.  A 
stack  gas  temperature  of  150?F  is  artificially  low  to  account  for  stack  casing 
cooling  which  cannot  be  considered  in  the  present  models. 

The  work  of  this  appendix  shows  that  the  temperature  correction  factors  for 
ships  with  cooled  stacks  are  relatively  ‘■mail  compared  to  the  expected  overall  ship 
temperatures.  Hence,  ship  temperatures  calculated  by  the  Single  Element  Method 
using  the  Patrol  Frigate  Model  are  valid  for  all  similarly  constructed  ships.  The 
work  on  stack  correction  factors  shows  that  hot  stacks  contribute  significantly  to 
the  average  ship  temperature  hence  their  effect  should  be  included  in  Single 
Element  Method  calculations  on  a  per  ship  basis.  However,  hot  stacks  may  be  the 
dominant  feature  in  some  ship  signatures  hence  an  average  ship  temperature  may  not 
be  meaningful  in  these  cases. 
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APPENDIX  B 


COMPUTER  PROGRAM  FOR  CALCULATION  OF  SHIP  CONTRAST  TEMPERATURES 


A  computer  program  using  the  methods  of  this  report  was  written  to  calculate 
ship  contrast  temperature  and  the  component  temperatures  making  up  its  definition. 
Component  temperatures  are  the  sky,  background  and  ship  temperatures  and  these  are 
calculated  and  presented  as  an  aid  in  understanding  t tie  probabilities  of  occurrence 
of  contrast  temperatures.  Contrast  temperatures  are  calculated  from  data  supplied 
hourly  for  day  operation,  night  operation,  total  time  operation  and  for  the  hours 
between  10  a.m.  and  3  p.m.  These  calculations  are  gathered,  statistically 
analyzed,  and  printed  in  a  form  useful  in  systems  studies.  The  definition  of  the 
probability  of  occurrence  of  contrast  temperature  and  the  rationale  for  this 
definition  was  given  in  Section  V  of  this  report. 

The  ship  contrast  temperature  program,  named  AS1RCT  for  Approximate  Ship 
Infrared  Contrast  Temperature,  will  be  described.  This  oe5>cription  will  consist 
of  a  flow  diagram,  a  listing  and  a  written  description  of  the  operation  of  the 
main  routine  and  subroutines  of  this  program.  These  are  given  to  show  how  weather 
data  is  used  and  how  calculations  are  performed  but  not  to  provide  a  user's  manual 
for  the  program.  However,  the  program  and  Instructions  for  its  3et-up  and  use  may 
be  obtained  from  the  Naval  Surface  Weapons  Center.  The  ASIRCT  program  was  written 
specifically  for  the  Patrol  Trigate  Model  and  for  the  8-12i<m  waveband.  However, 
this  program  can  be  modified  to  include  other  ships  or  wavebands  rising  the  basic 
format  of  this  program. 

The  ASIRCT  Program  consists  of  a  main  routine  ar.o  eight  subroutines.  A  flow 
diagram  for  this  program  is  given  in  Figure  (5-1)  and  the  indi  virtu:.’  routines  are 
listed  in  Figures  ( B-2)- ( B- 10} .  The  program  requires  three  punched  cards  and  a 
magnetic  tape  containing  the  weather  ship  data.  The  three  cards  arc*  initially  read 
in  to  provide  program  controls,  the  Single  Element  Method  tnermsl  model,  and  the 
ship  target  conditions.  These  are  the  only  input  data  cards  required  as  shown  by 
the  flow  diagram  in  Figure  (8-1).  The  Program  Control  Card  provides  options  over 
printing  or  plotting  the  temperatures  generated  and  the  option  of  calculating 
statistical  results  either  monthly  or  for  a  season  of  several  months.  Three 
standard  seasons  were  chosen  for  data  analysis.  These  are:  (1)  the  winter  season, 
November,  December,  January  and  February;  (2)  the  spring-fall  season,  March,  April, 
September  and  October;  and  (3)  the  summer  season.  May,  June,  July  and  August.  The 
Single  Element  Method  Card  provides  the  thermal  capacity  and  internal  temperature 
model  for  a  ship  of  three  thermal  capacities  and  three  internal  temperatures. 
Lastly,  the  Target  Conditions  Card  provides  the  speed  and  heading  of  the  target 
ship,  its  emissivity  in  the  waveband  of  interest  (8-12, .m  here),  its  solar 
absorbtance  and  its  stack  model. 

Program  ASIRCT  first  calls  subroutine  SKRAD  which  provides  the  clear  and 
cloudy  sky  radiances  in  tabular  form  for  future  use.  This  subroutine,  which  is 
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listed  in  Figure  (B-3),  contains  both  the  Midlatitude  Winter  and  Midlatitude  Summer 
Atmosphere  Models.  Midlatitude  Winter  is  arbitrarily  used  for  the  months  November- 
March  and  Midlatitude  Summer  for  the  remaining  months.  The  program  now  enters  a 
large  loop  which  begins  by  reading  the  weather  data  tape  as  seen  in  Figure  (B-l). 

The  weather  tape  provides  a  myriad  of  information,  only  a  portion  of  which  is  used 
by  the  program.  Useful  information  includes  day,  month  and  location  of  the  weather 
ship,  hour  of  day,  sea  level  conditions  of  temperature,  dewpoint,  visibility  and 
wind  velocity  and  cloud  height  and  conditions.  This  information  was  recorded  hourly 
by  Weather  Ship  J  hence  it  will  be  used  to  calculate  contrast  temperatures  hourly. 
The  large  loop  of  Program  as  I RCT  takes  an  hourly  reading  of  weather  data  and  cal¬ 
culates  the  contrast  and  component  temperatures  for  this  sample.  It  then  returns  to 
the  beginning  of  the  loop  and  repeats  the  process  for  another  hourly  reading  of 
data.  The  temperatures  calculated  are  all  saved  for  the  statistical  analysis  which 
is  performed  after  completing  the  loop.  A  description  of  large  loop  operation 
follows. 

After  reading  the  weather  data,  the  program  next  determines  if  subroutine 
SKRAO  should  be  called  to  provide  new  sky  radiances.  SKRAD  is  called  only  if  the 
date  indicates  a  change  in  season  from  midlatitude  winter  to  summer  or  vice-versa 
has  occurred.  This  happens  on  e'ther  April  1.  or  November  1  for  the  two  possible 
seasons.  The  pronram  next  calls  subroutine  CLDTR.  This  subroutine,  which  is 
listed  in  Figure  (b-4),  provides  the  cloud  conditions,  cloud  transmission  and  rain 
convective  coefficients.  Cloud  conditions  are  taken  from  the  weather  data  and 
categorized  into  7  general  conditions  for  later  printing  for  informational  purposes. 
From  cloud  type  and  altitude,  the  transmission  of  solar  energy  is  determined  using 
information  from  Reference  B-l  as  a  guide.  Rain  conditions  were  subdivided  into 
light,  medium  ana  heavy  drizzle  and  light,  medium  and  heavy  rain  for  use  with  the 
rain  convective  coefficients  derived  in  Section  II  of  this  report.  Control  of  the 
progran  returns  to  the  main  routine  which  precedes  to  call  subroutine  NEWOAY  if 
the  weather  data  indicate  the  start  of  a  new  day.  This  subroutine,  which,  is 
listed  in  Figure  (B-S)  finds  the  hour  of  sunrise  and  sunset  of  the  new  day.  The 
main  routine  then  calls  subroutine  SOLAR  which  is  listed  in  Figure  (8-6).  This 
subroutine  calculates  the  convective  cooling  coefficient  due-  to  relative  ship  motion 
which  is  given  by  Equation  (16).  Then,  if  the  sun  is  present  for  the  hour  under 
consideration,  the  indirect,  direct,  and  absorbed  solar  energies  which  are  given 
by  equations  ( 1 3 M 1 5 )  are  computed.  Control  of  the  program  returns  to  the  large 
loop  where  the  ship  temperature  including  correction  factors  are  calculated  by  the 
Single  Element  Method  given  in  Section  V  of  this  report.  Subroutine  TRSKY,  which 
is  listed  in  Figure  (B-7).  is  then  called.  This  subroutine  first  determines  the 
sky  radiance  using  Table  7  and  the  wind  speed,  cloud  height  and  cloud  cover  supplied 
by  the  weather  data.  Subroutine  TRSKY  includes  the  variation  of  sky  radiance  with 
relative  humidity  although  this  has  only  a  small  influence  as  previously  explained. 
This  was  included  before  it  was  realized  its  effect  was  small  and  is  retained  in 
case  future  analyses  depending  on  relative  humidity  are  developed.  The  sky 
radiance  found  is  related  to  sky  temperature  using  the  blackbody  relationship  given 
in  Figure  5.  Using  the  just  found  value  of  sky  radiance,  the  background  and  con¬ 
trast  temperatures  are  computed.  These  are  stored  for  future  use  in  determining 
probabilities  of  occurrence  by  statistical  analysis.  Next,  subroutine  CONHR  is 
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is  called  If  the  hour  of  day  Is  between  10  a.m.  and  3  p.m.  This  subroutine,  which 
is  listed  in  Figure  (B-8),  stores  contrast  temperatures  for  the  above  hours  when 
they  are  likely  to  be  at  their  maximum  values.  Control  of  the  program  returns  to 
the  main  routine  where  the  large  ioop  ends  by  determining  if  the  hourly  data  cal¬ 
culated  is  to  be  printed.  If  printing  is  desired,  many  quantities  in  addition  to 
the  contrast  and  component  temperatures  are  printed  for  information.  These  include 
the  indirect,  direct  and  absorbed  solar  energy,  the  wind  speed  and  direction,  the 
cooling  coefficient,  sky  radiance,  rain  coefficient,  air  and  sea  temperature  and 
the  generalized  cloud  condition.  The  program  now  either  returns  to  the  beginning 
of  the  loop  for  more  weather  data  or  leaves  the  loop  to  statistically  analyze  the 
temperatures  stored . 

Data  are  statistically  analyzed  after  1  month  or  a  specified  number  of  months  of 
contrast  temperatures,  computed  hourly,  have  been  stored.  This  is  done  by  calling 
subroutine  ANAL  whose  main  purpose  is  to  call  subroutine  SORT  which  actually  per¬ 
forms  the  statistical  calculations,  subroutine  ANAL,  which  is  listed  in  Figure  (B-9) 
calls  SORT  once  for  every  contrast  temperature  or  component  temperature  calculated 
in  the  main  routine.  These  are,  the  component  temperatures  for  ship,  sky  and 
background  and  the  contrast  temperatures  for  daytime,  nighttime,  total  time  and  the 
six  hours  between  10  a.m.  and  3  p.m,  SORT  is  also  called  to  find  probability  of 
occurrence  for  the  sum  of  the  three  correction  factors,  TCI,  TC2  and  TC3.  After 
these  calls  to  SORT,  subroutine  ANAL  prints  a  summary  table  of  the  probability  of 
occurrences  found  in  SORT.  These  tables  will  be  reproduced  and  given  in  Section 
V  of  this  report.  Subroutine  SORT,  which  is  listed  in  Figure  (B-10),  calculates 
the  probability  of  occurrence  of  temperatures  which  are  given  as  a  list  of  standard 
temperatures.  These  probabilities  are  printed  and  as  an  option  may  be  plotted 
using  a  Gould  electrostatic  plotter,  for  example.  After  the  summary  table  is 
printed,  control  of  the  program  returns  to  the  main  routine  where  it  either  ends  or 
returns  to  the  beginning  of  the  main  loop  to  read  more  weather  data. 


ANAL  SORT  ►—  PRINT  STATISTICS 


PLOT  STATISTICS 


FIGURE  B  1  FLOW  DIAGRAM  OF  ASIRCT  PROGRAM 
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PROGRAM  A5IRCT (INPUT,  TAPFt , OUTPUT) 

CeHMON/PAYH/IPCNTH,IOAY,IYFAR,ALOC,ALAT,ILCNG.nUAO,IPAGe 
ftCNHON/SOLH/TAIP,TOPT,ILCCAL,ZLAT,NOIR,MVFL,T$EA,VlSR,CFRA,ENISS, 
!01,tC'',TAS,TSC,TTRN,G0IR,0$C,C0NV,  WRFi , AZH ,CTRA NS , SOL ARA , ZENITH, RH 
(“CNHOK/PR PL/If K,KH0NTH,  AYF AR,NCPLT 

CCHHON/  STATV/Hl,M?»M3,M4,M5,Hfe,NA,NB»NC,NF,NF,NS»NT,LCN,LQ»  IN 
CCHNON/SKPR/ISFAS.ChIT.RAC, $KYT,CRAl,CRA2,CRA3,CL01,CL02,CL03 
OTMFNCTON  CRAl  (6)  ,CPA?(M  ,CRA3<6)  ,CLC1(  JC),CL02<30> ,CLOJCZO> 
DINEN^I  ON  CCOV(tO)  .TWIN  <41  ,ISPF(4)  .  ISUH  (41 . 1  NON  (121  ,KHON(12) 
OTH^NCION  Nt(*9>,*2(7  9)  ,M3<79)  »M<*  (79)  »H5  (79)  , MS <791  ,NAC79I  ,NB  (791  • 
1NCI7  9),NF  <’9>  ,NF(7  9),thoo  (3  6)  ,NS(79)  ,NT  ( 7  <*)  ,  TEMP  (791  ,  TU$T  (7  91 
CATA  P/0. C 1 745/ 

HAT A  CCOV/O., 0.1 0,0. 7 5, 0.4 0,0. 50, 0.60, ft. 75,0.90. 1.0, 1.0/ 

OATA  Tk OD/10. ,5. ,4. ,3. 6, 3. 7, 7.8, 7. 4,2. 0,1. 6, 1.6, 1.4, 1.2, 1.0, 0.9, 
14.8,P •*, 0. 6,0. *,0.4, 0.3,0. 2,0.1, 0.0.-.2, -.4,-  .6  ,*.6, **1.0, *1.4, 

1 -1. 0, -2. 7, 6, 0,-4. 0, -5. 0,-1 0.0/ 

OATS  TLI5T/60. ,5*.,50 .,45.,40.,35. ,30.,?7.,?4.,?l.,t0.,16.,15.,t4. 

1.13.,  12.,  11.,  10.,  9. 5,  9.  %*.5,  A.  0,7. 5, 7.  1,6. 5, 6. 0,5. 5, 5. 0,4. 6, 4. 2, 
13.8.3. 4, 3.  0, 7. 8,  "*•  6, 7. 4, 2. 7,  ?•  0,  i.  8, 1.6,1. 4, 1.2,  l.Q,  0.0, 0.6, 0.4, 
10.2,0.0,—  .2, —  .4, -.6, -.3,— 1.0,  — 1.2, -1.4, -1.6,  — 1.8, -2*0,  —  2.2, —2.4, 
i“?.5»-7.8,-3.0,-3.4,-Y.3,«4.2»-4.6,-5.0»-5.5»-6.0»-6.5»-7.0»-0.0» 
1-e.P, -10. 0,-1 2. «, -15.0, -70.0,-40.0/ 

PAT A  TFHD/J5. ,74., 23.  ♦ 7? . , 7 1 . , 70 . , 1 9 . , 1 5 . , 1 7 . , 16, , 15. , 14 . , 13. ,12 . , 

111.. 1,«,9.«*.  ,7.,6.,c.,4.,3.,7.,l. ,0.,-l.,-?.,-3.,-4«, —5 . , —ft . ,— 7 . , 
1-8., -Q. , -10., -11., -17. *  - 1  ’ . , -14. ,-15. ,-16., -17., -10., -19., -20., -21 

1..  3. ,-74., -?6., -?7., -? 8., -79., -30., -31., -32., -33. ,-34. 
1  •  —  35  «  »-36*»  —  J,«,-78»»— 39.  ,—40.  ,-41  .,-4?  ,  ,-43.  ,  —  **4.,  —45  .  ,  — 4ft  .  ,  —  47  .  , 
1-48. , -49. ,-50 . ,  -  e  1 .,-5 7. ,-?3./ 

PAT  4  IW IN /I, 7, 11,12/,  IS »F /3, 4, 9, 10/, ISUN/5,6.7, 5/ 

CALL  fPRSFT  (T  COUNT,  i--  ’, 

5*1 

UI*<r*r 

TCOmiT=Cl 

nTAM*ft. 

TOAGF=l 

JCA*«ft 

UAST*r 

C  CONTROLS  CVFR  CALCULATION, PRINTING  ANO  PLOTTING 

,r  AO  m,N0PBT,NOPLT,  JH0NTH,LH0NTH,ICPT,KSCAS,NN,  (KHON<  J)  ,J«1,NN1 

101  FO«HAY(7I'=,ll  141 

C  SINCLF  r^FHFNT  MFTHCO  ThFRHAL  HCOri 

PF AC  IP?, Hi, ROPY ,aCP* ?, °CPY3, PI, R7, TINT, TINT 2, TINTS, T1,T2 

102  FOR  v'A T( 4F 8.2,  2F6.7,  3F8.  2,  ?F6.  ?| 

C  SHIP  VFLOCITV, SOLAR  PPOPFRTTF*  ANO  STACK  HOOFL 


FIGURE  8  2  MAIN  ROUTINE  OF  PROGRAM  ASIRCT 
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»FAO  1C  3,  S»FFC.HCr.,FHiSSfSCL*«A»H?K,TSK,Sl 
1C3  *C**kr(T M?.?1 
»FIA?*O0. 

TFM  «. 

ZFNITH«90  . 

WSLO*le.O 

MsiPi^MSic-.e  *no.-7°» 

?73=?P-H?L01 

M*.0?*.Se*<l  .-C0S{P*?Z71  1  ••« 

A?M*RFlA?*MpG 
OC  3 7  !«!,<, 

If  (»CSrAS.GT.nfO  to  7i 
IRONtlMlWlNtll 
GC  TC  37 

31  IF  <*SFAS«  GT,?>  GC  T0  3? 
iHOsm*is°f  ( i> 

GC  T0  »7 

3?  If (<Sf4S.GT.7»G0  ’0  3fe 
IHOSI  I»  «ISUH(  T> 

33  CONTINUE 
GO  T0  37 
7fc  no  T6  T*i,AN 

36  TXONI  I)  *Kf*CN|  I) 

77  ooimt  300  t«t*Fro,M«r,f  TfHtTK,  Rf  LAZ 

?QC  FCRM AT  I  1H1 » •  TARGET  S°f  F  0  •  *F  *? .1 1  •  *NCTS  AT  HF  AO  ING»  •»  f  5  •  0  *  *OEG  FPC 
1H  NORTH  3M I P  FLFHFN*  VA«IA«LF$  7r  MTM*  •  ,F  $  •  0  *  •  A  7  IHUTM  FRCR  BOX* 
l*,F*.0,*OfC»FF5*l 

PRINT  ?01 tTINT ,HT,RCPKtfHISS,SOLAP A 

201  FORMATflM  ,•  INTtonaL  TP HF* • f r«j, i , *0EG  C  INTERNAL  M.T.COEFF **,F5. 
1  7t»MATT$/ w-C  Rro»s*,r6. 3 ,• M ATT- HP /H?*C  SURFACE  FMISSIVITT** 

1  »f X.  ?  »•  SOLAR  A^SORPT  ANCF**,f «..?) 

00  3  T* 1 i 7P 
Him  *c 
H?IIh5 
H7|l) *0 
NUT»*P 

Hsm  ’C 
Nft<n  *0 
NHIMt 

Npm  *o 

Nf (t»*P 
NS <11 *3 
nt  m  «o 
3  NCCT)-»0 


<nf'  •' 


T.",iW.'ljjPauU4im<MWy5.*!5 
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icf's-? 

TF(  JMCNTh  .OF.*  .A^O.  J‘40»!*M.LF.l  t>  ISFASsl 

CALL  ?KRA0 

J<FAS*TSF#< 

DC  u  :=i,3* 

4  NF  1 1  *  *«* 

LCN-0 

LF** 

lc«o 

«  pFA"(l,11  *UH, ALCP, ISFA,ITlMf , IOU A  0 , ILAT , TLONG, I  YEAR , I MONTH , IOAY , 
ITCH’’ ,  rtOCAU,  IMsrO,  AWIMDC  ,Vlf , IWFAT,AIRPRO,TAIR,TOPT,TSEA,ICOV»ICL 
if ,  ilc  ♦:  *r ,  i  hc  ,  cm  :  * ,  *a  vf  m  ,  c  hfll  h,  a  junk 

9  POR**ftT(lfe,lXf  A?,M2, !«.,!*,  IS ,4 II,  IX,  I2,Fi..0,iX,FF.2flX*I2*lX,F6.1, 

1  YF5.1 ,1*,?T?,  31?vF6.C*?tl«»‘r*.0l  ,l*,Ain» 

IFrCO'NT  .FO.  0»r-0  TO  10 
ICOUNTr p 
II I  NFS'*  7 
OC  tP  c 
10  CCNTIM«F 

IF (F  Oc I  1) ) 190,7 

7  IF (THONTH.FO.l^ONTMinO  TO  1R0 
IFdHCNTH.NF.  ICCN(tfn  CO  TO  ? 

6  if p*eu ,nf.o» > c.c  Tf)  11 
HCL^slLOCAi-l 7 
CTA'i=i.f) 

OC  Tn  « 

U  AVF40:TVF ao,i 9C0 

TFFfiOr'* 

IFt!HC»lTM.GF.5.8NO.TMCN,,H.l.F.10>  I  jF  A?*l 

I t (  ISFtO.FO.ISFASlGO  T0  17 

call  r»PAr 

JSF  A o  =  T  f  r 
t?  ala^mlat 
TLAT» TLAt 
MVFL=  AWTHCS 
wriD,  d*IwlMpo 
*COV  = ICOV*l 
CFPIsCCOV  ftfCOVl 
V!S3  =  Vi S 

CALL  CLnTP(lLC,T^C» imp, I«FAT,CTRANS,CONO,RCOFF) 

IllNF=ILIKr*l 
IF(TOAV.fo.jOAV»GO  tc  *»o 

OSUN^O. 

f'SO"*!).  ? 


B-7 


?\  MCUP«IL0CH-1? 

LCN*LCN*1 
DT  AH* MO UP -HOLD 
MCIO* WOUP 

IF<r,T»U.lT.«.  )OTflU*^T»UM?.C 
IFC^TAU.EC.O.  »  GO  T0  * 

CALL  «OUP<HOUP,t<;oi,?PFFC.HDG,SUNB,SN$FT,QSUN,OTAU,OTIN) 

QSOl*  (nSUN*QSOC>*OTTM/(?.0*PrPX> 

O^OOsO^UN 
WINT*WI •OtlU/»Cp» 
m«A!N  =  BCOFF*OTAI)/PC®* 

MCOL*  CONV*CTA'J/PCpX 
T  T  *  T  I  N 

TT»(TIN*QFCL*mCCL*TATP*mBAIN*TAIS*MINt*TINT)  / t 1 . ♦HCOl *HB A IN*NlNT » 
»l*nT4U/oro* 

A?*OT AU/PC®*? 

A**fPAU/RC»X3 

F»  (QS'J^nSCm  /■>.  ♦  (CONV*®COFF>  •TAIP*hI*t  I  Nt 

G=C^NV*«CCFF*NI 

M<5T<=w«*if*0T»U/cro» 

TCI = Pi*  <A?-A1 »*(P-G*T!>/<* .♦G* C A 1 » A  2 > ♦ Aj • A ?• G*Gi *B3* t A  3- A1  > •  IF-G* 
1TI) /  tt. ♦G*(A1 *A*» ♦Al*A1*G*G) 

«P  =  mikt  /  { j .  ♦hcol  ♦MffAiKmXK*| 

TC2=Ti*88*  (T!Kr>-TX»jT)  ♦T?*ee*  (TINT  3-T INT) 

TC3*S1*mSTK*<  »Hr.CL*HPATN*i  .)*  tTSK-TAIF>  -QSOL» 

1/ ( IMCOl ♦HP»IS*H^TK) • CMC cl* MBA  in*! . H 
TC=rCl*TC?*TC‘* 

TTbTT  *TC 

TTsf mtjs* TT*( i ,-fmisSI •  T  A  I  ® 

T IN*  TT 

CC  1=1  ,*f 
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T©on=ri  *S»(VT4  n  ,-FD  »TSTA 
OC  1  =  1, 'R 

J  =  I 

If  (toGO.GT.TI  ISTtinpO  r0  'n 
?6  CCNMf.W 

H?(  II  =  N$<  JIM 
ICON*  TT.trgo 

OP  t,  *>  1  =  1  , TO 

J  =  T 

TP(TCCM.r.T.TL  ISTfT,|GO  T0  41 

4**  rcNTlNUF 

4  ?  NT  (  l»  rST  {  J) 

Tf<HOUc>.lT.Sl'Ne.?P.MOU<>.G’r.SNSFT  IGC  TC  TO 

ir«n*< 

re  44  i*i, 

J“  I 

IF(Tpo*  .'••T.TLTMmiGfl  re  45 
44  CCNTIV’F 
4*  N  C  (  J  )  =»  Ct  JIM 

IF(M0UR.L,’.-'».Co.M0UO.GT.,)f,0  TO  47 
CALI  '"PNMC  IJ,  KC,MOU»t  M1 
GC  T  n  4  7 
IT  LN*IN*1 

DC  ic  »*lt»q 
J-I 

if roo  .r.T.ri,  i?t (inno  *o  if 
15  c  0  N  T  T  Nt 1  c 

15  NC  (  j)  Tf.C(  J»  *1 

CnwycPT  PA^IAT-CF  cUNCTICN  tp  PACT  ANCF 

47  OAO=%4'»OA0 

?f  (MOOOT  ,  f  ■)«  1 )  GO  T0  4* 

TF(MPCs.LT,s»iNC,oo.Mru®.GT.SNSPT»r,0  TC  46 
I  F  ( mom©  .L  f  ,0.  )  mu©*M0U©M  ?.  0 

C  c  J  »j  T  **!•»  ,^OU»,TCf  ,TAS,  TSC,TTeN,OOIR,OSC,OSUN,NVn.,NOIR,MQFt,CONV, 
KFNO,  PCOPF.TAIo^SFA,  iT,TC,PAO,SKYT,TRGD,TCON 
?I0  ©CP*4  A T 1 1H  ,P4.0f4C6.i,P7.?,?F7.?,r4.0,F6.0,F5.0,F7.2,lX,#6,F6.0, 

t if6.?.r5. ?,F7.?,1Ffe.7) 

OC  TO 

46  IF  .LF.C.  >HCI’©«MOU®*1?.0 

°  u  I N  T  -»n  ,  mou©  ,4vrl,wOIP,  kPFt,CONV,CONO,RCOFPfTAIRtTSEA,TT,TC,RAD, 
IFKYT,  ">«r,n,TCCN 

?11  CC©M AT  f 1H  ,F4.%4?X,©t. 5,r6.C ,F5.0,F7.?, iX,A8,F6.3, 3F6.2,F5.2, 

lc7. Vf,.  ?> 

48  If  <!NC*'TM  ,FC.  J*CNTM)r,0  TO  5 


M 
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jfcQMu::  IHONTM 
K  *K  ♦  1 

IFlIOP’.CC.e)  so  TO  ? 

00  CCNT  I  Nl'F 

KHO^TWiTHCS^fH-l 

f»U  AN  At  tHOO,mT,9eP¥tT*G0,TlIST,TFHP> 

iFcitftST.Fo.nco  To  199 
cc  r0  f 

199  Il*cT=l 
r.C  TO  90 
199  CCNT I suF 
STD® 
cnO 


*W> 


■UM 
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SUP9PUT I Kf 

common/ cKPP/i^ra^^HiT^oao^KYTtC^Ai  ,C9A?,CRA3,cirn  ,clo?  ,ric3 
o 7  nr  ns  ton  c°ai  f8»  ,C9A?<fe)  ,r«?A3  <f,»  ,rtoi  no*  301  <3  01 

nTMr  n,c  |pi.  C5UM 1  (6),C$UM?(ftl,CSU'43<6>  ,  r.  W I  N 1  (61«rH!N?(ft|  ,CWTNM61* 
ir9SMi  i  *o  >  ,cns»2  ( 301  .c.-ismac*)  ,  roMM  no»,  cown?(  30  > ,  count i 30  > 
nata  c$‘jni/40.6516,*,i.7i  68,43.7315,47.4141  .59.7P9°»78.tf>oo/ 

"AT  a  CSJN7/40.  6?1  3,  <*1.  1  SI  P.4 '.0577, 4  7.  09  70  ,58.  87??, -*8.  06  7*7 
OAT  A  :SUWV4Q.  58  6  A  «4t.O  7  0S,  4?.  <*6 7  4, 4  6.  5 4 70  ,9  7 .  a06  0 ,  7  7.  90  37/ 
rU^Ann^l  /6A.  *66Q»ft6.  9?04«67.1  301  »67.715C,7n .1 194 , 77.565  9, A9.461 e, 
l6',.rC?,»^<).f7i  7,70.0096,  71.6961,  77»5484»  ?7  » 1  ft9A»7?.1709»7‘,.7379» 
t  ?’..*»c4, 73.  30  3?,7  7.4  79?,Ts..q49fc,7<,.956fl.  74.9  766,  rr  .  037 1»  75.7176  * 
1"»‘',n9‘;4t79,9ir'>0»7*'.9Mont74.«>43?»T6.9‘>?*,  76.4  9o1f7‘r.66.36/ 
r'ATlCOSM?/66.9616*66.>1^4»67.0999f6?.63?3,69.99»6»77.4  21  4t69.4593t 
l  ft  4. 4Q47,  f  Q.60  * 6, ft e. 9480  ,71 . 7941, 7 -.4  17fe,  7’ .  I  4  A  0  *  7?  .  1  f  64,  7?,  7  747, 

17  ?.  *26b,  ?’.16r,4»7?.4474,74,94o0,74.<>5'.9,74.O71Tt7p.0?11»75,?ft4Ot 

t  43*4, 7ft.  9  3  8  9 ,  76.  9’  -)7  *  7fc  •  °4?4,  7  6.9604,  ?A.f*89A, 7’,5cie*4/ 

nftTSCO"wV6*.*,56',»66.o1*7,67.06'f4»67.6  4lt:,69.95A*>,7*.79M4,6s.4569, 

16-9.4  871  ,<9.c  86  0,69.8*9 1,7 1,1  690  »7  7.7506,  77.1461,7  7,1619,  77.' 10  7, 

1-»;>.<568,?3.0?75.77,?0''7,74.94  86,74.9  531.74.9  67l,75.8tOO,7K.713», 

17  ’.^4?,  7A.9  3M*,7e.ot-»6.-ft.0  4l7,  'ft.'UA7,  76 .0  •  0 9,  7  7 . 57u  6/ 
piTC  T f  l/?1.957?,’7.OT~4,,f.6491 ,33. 84ft*. 53. ’8.1801/ 

PlTf,  rwTt>7/?l. A6AP,?7. 80  Ou,  76.  09  30,  3’.8??c,5?,OO7f,7*.0i»ft37 
OATA  CH!»  3/71. 7841,  ->?.6M  %76.6RP»,M.8'»°3, 60. 0894, 77.80  *4/ 

oaTar''HNt/6  4.4360.64.t:7ft0»64.*16T,6r-.6  40  4,6A«O‘»4?,7'f,5  9?9,6,.99?4, 
16*.0  449,f.A,?^l9,fA.77  49,70.»041,’7.6,)6,71.4  3?5t71.4  6O9.?1.t441, 
t,l.'rQ70,77.*631,7?.  5975,  ’4,  7ifti,7c,.T?7ft,  76.7478,74.3193,  ’5. 1 470» 
17^.  •" ’65,  76. A  94  9. 76.*  940, 7ft. A  Q9  fc, Tft.oiO  3,  7ft. 0617, 77. 8  79  9/ 
'V\Taft''WV?/64.4?6Q»A4.ft067,A4.7ft’9,ft5.**706»ft3.405A»77.388**t*7.988,» 
16fl. 03?ft, 6". 17*3,ftA. 6043,70. 4175, 7’.3A?0, 71. 4791, 71. 4508, 71. 5193, 

1T  t. ■*  7 ’1,77.6716, 77. 4136, 7^. 7160, 74.7?10,  74. 7  40  8, 74. 7965, 7^.0605, 

1’  ‘.  ^Of-4.  7ft.  *94  7,76.  8  066, 76.  8 08 4, 7 6. 90ft 8*  74.9  474,77.5776/ 

04Ta  C  PW  *' 3/ftU. 4185.64.. *  6C  »r>4. 6  96 1. 66. »1«»6»67.919*»77. 1674,67.0  *90, 

lft,.n7',ft,f*. 1417, ft*. 5039, '’n«0734,77.l70  1,71  .476?, 71. 4441,  71.4984, 

1  ’  1. ft ?  31  ,  77.4  19 7, ’7. ’O'*.,  ’4.  714  3,  74. 7  1 91, 74. 7 7  39, 74. 77*8,  74.6970, 

1 7 ?4  3, 76.  *94  6,76. *036, 76, 8 97ft, 7ft . 9 043 , 7ft .9 *73 . 77 . 4079/ 

I p  (166  43 .50 . ?»  GO  TO  71 
00  10  T : 1 ,6 
r»at  ( T ) =rftu“t ( T > 
m’  ci  irsu^’tTi 
10  r^’in  =rft!|M3(ii 


FIGURE  8-3  SUBROUTINE  SKRAOOF  PROGRAM  ASIRCT 
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on  i *i  i  =  i*?o 
olci <i>  *rns*M m 

no?  (TI=C0SM?(l> 

15  ci.ni m  =ccsh?ci» 
*.r  tt  40 
?o  pn  io  i*i*6 

r=>M  C!)=CWIN1  f  T  * 
CCfi?(T»*fWlN?(TI 
10  <!l  =fWlNlf  II 

no  ?5  i  =  t.m 

cmi  (  T  1  =COMNl  (II 

Cio? (Ti srown?(i» 

n  m (Ti eroMNi ( ii 
40  rO^'T  J*Mjf 
or  T(|CM 
r  »ip 


X 

t 


| 


B*12 


-  7 
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SU^OUTINF  ClOTRITlO.TNC,  TMC,lWF-AT,cmNS,C0N0,RCOFF> 
"TENSION  l  ITR<15> , IMr0  (141 ,  IMFVC  1 11 ,LIT0(4>  ,JMFP(  31  .  JHEVOI 
nTM^NctlON  WF  ATM ( 7  I 

FATA  WFATM/8HN0  PRFC  ,8MPRFC  WHR, 6HF0G/TCFF, 8HOR * 77LF  ,8MRATN 
16M5N0W  .8M5M0WFRS  / 

haTA  L  I  re/60, 61,  Ff,  66,  70 , 7t , 76 , *6 ♦  60 . 6  3, 85 ,6  7, 69. 91 , 93  / 

OAT  A  TMFr/6?,6  3,7?,  73,  77,61. 64,66, 68 ,90, O’, 9  4.95,06/ 

DATA  IMF V/64,65,67, 69, *4 . ”5 , ’9 , 6? , 97 , 98. 99/ 

DATA  LIT 0/5 0,61  ,56,58/ 

DATA  JNFp/52,53,5?/ 

OATA  JHFV/64 ,55 ,5°/ 

TFUir.FR.OJG'l  TO  50 
te  (HO.GT.21Gr  TO  20 
CTO,*  KS  =  P  .6 
G"  TO  «(* 

50  T*«IlC.Gt.3ir,P  TO  25 

C Ts A N5*  0  .c0 

r. 0  rO  60 

25  ir  (UC.GT.5l  O'*  TO  30 
CT  PA  6*5*0 ,33 
G"  rr  60 

TO  IF(TLr.FT.7>5G  TO  .35 
rTPAN9aP.?5 
GO  TO  60 

35  IF  (TlG.r-T.8ir."  TO  40 
OTPA6S*0.? 
r,0  TO  6  f* 

4  0  rrot.ns-0.? 

GO  TC  80 

50  Tr  (T^O.FQ.O) GA  TO  60 
IFf  IMC.GT.51GO  TO  55 

C  T  o  A  ‘!Gs  n  .  41 
GO  TO  60 
5*  rTOAN9*0.49 
GO  tc  60 

60  Tc (TMO.CT.OIGO  TO  65 
CT»ANSM  .00 
r,o  TO  60 

69  IF (IMC.GT.7IGC  TO  70 
(>5PANS*0.  *3 
GO  TC  *0 


FIGURE  8-4  SUBROUTINE  CLDTR  OF  PROGRAM  ASIRCT 
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•»0  IP (IHC.GT.8) GO  TO  75 
CTPANS=0.74 
GO  TO  50 
▼5  CTRANS*0.65 
40  OONTINUF 

ic  (TWFAT.GT.19)  GO  TO  85 
CONO*WFATW(l ) 

GO  TO  1?0 

*5  IP(IWPAT.GT.39»C0  TO  90 
CO’IOsWF  A  TM<  ? 1 

on  to  i?o 

0(5  TP  ( IMF  AT  .GT.491G0  TO  05 

cono=wfathm» 

GO  TO  1*0 

op  TF<IWFAT,r-T.59>  GO  TO  105 
C0W0=WPATM(4) 

G r>  TO  170 

105  IF(JWFAT.GT.5OIG0  TO  110 
CONO-MF  ATM(5> 

GO  TO  120 

110  IPIIWFAT.GT. 79»GO  TO  115 
CONO=WF  ATWf 51 
GO  TC  120 
115  CONO=WF  A  TM( 71 
1?0  OC0PP=0. 

TP <IWPAT,t7.501 GO  TO  150 
IP  (IWFAT.GT.50I  GO  TO  140 
no  1?1  T  =  l«4 

TP(TWPAT.FQ.LIT0(T» IGO  TO  12? 

;n  contiwup 

GO  TO  123 
122  prOFF«l?. 

GO  TO  1*0 
1*3  00  124  I  =  1 1  3 

TP(IWPAT.FQ.JMFO(inGO  TO  125 
1?U  oontinuf 
r.o  to  l?f 
125  PCOPPs?0. 

GO  TO  150 
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1?6  DO  1?T  1'1*3 

IMIWPAT  ,F0. JWFVfl)  >60 
\?7  CONTTNUF 
i?*  ®rorF=3o. 
r,n  to  150 

130  on  13*;  j*i»15 

TF(THrar.FO.UTO(J))GO 

l 35  rONT’NUF 
GO  TO  139 
1  36  prOFF=»0. 

GO  *0  150 
130  00  140  J*t*t4 

IC(IHFAT.Fa.THfCU>>GO 


TO  1?8 


TO  136 


TO  14? 


140  O0NTTMUF 
m  TO  144 
14?  pr,orf  =  i35* 

GO  TO  15P 

144  00  145  J-l*ll  ^ 

TPf IWFftT.FO.TMFVCJI >00  TO  14* 

145  C0NT7HUr 

1 48  PFOFFsPOO. 

150  roMTi*iuF 

prTUPM 


F  NO 
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SUgcOOTlN-  NfKOA Y< ILINE  ♦  7S0L  ,SUNR,SNSfT1 

COMCN/OAYH/InnNTM,IOAY,  IVEAR,  ALOC,  AUT  »  ILONG,  IQUAO,  IPAGf 
C0eM0N/P9Pl/KK ,KHONTH,AVEAR,  NOPLT 
ntMPNSION  Of  CL  <•♦  91 
CATA  P/C. 01*4$/ 

OAT  A  0ECL'-1388.  ,-l34C.  ,-1275.  ,-1193.*-i036*»-814.»-775.»  -62  7.  *-47 
11.  ,-’l  3., -145., 21., 256.. 41 6.  ,5 70 . ,7 1 7. , 891. , 1013. ,1121., 121*. ,1317 
l.« 1367. • 1*97. • 14 C7.« 134:. ,1354. v 1299. .1225. , 1 C 92. ,981 . ,857 . , 722. , 
191  3.  ,758.. 199.  .3  6.  ,-l7«*.  .-?3  6.  ,-495 . ,-648. . -852. ,-932  .*-1098., 

1-1  199., -l?7i*., -t  358  .,-1  3  94.,  -1417., -138  *./ 

II=4*(lHONTH-l» 

I  J=  ( ID  AY  -1 »  /  7 
IK=II*IJ*1 

If  CIOAY.Gf  ,?°1  IK*IK-1 
PR  AC-IOAY-7MJ-1 

7$OL=OSCL<IKU  lOtCLClK^U -Of  CL  CIK1I  »  (fRACI/7. 
t50L  =  ZS0L/8C  . 

7l=SPl|f  *4LAT1 
77=rcS (P*ALATI 

«*-SIN  (P*7SJU*SIN  lP*ALA7»/<CCStP*7SOL»  *C0S  IP* ALATI 1 
Aa*ACOS<  A>  /p 
S*JN"  =  -A0/15.  *1  2. 

^NSr  T=  ?*4  .  -SUN’S 

jlat=alat 

If  (NCPLT.f  0.  II  GO  TO  20 

°R  I  NT  ILC.  IfiONTM  ,IOAY,IYf  AP.  ALOC,JL  AT,lLONG,IQUAO,SUNR,SK‘3Er 
150  cORMAT  <1M]  f*  0  AT  f  IS*, 12,*/*, 12,*/*,  I*,,*  HfATHfR  SWIP*,A2,3X, 

1  *  A  T  LATIT'IQp  • ,1 2, 3X  ,*ANP  LO  NG  IT  UOF  *  ,  I  *  ,  *  IQ  UA  GRANT*  ,  1 2  ,*  )  SUN  RIS 
If.  AT*, F5.2, 'HOURS  SUNSFT  A  T*  ,f  5 .2  ,  *HOURS*1 

If  <ILINP.GT.371IPAGF*1 
If  <IPAGE.fO.OGO  TO  20 
P°INT  215 

215  fORMATIlMl,*  CLOUD  AfRO  SCAT  TOTAL  OIR.  INOIR  SOLAR  WlNC 
1 01  Of CT  Rf L  CONVT  WEATHER  RAIN  AIR  T  SEA  T  SHIP  CORR.  SKY  8  SKY 
1  T  ° ACKT  C 0 N T R • I 
PRINT  22  2 

22  0  f C°NAT  <1H  ,*  TIM  TRANS  TCANS  A0SOR  T PANS  SOLAR  SOLAR  A»SOR  SPEED 
ICG-NTH  mind  COEff  COOE  COFfF  Off,  C  OEG  C  CMP  DEG  C  f  AD  12  OEG 
1  C  DEC-  C  OEG  C*> 

IUNfsO 
I°AGf =  0 
20  CONTINUE 

SNSET=SNSCT-12. 

SUNOaSUNR-1? . 

Rf TURN 
f  N  C 

FIGURE  B-6  SUBROUTINE  NEWDAY  OF  PROGRAM  ASIRCT 


All 
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5U9R0UTINP  SOL  AR (HOUR. ZSOL  »S  Pf  EQ  «HDG  »SUNR»  $NSET»  QSUNt  OTAU  »  OT IMI 
CO HMON/SOLH/TAIRtTO°T,T LOCAL  »ZLAT »HOIR«  WVEL ♦ T5EA* VISR »CFR A*ENISS » 
101  ,TCC«  TAS,TSC*TTRN,Of  !8  ,OSC  ,CON  V  ,HREL  »  A  7*»  CTRANS*  SOL  ARA.  ZENIf  M,  fin 
OAT  A  P/C.J174,*;/ 

AAa?73.l5/  1273.1  *5*1  AIR* 

HS  AT  =  AA*E<P(  13.9  76  6-1**.  9595*  AA-2  .4! j ft*AA*AA> 

AA=?7*. 15/(273. 15*T0PTl 

4C0N=AA*f  XP<  13,9 766-1 A.  9595*  AA-2  .4788*AA*AAI 


9HS10C .•W3GN/WSAT 
MUHIO=HSAT*9H/10  CO  . 

72  =  ?rN(P*7SCD  CALT02 

74  =  COS (3*/S0L1  CALT32 

71 =SIN(P*7LAT> 

73  =  COS  |P*  7L  AT) 

75  =  COS  <°*15 .  *MOU«» 

76  =  7 1*  72*  7  3*  74  *7 5  CALT03 

SOL 7  =  ( AC’S  1  76 )l /*  CALT03 

HI  =  h  VC  L*SI N<°*  W01RI  CALT10 

H2  =  HVf  L*COS(P*HOIR1  CALT10 

W3  =  ?  PF  FO*S IN (°  *MQO) 1  CALT1C 

=  SP^F0*CO3 (P*hOG»*H?  CALT 11 

HPPL  =  SQPMw.3**?  ♦  H4**?1  CALT11 

C0NV=3.6P?5*K*FL  •*  3  «  6 
IF  (CONV.  LT.bICCN  V=fe. 


C  CO  NV  ?  b  T  WIND  SP^tOS  FCOH  KNOT*  TO  ^STtRS/SEC 
HVSL  =  hvEl/2.  C 
H?fL=w*tL/2.  C 

?T  IM  =  CTA'J 

Ir  (S0L7.LT. «Q.»  ?:,IC3  CALT03 

33  Ylsr0S(P*7LAT)*SINIP*7S0L> 
ve  =COS  (p,7sou  •SIN(P*7LAT» 

Y»=vi-Y?*''5 
*=~OS(P* (9C.-S0L  71  I 
PJ  S0LA7=  ACOS  (73/XI /P 

IF  (WOUR.C  *.:•.»  SOLA  Z=  363  .  -SCL  A7 


C  FOLA®  HFST  iciq  CALT9<t 

C  CALCULATION  CF  StANT  path  EQUIVALENT  tl°  “ASS  CALTOV 

IP  <  S0L7.LT. 6C  .»  5 1  ♦  5?  CALTO* 

c  1  A  =  l./CO?  (°*SOLZ»  CALT  o* 

A1  =  A  CALT05 

00  TO  53  CALT05 

52  A  s  l, /COS <P*( 30L7  -  .011  CALTO? 

61  =  i  ./COSCP*  (SOL  7-.-.1  »  CALT05 


FIGURE  B-8  SUBROUTINE  SOLAR  Or  PROGRAM  A8IRCT 
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53  CONTjsUf  CALT05 

mw«1.j*mu1I0  CALT95 

TM5*K*»C  •CFXP|-.i93*A>-l.M  CALT05 

A3 w*FXtM -2.5/1  IWW*  Alt  ••.  D975I1  CALT05 

TSC*TMS-A«W 

TA<  s  -:«P<  -5.  •U/VISRI  CALT05 

T<:  =  1  UTms*TAS»/2.»  CALTQ5 

77  =  5IU(P*5ot?l  CALTtte 

e  *»7t  ••CFPA1*  (1  .-CFPA»*<>1  CALT06 

T55  =  1  .-<1  .-R>»<  l.-CTKANS*  CALT05 


TTRN  =  TCC*TA5*TSC 
CO  I  S*1  35  j .  *T  TftN 
GcC=Qn  fi<*75 

5 1  OX  *COS  <*»•<  A7H-S0LA7}  » 
l)’s3IMP,7IFSirHl 
ni  *005  (0*7FNITM| 

n*.  =  7  7*  *0?*  C3  •  7b  CALTC8 

IF  UK.UT.O.I  3fc*C.  CALTC8 

A2  =  T  (A  7^ -SOL  A  71 

ns  s  .  >5M1.  ♦  2.*77*UC0S«n*aC»  CALT09 

lr <  o  sun. ct . o. 150  ro  i? 

0  T  I  msm  O’J  3  »S'JKr 

*?  CSU7'sT-0lA5  A*  (L  *•  cni»»r5  'O^CI 
cr  T  1 1  *  ►> 

SUN  RTLOW  m05  I  73N  CALT09 

133  T  C  C  *  3 . 

TS  5*;. 

Tc r=c. 

T” =  C  . 
r-)io«;  . 

or>c=c. 

if  cnsuN.po,  i.»  on  to  tu 
CT  I^2M  3M-?-SN«£  t 
ovjn*:  .3 

no  roNtp.ijc 
P?  TU°N 
FNO 
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SU«9CUTINF  TPSKt  (OFRA.Rh.wS.TAIR  »TpKP»NRI* 

Cf'RHON/0«oF/lSFAS,CMlT,KA0,SXYT,CaAl,CRA2,CRAJ*CL0t«et.02,Cl.03 
01 Hr NO  ION  09 A1  (6  I «  CRA?  <  6 ) .  C9  A  X tb) ,CL01< 3Q> *CLD?UO» ♦CtOSf 33»  ♦ 
ISOHUt  ,ACmT<5>  ,AWSI6».RAOL((>».9AOUt6I.RORLCfe)  *RCRU<6» 
r-  hfksiqh  R6«  0  17  9*  «  Tf  K°  C  79 1 ,  CC  OY  1 ID  >  .NR  <79> 

DATA  u  CO  V/  0.  «  .  1*  t  25»  «5*.b*.75*»9»l.O  •  1  •  0  / 

OATA  A?m/10:  ..35.,b0./ 

CATm  ACMT/17!>U.»  i?  sc.  ♦  *3  <).•<•?:  •  •  1  ./ 

C  A  T  A  A  m  0  /  'i  «  0  « 5 « C  .lC.j.iS.C.^D.O.  T  C  *  0  / 

cat  A  n^o/m.  f,fc76,9  9.  3  8  75.9  8.  32  N6*  9  f .  b  794 , 9*> .  Q  512  ♦  93  .  <*hC  2  « 9 1 .  H  54 
1  .90.25  9?  .  -4ft.  71  CC  ,* 7.157  *•  A5.  b-*l7»6W.  H?9.3?.5Ml*8l  .1660*  79. 7978 ♦ 

1  78. ’f  62.  ’S.AM*.  75.4  <31 . 74.0  415,  72.5  55  3 ,7l.3C  7S.59.9653.68.6393* 

16  7.  *29  6,  SS.OUl*  64.  75M.  63. 4  975,  S'*.  2521. 61.  0230 , 59. 80 97. 5 6, 6 i2t ♦ 
lc7.*»J:-3,  56. 2 b*2,  55.1137,53.9787,  52.  3592*51. 755 1 , 5C .  6663.  «•  9 •  5927 * 
148.51UT,.47...9l0.46.4f?r  ,.,5.4Si9?,*4.-*50  7,41.4663*42.4977.4i.5431* 
ii*C  .  6  13  D,  H.6V3, 38.  2218*  M.7o5<»,  37.  9688*  86.9857.36.1168*  35.261?* 

1  ?4.  4  ::  5,  !«.  59  51,  72.7  743*  Jl.  9*oi  , »i  . 19? 3 , 30 . 41 69. 29. 6588 , 28. 91 18  . 
t?M7M,  ■>?.  +  5k9,  ’6.  C  5  3*.  25. 3  7J  6.2-..  7  0  24, 24.  C426.  ?  3.  3971. 2  2 . 7637. 
122  .  1h?  >,  n  .5  33?,  ?3. 9  36  A  ,  ?C.3  5C  1*  19.  ’75  l  *  1  9.  21 3  6*  18 . S6?4*  l  3. 1  226/ 

to  ir  T  =  i « 5 

1 1  =  T 

IF  <C«I  T.  '“•c.ACmT<  I>  150  TO  15 
ID  CONTINUE 
IS  T  a  1  *  *>•  (I  I  -  II 
U  =  I*6 

:c <»w.LT.^5. ICO  TO  ? 5 
Tr  (9M.  t^.351  r.o  TC  30 
IF (Pw.FQ.iCC .100  TO  35 
J:  0 

T  ?C  N= I . I 1 

J*  J«1 
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